Transcriptomic analysis of yak longissimus dorsi muscle identifies genes associated with tenderness.

Meat tenderness is an important sensory index when consumers choose meat products, which determines the value of meat products and consumers' buying intentions. Yak meat is rich in nutrition and unique in flavor, which is favored by consumers. However, its meat has the deficiencies of low tenderness and poor taste, which has a negative impact on the value of its meat products and customer satisfaction. To identify the genes affecting the yak meat tenderness, we used RNA-seq to analyze the longissimus dorsi muscle of yaks with different tenderness, screened a total of 1120 differentially expressed genes (DEGs). Meanwhile, 23 pathways were significantly enriched. By further analysis, we identified eight genes related to yak meat tenderness (WNT5A, ARID5B, SERPINE1 KLHL40, RUNX1, MAFF, RFX7 and ARID5A). Notably, SERPINE1 was involved in the significant enrichment pathways of 'complement and coagulation cascade pathway', 'HIF-1 signaling pathway' and 'AGE-RAGE signaling pathway in diabetic complications' which can regulate meat tenderness. This implies that SERPINE1 may play an important regulatory role among them. The DEGs associated with yak meat quality screened in this work will be helpful to identify potential biomarkers related to meat tenderness.

Plant community-mediated effects of grazing on plant diseases.

Grazing is one of the most important management practices for grasslands. To date, most studies on how grazing affects plant diseases have focused on a single plant species, ignoring plant community characteristics and phylogeny. We used data from a 6-year yak grazing experiment (0, 1, 2, and 3 yak(s) ha - 1 treatment) in an alpine meadow ecosystem of Qinghai-Tibetan Plateau, from which we tested grazing effects on foliar fungal diseases at both population and community levels. By measuring plant community variables (including richness, evenness, phylogenetic diversity, and composition) and disease severity, we evaluated the relative importance of plant community-mediated effects of yak grazing on community pathogen load with a multi-model inference approach. We found significant differences in pathogen load among different grazing treatments; we recorded the highest and lowest pathogen loads in the 1 yak ha - 1 treatment and in the 3 yaks ha - 1 treatment, respectively. Pielou's evenness index and community proneness (i.e., an estimate of the capacity of plant communities to support diseases) best explained variation in pathogen load, indicating that plant community-mediated effects (through evenness and proneness) of yak grazing determined pathogen load. Our study provides empirical evidence that grazing influences foliar fungal disease prevalence through plant community evenness and composition, which demonstrates the necessity of incorporating host plant community characteristics into disease load prediction frameworks.

Contrasting effects of mammal grazing on foliar fungal diseases: patterns and potential mechanisms.

Plant pathogens and their hosts often coexist with mammal grazers. However, the direction and strength of grazing effects on foliar fungal diseases can be idiosyncratic, varying among host plant species and pathogen types. We combined a 6 yr yak-grazing experiment, a clipping experiment simulating different mammal consumption patterns (leaf damage vs whole-leaf removal), and a meta-analysis of 63 comparisons to evaluate how grazing impacts foliar fungal diseases across plant growth types (grass vs forb) and pathogen life histories (biotroph vs necrotroph). In the yak-grazing experiment, grazing had no significant effect on disease severity, and grasses experienced a higher disease severity than forbs; there was a significant interaction between pathogen type and grazing. In both the yak-grazing experiment and meta-analysis, grazing decreased biotrophic pathogens (mainly rusts and powdery mildew), but did not affect necrotrophic pathogens (mainly leaf spots). The clipping experiment suggested that grazers might promote infection by necrotrophic pathogens by producing wounds on leaves, but inhibit biotrophic pathogens via leaf removal. In conclusion, our three-part approach revealed that intrinsic properties of both plants and pathogens shape patterns of disease in natural ecosystems, greatly improving our ability to predict disease severity under mammal grazing.

Transcriptome profile and unique genetic evolution of positively selected genes in yak lungs.

The yak (Bos grunniens), which is a unique bovine breed that is distributed mainly in the Qinghai-Tibetan Plateau, is considered a good model for studying plateau adaptability in mammals. The lungs are important functional organs that enable animals to adapt to their external environment. However, the genetic mechanism underlying the adaptability of yak lungs to harsh plateau environments remains unknown. To explore the unique evolutionary process and genetic mechanism of yak adaptation to plateau environments, we performed transcriptome sequencing of yak and cattle (Bos taurus) lungs using RNA-Seq technology and a subsequent comparison analysis to identify the positively selected genes in the yak. After deep sequencing, a normal transcriptome profile of yak lung that containing a total of 16,815 expressed genes was obtained, and the characteristics of yak lungs transcriptome was described by functional analysis. Furthermore, Ka/Ks comparison statistics result showed that 39 strong positively selected genes are identified from yak lungs. Further GO and KEGG analysis was conducted for the functional annotation of these genes. The results of this study provide valuable data for further explorations of the unique evolutionary process of high-altitude hypoxia adaptation in yaks in the Tibetan Plateau and the genetic mechanism at the molecular level.

Large herbivore grazing accelerates litter decomposition in terrestrial ecosystems.

Plant litter decomposition is critical for carbon and nutrient cycling globally. However, the effect of large herbivore grazing on litter decomposition and its mechanisms remain less explored. Here, 1203 paired observations and 381 independent experiments were analyzed to determine how litter decomposition and nutrient cycling respond to changes in grazing intensity. Grazing significantly increased litter decomposition rate by 14.08 % and litter carbon release by 5.03 %, and this effect was observed in grasslands and croplands but not in forests. The positive grazing effect was also found under sheep and cattle/yak grazing. Moderate grazing advanced the home-field advantage effect but inhibited under heavy grazing for grazed litters. The grazing effect was larger for high quality litter than for low quality litter. Litter decomposition slowed under >10 years heavy grazing but accelerated under moderate grazing. The effects of large herbivore grazing on litter decomposition were jointly influenced by grazing intensity, livestock type, climate condition, decomposition duration, litter quality, and soil properties. Our results demonstrated that large herbivore grazing accelerates litter decomposition globally and emphasized the significance and importance of grazing intensity on litter decomposition, which should be integrated into terrestrial ecosystem models.

Pangenome and multi-tissue gene atlas provide new insights into the domestication and highland adaptation of yaks.

BACKGROUND: The genetic diversity of yak, a key domestic animal on the Qinghai-Tibetan Plateau (QTP), is a vital resource for domestication and breeding efforts. This study presents the first yak pangenome obtained through the de novo assembly of 16 yak genomes. RESULTS: We discovered 290 Mb of nonreference sequences and 504 new genes. Our pangenome-wide presence and absence variation (PAV) analysis revealed 5,120 PAV-related genes, highlighting a wide range of variety-specific genes and genes with varying frequencies across yak populations. Principal component analysis (PCA) based on binary gene PAV data classified yaks into three new groups: wild, domestic, and Jinchuan. Moreover, we proposed a 'two-haplotype genomic hybridization model' for understanding the hybridization patterns among breeds by integrating gene frequency, heterozygosity, and gene PAV data. A gene PAV-GWAS identified a novel gene (BosGru3G009179) that may be associated with the multirib trait in Jinchuan yaks. Furthermore, an integrated transcriptome and pangenome analysis highlighted the significant differences in the expression of core genes and the mutational burden of differentially expressed genes between yaks from high and low altitudes. Transcriptome analysis across multiple species revealed that yaks have the most unique differentially expressed mRNAs and lncRNAs (between high- and low-altitude regions), especially in the heart and lungs, when comparing high- and low-altitude adaptations. CONCLUSIONS: The yak pangenome offers a comprehensive resource and new insights for functional genomic studies, supporting future biological research and breeding strategies.

Effects of Climatic Change on Phylogeography and Ecological Niche of the Endemic Herb Elymus breviaristatus on the Qinghai-Tibet Plateau.

Past climatic and topographic variations have created strong biogeographic barriers for alpine species and are key drivers of the distribution of genetic variation and population dynamics of species on the Qinghai-Tibet Plateau (QTP). Therefore, to better conserve and use germplasm resources, it is crucial to understand the distribution and differentiation of genetic variation within species. Elymus breviaristatus, an ecologically important rare grass species with strong resistance, is restricted to a limited area of the QTP. In this study, we investigated the phylogeography of E. breviaristatus using five chloroplast genes and spacer regions in natural populations distributed along the eastern QTP. We identified a total of 25 haplotypes among 216 individuals from 18 E. breviaristatus populations, which were further classified into four haplogroups based on geographical distribution and haplotype network analysis. Notably, we did not observe any signs of population expansion. High genetic diversity was exhibited at both species and population levels, with precipitation being the main limiting factor for population genetic diversity levels. Higher genetic diversity was exhibited by populations located near the Mekong-Salween Divide genetic barrier, suggesting that they may have served as a glacial refuge. The significant pattern of genetic differentiation by environmental isolation highlights the influence of heterogeneous environments on the genetic structure of E. breviaristatus populations. Additionally, the results of ecological niche models indicated that the geographic distribution of E. breviaristatus populations has decreased rapidly since the Last Glacial Maximum but is not threatened by future global warming.

Moderate Grazing Promotes Arthropod Species Diversity in an Alpine Meadow.

Livestock grazing is an important tool used in grassland land management practices. Studies have substantially addressed the effect of grazing on plant species diversity, revealing that moderate grazing increases plant species diversity. However, few studies have dealt with the relationship between grazing and arthropod species diversity, which remains unclear. Here, we hypothesize that moderate grazing promotes arthropod species diversity because arthropods are directly or indirectly dependent on plant diversity. In this study, we conducted a two-year plant and arthropod survey from 2020 to 2021 at four levels of grazing intensity, i.e., nongrazing (as a control), light grazing, moderate grazing, and heavy grazing, of the long-term grazing experiment starting in 2016. The data show that plant species diversity peaked in the moderate grazing treatment, and herbivore species diversity was positively correlated with plant species diversity (and hence peaked in the moderate grazing treatment). Moderate grazing promoted parasitoid species diversity, which was positively correlated with herbivore species diversity. However, predator species diversity did not significantly differ among the four treatments. In addition, saprophage species diversity decreased, whereas coprophages increased with increasing grazing levels, such that species richness (but not species diversity of detritivores statistically) was highest in the moderate grazing treatment. Consequently, the species diversity of arthropods as a whole peaked at the moderate grazing level, a phenomenology that is consistent with the intermediate disturbance hypothesis. Considering that moderate grazing has been found to increase plant species diversity, facilitate soil carbon accumulation, and prevent soil erosion, we suggest that moderate grazing would maximize multi-functional ecosystem services.

Maternal phylogeny of a newly-found yak population in china.

The Jinchuan yak is a new yak population identified in Sichuan, China. This population has a special anatomical characteristic: an additional pair of ribs compared with other yak breeds. The genetic structure of this population is unknown. In the present study, we investigated the maternal phylogeny of this special yak population using the mitochondrial DNA variation. A total of 23 Jinchuan yaks were sequenced for a 823-bp fragment of D-loop control region and three individuals were sequenced for the whole mtDNA genome with a length of 16,371-bp. To compare with the data from other yaks, we extracted sequence data from Genebank, including D-loop of 398 yaks (from 12 breeds) and 55 wild yaks, and whole mitochondrial genomes of 53 yaks (from 12 breeds) and 21 wild yaks. A total of 127 haplotypes were defined, based on the D-loop data. Thirteen haplotypes were defined from 23 mtDNA D-loop sequences of Jinchuan yaks, six of which were shared only by Jinchuan, and one was shared by Jinchuan and wild yaks. The Jinquan yaks were found to carry clades A and B from lineage I and clade C of lineage II, respectively. It was also suggested that the Jinchuan population has no distinct different phylogenetic relationship in maternal inheritance with other breeds of yak. The highly haplotype diversity of the Pali breed, Jinchuan population, Maiwa breed and Jiulong breed suggested that the yak was first domesticated from wild yaks in the middle Himalayan region and the northern Hengduan Mountains. The special anatomic characteristic that we found in the Jinchuan population needs further studies based on nuclear data.

Population genome of the newly discovered Jinchuan yak to understand its adaptive evolution in extreme environments and generation mechanism of the multirib trait.

The adaptation and diversity of animals to the extreme environments of the Qinghai-Tibet Plateau (QTP) are typical materials to study adaptive evolution. The recently discovered Jinchuan yak population has many individuals with multiple ribs. However, little is known about this yak's origin, evolution, and the genetic mechanisms that formed its unique multirib trait. Here, we report a valuable population genome resource of the Jinchuan yak by resequencing the whole genome of 150 individuals. Population genetic polymorphism and structure analysis reveal that Jinchuan yak can be differentiated as a unique and original yak population among the domestic yak. Combined with geological change, the Jinchuan yak's evolutionary origin is speculated to be about 6290 years ago, which may be related to the unique geographical environment of the eastern edge of the QTP during this period. Compared with other domestic yaks, this new population has 280 positively selected genes. The genes related to skeletal function hold a considerable and remarkable proportion, suggesting that the specific skeletal characteristics have been enhanced in the adaptive evolution of Jinchuan yak in the extreme plateau environment. The genome-wide association study has revealed that TUBA8 and TUBA4A, the genes that regulate the cytoskeleton, are potential genes associated with the multirib trait. Our findings provide a basis to further understand the generation mechanism of the adaptive evolution of this new population in high-altitude extreme environments and the multivertebrate trait of domestic animals.

Productive Overcompensation of Alpine Meadows in Response to Yak Grazing in the Eastern Qinghai-Tibet Plateau.

Understanding the interaction between large herbivores and pasture production, especially with respect to the grazing optimization hypothesis, is critical for pasture management and generating theoretical and testable predictions. However, the optimization hypothesis remains contradictory in alpine meadows on the Qinghai-Tibet Plateau (QTP). In this study, we tested the grazing optimization hypothesis using four yak-grazing intensities (no grazing, light grazing, moderate grazing and heavy grazing) in alpine meadow habitats from 2015 to 2017. The results indicated that species diversity did not differ significantly among grazing regimes during the experimental period. However, the aboveground net primary production (ANPP) under moderate grazing consistently significantly exceeded that in control enclosures over 3 years, confirming the grazing optimization hypothesis. Levels of overcompensation varied among grazing intensities and years, and grazing-induced plant compensation may only occur in the short term. The enhanced regrowth of Poaceae and Cyperaceae induced by yak grazing might contribute to the overall level of overcompensation by plant community. Our results strongly support the grazing optimization hypothesis in the context of alpine meadows grazed by yaks, emphasizing the complex interactions between ANPP, herbivores and other ecological factors in alpine meadows on the QTP. These findings provide new insights for the development of an ecological conservation strategy that will help restore this fragile ecosystem and balance the seemingly incompatible requirements of animal husbandry.

Genetic Diversity, Molecular Phylogeny, and Selection Evidence of Jinchuan Yak Revealed by Whole-Genome Resequencing.

Jinchuan yak, a newly discovered yak breed, not only possesses a large proportion of multi-ribs but also exhibits many good characteristics, such as high meat production, milk yield, and reproductive performance. However, there is limited information about its overall genetic structure, relationship with yaks in other areas, and possible origins and evolutionary processes. In this study, 7,693,689 high-quality single-nucleotide polymorphisms were identified by resequencing the genome of Jinchuan yak. Principal component and population genetic structure analyses showed that Jinchuan yak could be distinguished as an independent population among the domestic yak population. Linkage disequilibrium analysis showed that the decay rate of Jinchuan yak was the lowest of the domestic yak breeds, indicating that the degree of domestication and selection intensity of Jinchuan yak were higher than those of other yak breeds. Combined with archaeological data, we speculated that the origin of domestication of Jinchuan yak was ∼6000 yr ago (4000-10,000 yr ago). The quantitative dynamics of population growth history in Jinchuan yak was similar to that of other breeds of domestic and wild yaks, but was closer to that of the wild yak. No significant gene exchange between Jinchuan and other domestic yaks occurred. Compared with other domestic yaks, Jinchuan yak possessed 339 significantly and positively selected genes, several of which relate to physiological rhythm, histones, and the breed's excellent production characteristics. Our results provide a basis for the discovery of the evolution, molecular origin, and unique traits of Jinchuan yak.

Author Correction: Correlations between gut microbiota community structures of Tibetans and geography.

A correction to this article has been published and is linked from the HTML version of this paper. The error has been fixed in the paper.

Correlations between gut microbiota community structures of Tibetans and geography.

Microbial communities of human gut directly influence health and bear adaptive potential to different geography environment and lifestyles. However, knowledge about the influences of altitude and geography on the gut microbiota of Tibetans is currently limited. In this study, fecal microbiota from 208 Tibetans across six different locations were analyzed by MiSeq sequencing; these locations included Gannan, Gangcha, Tianzhu, Hongyuan, Lhasa and Nagqu, with altitudes above sea level ranging from 2800 m to 4500 m across the Tibetan plateau. Significant differences were observed in microbial diversity and richness in different locations. At the phylum level, gut populations of Tibetans comprised Bacteroidetes (60.00%), Firmicutes (29.04%), Proteobacteria (5.40%), and Actinobacteria (3.85%) and were marked by a low ratio (0.48) of Firmicutes to Bacteroidetes. Analysis based on operational taxonomic unit level revealed that core microbiotas included Prevotella, Faecalibacterium, and Blautia, whereas Prevotella predominated all locations, except Gangcha. Four community state types were detected in all samples, and they mainly belong to Prevotella, Bacteroides, and Ruminococcaceae. Principal component analysis and related correspondence analysis results revealed that bacterial profiles in Tibetan guts varied significantly with increasing altitude, BMI, and age, and facultative anaerobes were rich in Tibetan guts. Gut microbiota may play important roles in regulating high-altitude and geographical adaptations.

Toward Understanding the Genetic Basis of Yak Ovary Reproduction: A Characterization and Comparative Analyses of Estrus Ovary Transcriptiome in Yak and Cattle.

BACKGROUND: Yaks (Bos grunniens) are endemic species that can adapt well to thin air, cold temperatures, and high altitude. These species can survive in harsh plateau environments and are major source of animal production for local residents, being an important breed in the Qinghai-Tibet Plateau. However, compared with ordinary cattle that live in the plains, yaks generally have lower fertility. Investigating the basic physiological molecular features of yak ovary and identifying the biological events underlying the differences between the ovaries of yak and plain cattle is necessary to understand the specificity of yak reproduction. Therefore, RNA-seq technology was applied to analyze transcriptome data comparatively between the yak and plain cattle estrous ovaries. RESULTS: After deep sequencing, 3,653,032 clean reads with a total of 4,828,772,880 base pairs were obtained from yak ovary library. Alignment analysis showed that 16992 yak genes mapped to the yak genome, among which, 12,731 and 14,631 genes were assigned to Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, comparison of yak and cattle ovary transcriptome data revealed that 1307 genes were significantly and differentially expressed between the two libraries, wherein 661 genes were upregulated and 646 genes were downregulated in yak ovary. Functional analysis showed that the differentially expressed genes were involved in various Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. GO annotations indicated that the genes related to "cell adhesion," "hormonal" biological processes, and "calcium ion binding," "cation transmembrane transport" molecular events were significantly active. KEGG pathway analysis showed that the "complement and coagulation cascade" pathway was the most enriched in yak ovary transcriptome data, followed by the "cytochrome P450" related and "ECM-receptor interaction" pathways. Moreover, several novel pathways, such as "circadian rhythm," were significantly enriched despite having no evident associations with the reproductive function. CONCLUSION: Our findings provide a molecular resource for further investigation of the general molecular mechanism of yak ovary and offer new insights to understand comprehensively the specificity of yak reproduction.