Endangered wild yak (Bos grunniens) in the Tibetan plateau and adjacent regions: Population size, distribution, conservation perspectives and its relation to the domestic subspecies

Identifying climate refugia for wild yaks (Bos mutus) on the Tibetan Plateau

Climate change is threatening fragile alpine ecosystems and their resident ungulates, particularly the wild yak (Bos mutus) that inhabits alpine areas between the tree line and glaciers on the Tibetan Plateau. Although wild yaks tend to shift habitats in response to changes in climatic factors, the precise impacts of climate change on their habitat distribution and climate refugia remain unclear. Based on over 1000 occurrence records, the maximum entropy (MaxEnt) algorithm was applied to simulate habitat ranges in the last glacial maximum (LGM), Mid-Holocene, current stage, and three greenhouse gas emission scenarios in 2070. Three habitat patches were identified as climate refugia for wild yaks that have persisted from the LGM to the present and are projected to persist until 2070. These stable areas account for approximately 64% of the current wild yak habitat extent and are sufficiently large to support viable populations. The long-term persistence of these climate refugia areas is primarily attributed to the unique alpine environmental features of the Tibetan Plateau, where relatively stable arid or semi-arid climates are maintained, and a wide range of forage resource supplies are available. However, habitat loss by 2070 caused by insufficient protection is predicted to lead to severe fragmentation in the southeastern and northwestern Kunlun, Hengduan, central-western Qilian, and southern Tanggula-northern Himalaya Mountains. Habitat disturbance has also been caused by increasing anthropogenic effects in the southern Tanggula and northern Himalaya Mountains. We suggest that sufficient protection, transboundary cooperation, and community involvement are required to improve wild yak conservation efforts. Our combined modeling method (MaxEnt-Zonation-Linkage Mapper-FRAGSTAT) can be utilized to identify priority areas and linkages between habitat patches while assessing the conservation efficiency of protected areas and analyzing the coupled relationship between climate change and anthropogenic impacts on the habitat distribution of endangered species.

Impact of infectious diseases on wild bovidae populations in Thailand: insights from population modelling and disease dynamics

The wildlife and livestock interface is vital for wildlife conservation and habitat management. Infectious diseases maintained by domestic species may impact threatened species such as Asian bovids, as they share natural resources and habitats. To predict the population impact of infectious diseases with different traits, we used stochastic mathematical models to simulate the population dynamics over 100 years for 100 times in a model gaur (Bos gaurus) population with and without disease. We simulated repeated introductions from a reservoir, such as domestic cattle. We selected six bovine infectious diseases; anthrax, bovine tuberculosis, haemorrhagic septicaemia, lumpy skin disease, foot and mouth disease and brucellosis, all of which have caused outbreaks in wildlife populations. From a starting population of 300, the disease-free population increased by an average of 228% over 100 years. Brucellosis with frequency-dependent transmission showed the highest average population declines (-97%), with population extinction occurring 16% of the time. Foot and mouth disease with frequency-dependent transmission showed the lowest impact, with an average population increase of 200%. Overall, acute infections with very high or low fatality had the lowest impact, whereas chronic infections produced the greatest population decline. These results may help disease management and surveillance strategies support wildlife conservation.

Mapping conservation priorities for wild yak (Bos mutus) habitats on the Tibetan Plateau, China

The wild yak (Bos mutus) is a cold-tolerant herbivore native to the Tibetan Plateau and has been categorized as vulnerable by the International Union for Conservation of Nature and Natural Resources. Low population densities within currently fragmented habitats and unclear landscape conservation priorities warrant attention. Herein, we employed the maximum entropy (MaxEnt) model using over 900 wild yak occurrence records to model wild yak habitat suitability. Our analysis revealed unprotected wild yak landscapes covering 30.79 % of the habitat area, indicating a conservation gap between protected areas (PAs) and wild yak habitats. To protect metapopulation dynamics and mitigate high risks of poaching, habitat degradation and fragmentation, resource competition, and degenerated genetic characterization of wild yaks in fragmented and degraded habitat, we identified eight habitat patches as landscape conservation units (LCUs) and 14 linkages among the LCUs, enhancing the connectivity between LCUs to decrease negative effects of genetic threats. A centrality analysis demonstrated that Changtang, Arjinshan, and Hoh Xil national nature reserves and their linkages are all critical for the maintenance of habitat connectivity. Here, we suggest that habitat- and LCU-specific conservation strategies should be highlighted during the establishment of PAs and transboundary cooperation. Ultimately, our results can assist conservationists and land managers in comprehending wild yak distribution, movement, and habitat requirements, as well as for the development of effective protection strategies. Furthermore, the combined modeling method (MaxEnt-Zonation-InVEST) could be utilized as a component for identifying conservation priorities and linkages between core patches for species and assessing the efficiency of PAs, core habitats, and corridors in achieving conservation goals. Our study can provide a framework in identifying priority conservation and connectivity between habitat patches to facilitate effectively conservation and genetic resilience for endangered species in fragmented habitats.

Establishment of SV40 Large T-Antigen-Immortalized Yak Rumen Fibroblast Cell Line and the Fibroblast Responses to Lipopolysaccharide

The yak lives in harsh alpine environments and the rumen plays a crucial role in the digestive system. Rumen-associated cells have unique adaptations and functions. The yak rumen fibroblast cell line (SV40T-YFB) was immortalized by introducing simian virus 40 large T antigen (SV40T) by lentivirus-mediated transfection. Further, we have reported the effects of lipopolysaccharide (LPS) of different concentrations on cell proliferation, extracellular matrix (ECM), and proinflammatory mediators in SV40T-YFB. The results showed that the immortalized yak rumen fibroblast cell lines were identified as fibroblasts that presented oval nuclei, a fusiform shape, and positive vimentin and SV40T staining after stable passage. Chromosome karyotype analysis showed diploid characteristics of yak (n = 60). LPS at different concentrations inhibited cell viability in a dose-dependent manner. SV40T-YFB treated with LPS increased mRNA expression levels of matrix metalloproteinases (MMP-2 and MMP-9), inflammatory cytokines (TNF-α, IL-1β, IL-6), and urokinase-type plasminogen activator system components (uPA, uPAR). LPS inhibits the expression of tissue inhibitors of metalloproteinases (TIMP-1 and TIMP-2), plasminogen activator inhibitor-2 (PAI-2), fibronectin (FN), anti-inflammatory factor IL-10, and collagen I (COL I) in SV40T-YFB. Overall, these results suggest that LPS inhibits cell proliferation and induces ECM degradation and inflammatory response in SV40T-YFB.

Effects of Maize Varieties on Biomass Yield and Silage Quality of Maize–Soybean Intercropping in the Qinghai–Tibet Plateau

Forage deficiency is the bottleneck that restricts the development of plateau animal husbandry. Maize (Zea mays L.)–soybean (Glycine max L.) intercropping can improve the forage biomass yield and silage quality. This experiment was conducted in Ganzi Tibetan Autonomous Prefecture to explore the effects of four maize varieties (M1, Rongyu Silage No. 1; M2, Yayu 04889; M3, Demeiya No. 1; M4, Zhenghong 505) on biomass yield, nutritional composition, and silage quality in maize–soybean intercropping. The results showed that M1S had the highest total dry matter yield (18.03 t ha−1), M3S had the highest crude protein (CP) content (8.46% DM), and soybeans had the highest water-soluble carbohydrate (WSC) content (8.55% DM). After silage, the CP content (13.44% DM) of mixed silage in M3S was higher, and the contents of neutral detergent fiber (39.42% DM) and acid detergent fiber (25.42% DM) were lower than those in maize silage alone. The WSC content (4.45% DM) of mixed silage in M3S was higher and the pH value (4.46) and ammonia–nitrogen to total nitrogen (3.97%) were lower than those of soybean silage alone. The results of membership function analysis showed that M3S was the best in fresh feeding and silage utilization, followed by M1S. Therefore, M3S (Demeiya No 1. intercropped with soybeans) is recommended in high-altitude areas.

Possible Consequences of Climate Change on Survival, Productivity and Reproductive Performance, and Welfare of Himalayan Yak (Bos grunniens)

Simple Summary

Climate change is a global issue, with a wide range of ecosystems being affected by changing climatic conditions including the Himalaya. Yak are exquisitely adapted to the high-altitude conditions of the Himalaya and are thus highly likely to be affected by climate change. This paper reviews the evidence of how the reported impacts of climate change on the environment and ecosystem of the Himalaya are affecting the survival, productivity and welfare of Himalayan Yak. This review identified that we do not know how big the impact of climate change is on yak as very few papers have measured that impact and, in many cases, potentially climate-change-related effects (such as changes in feed supply) are principally driven by human factors.

Abstract

Yak are adapted to the extreme cold, low oxygen, and high solar radiation of the Himalaya. Traditionally, they are kept at high altitude pastures during summer, moving lower in the winter. This system is highly susceptible to climate change, which has increased ambient temperatures, altered rainfall patterns and increased the occurrence of natural disasters. Changes in temperature and precipitation reduced the yield and productivity of alpine pastures, principally because the native plant species are being replaced by less useful shrubs and weeds. The impact of climate change on yak is likely to be mediated through heat stress, increased contact with other species, especially domestic cattle, and alterations in feed availability. Yak have a very low temperature humidity index (52 vs. 72 for cattle) and a narrow thermoneutral range (5–13 °C), so climate change has potentially exposed yak to heat stress in summer and winter. Heat stress is likely to affect both reproductive performance and milk production, but we lack the data to quantify such effects. Increased contact with other species, especially domestic cattle, is likely to increase disease risk. This is likely to be exacerbated by other climate-change-associated factors, such as increases in vector-borne disease, because of increases in vector ranges, and overcrowding associated with reduced pasture availability. However, lack of baseline yak disease data means it is difficult to quantify these changes in disease risk and the few papers claiming to have identified such increases do not provide robust evidence of increased diseases. The reduction in feed availability in traditional pastures may be thought to be the most obvious impact of climate change on yak; however, it is clear that such a reduction is not solely due to climate change, with socio-economic factors likely being more important. This review has highlighted the large potential negative impact of climate change on yak, and the lack of data quantifying that impact. More research on the impact of climate change in yak is needed. Attention also needs to be paid to developing mitigating strategies, which may include changes in the traditional system such as providing shelter and supplementary feed and, in marginal areas, increased use of yak–cattle hybrids.

Phylogeographical analysis shows the need to protect the wild yaks' last refuge in Nepal

The wild yak Bos mutus was believed to be regionally extinct in Nepal for decades until our team documented two individuals from Upper Humla, north-western Nepal, in 2014. The International Union for Conservation of Nature (IUCN) seeks further evidence for the conclusive confirmation of that sighting. We conducted line transects and opportunistic sign surveys in the potential wild yak habitats of Humla, Dolpa, and Mustang districts between 2015 and 2017 and collected genetic samples (present and historic) of wild and domestic yaks Bos grunniens. We also sighted another wild yak in Upper Humla in 2015. Phylogenetic and haplotype network analyses based on mitochondrial D-loop sequences (~450 bp) revealed that wild yaks in Humla share the haplotype with wild yaks from the north-western region of the Qinghai-Tibetan Plateau in China. While hybridization with domestic yaks is a major long-term threat, illegal hunting for meat and trophy put the very small populations of wild yaks in Nepal at risk. Our study indicates that the unprotected habitat of Upper Humla is the last refuge for wild yaks in Nepal. We recommend wild yak conservation efforts in the country to focus on Upper Humla by (i) assigning a formal status of protected area to the region, (ii) raising awareness in the local communities for wild yak conservation, and (iii) providing support for adaptation of herding practice and pastureland use to ensure the viability of the population.

Physicochemical properties, digestibility and anti-osteoporosis effect of yak bone powder with different particle sizes

As a kind of promising resource, animal bone has been widely processed into functional foods. However, there is little research about the effect of particle size on the physicochemical properties and digestibility of yak bone powder (YBP), as well as its anti-osteoporosis activity. In this study, the YBP with median particle sizes (MPS) ranging from 19.68 to 128.37 μm were prepared, and their digestibility and anti-osteoporosis activity were investigated. The results showed that smaller MPS significantly increased water holding capacity and protein solubility without changing composition. The MPS reduction greatly promoted protein digestion, producing more peptides<3 kDa and free amino acids while decreased Ca²⁺ and P⁵⁺ release during gastrointestinal digestion. The in vivo results revealed the positive effect of YBP on ovariectomy-induced osteoporosis in rats. The bone mineral density of ovariectomized (OVX) rats was obviously improved by regulating bone turnover markers (B-ALP, OCN, S-CTX, ES and TRAP), thus potentially shedding light on osteoporosis remission. However, different MPS exhibited a weak effect on osteoporosis in OVX rats. Therefore, YBP could be produced in relatively large particle size without sacrificing food sensory quality, the processing time of which could also be shortened for higher productivity and lower cost.

The complete mitochondrial genome and phylogenetic analysis of Yanglong yak (Bos grunniens)

In this study, we assembled the complete mitochondrial genome of Yanglong yak (Bos grunniens) from Illumina sequencing reads. The mitochondrial genome is 16,323 bp long with an A + T-biased nucleotide composition, and encodes 13 protein-coding, 22 tRNA, and two rRNA genes along with a noncoding control region. In addition, its gene order is identical to those of the previously published mitochondrial genomes of its congeners. Phylogenetic analysis indicates that this breed is closely related to Datong yak, Pamir yak, Tianjun yak, polled yak, Seron yak, Sunnan yak, a series of Domestic Yak and wild yak, followed by Jinchuan yak and Gannan yak, and slightly far away from Bison and Bos taurus.

Domestication Shapes the Community Structure and Functional Metagenomic Content of the Yak Fecal Microbiota

Domestication is a key factor of genetic variation; however, the mechanism by which domestication alters gut microbiota is poorly understood. Here, to explore the variation in the structure, function, rapidly evolved genes (REGs), and enzyme profiles of cellulase and hemicellulose in fecal microbiota, we studied the fecal microbiota in wild, half-blood, and domestic yaks based on 16S rDNA sequencing, shotgun-metagenomic sequencing, and the measurement of short-chain-fatty-acids (SCFAs) concentration. Results indicated that wild and half-blood yaks harbored an increased abundance of the phylum Firmicutes and reduced abundance of the genus Akkermansia, which are both associated with efficient energy harvesting. The gut microbial diversity decreased in domestic yaks. The results of the shotgun-metagenomic sequencing showed that the wild yak harbored an increased abundance of microbial pathways that play crucial roles in digestion and growth of the host, whereas the domestic yak harbored an increased abundance of methane-metabolism-related pathways. Wild yaks had enriched amounts of REGs in energy and carbohydrate metabolism pathways, and possessed a significantly increased abundance of cellulases and endohemicellulases in the glycoside hydrolase family compared to domestic yaks. The concentrations of acetic, propionic, n-butyric, i-butyric, n-valeric, and i-valeric acid were highest in wild yaks. Our study displayed the domestic effect on the phenotype of composition, function in gut microbiota, and SCFAs associated with gut microbiota, which had a closely association with the growth performance of the livestock. These findings may enlighten the researchers to construct more links between economic characteristics and gut microbiota, and develop new commercial strains in livestock based on the biotechnology of gut microbiota.

Prevalence of bovine viral diarrhea virus (BVDV) in yaks between 1987 and 2019 in mainland China: A systematic review and meta-analysis

BACKGROUND

Bovine viral diarrhea is an infectious disease that causes symptoms such as bovine diarrhea and abortion. It can cause severe losses to the animal husbandry, and the overall epidemic situation of yak's BVDV in China is unclear. Meta-analysis can reveal the basic epidemic situation of BVDV in different yak distribution areas in China, and estimate potentially related factors, to pave the way for clarifying the epidemic situation of yak in the domestic scope.

METHODS

We proceeded to a systematic review and meta-analysis of data from papers on the BVDV incidence and prevalence in yaks in China by searching PubMed, ScienceDirect, Chinese Web of Knowledge (CNKI), Wanfang, and Chongqing VIP for publication from 1987 to 2019. We excluded reviews and duplicate studies, 24 studies denouncing the prevalence of BVDV in yak in China were selected upon our inclusion criterion finally. We estimated the pooled prevalence of BVDV infection in yaks by a random-effects model and evaluated its overall infection burden in China.

FINDINGS

In total, the pooled prevalence of BVDV in yaks in China was 36.0% (95% CI 25.6%-46.4%) based on the data obtained from 13,446 yaks, by detecting antigens and antibodies. The highest BVDV positive rate in yak reached 67.5% in Xinjiang province of China. The prevalence in the six provinces of China was validated to be quite variable (24.4%-67.5%) and reached 369% in yaks of northwest China. Besides, the BVDV antigen-positive rate was estimated at 13.8% (95% CI 8.6%-19.0%) based on 5 studies, comparatively, the pooled BVDV antibody-based on 18 studies was about 32.9% (95% CI 24.6%-41.2%) in China.

INTERPRETATION

This systematic review and meta-analysis firstly established an estimated prevalence of BVDV in yaks in China, as a whole, and estimates potential relevant factors, including geographic location, publication year, age, detection methods, etc., Our findings suggest that the scientific community or decision-makers can formulate corresponding prevention and control plans based on potential risk factors.

Reference gene selection and myosin heavy chain (MyHC) isoform expression in muscle tissues of domestic yak (Bos grunniens)

Domestic yak (Bos grunniens) is the most crucial livestock in the Qinghai-Tibetan Plateau, providing meat and other necessities for local people. The skeletal muscle of adult livestock is composed of muscle fibers, and fiber composition in muscle has influence on meat qualities, such as tenderness, pH, and color. Real-time quantitative polymerase chain reaction (RT-qPCR) is a powerful tool to evaluate the gene expression of muscle fiber, but the normalization of the data depends on the stability of expressed reference genes. Unfortunately, there is no consensus for an ideal reference gene for data normalization in muscle tissues of yak. In this study, we aimed to assess the stability of 14 commonly used candidate reference genes by using five algorithms (GeNorm, NormFinder, BestKeeper, Delat Ct and Refinder). Our results suggested UXT and PRL13A were the most stable reference genes, while the most commonly used reference gene, GAPDH, was most variably expressed across different muscle tissues. We also found that the extensor digitorum lateralis (EDL), trapezius pars thoracica (TPT), and psoas major (PM) muscle had the higher content of type I muscle fibers and the lowest content of type IIB muscle fibers, while gluteobiceps (GB) muscle had the highest content of type IIB muscle fibers. Our study provides the suitable reference genes for accurate analysis of yak muscle fiber composition.

Mitochondrial genome data confirm that yaks can serve as the intermediate host of Echinococcus canadensis (G10) on the Tibetan Plateau

Background

Cervids used to be considered the only animal intermediate hosts of the G10 genotype of Echinococcus canadensis. Yaks are often herded in the Qinghai-Tibet Plateau, China, where echinococcosis remains prevalent. However, no E. canadensis G10 cases have been recorded in yaks until now. The aim of our study was to identify causative agents of echinococcosis in yaks in this region.

Methods

Total genomic DNA was extracted from the germinal layer of one hydatid using a Blood and Tissue Kit. Full-length mitochondrial (mt) cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1) genes were amplified by PCR. All purified PCR products were directly sequenced in both directions. Then seven pairs of overlap primers were designed to amplify the entire mt genome sequence of a suspected E. canadensis G10 isolate. Phylogenetic analyses were performed based on concatenated nucleotides from the 12 protein-coding genes of mt genomes of Echinococcus species in a Bayesian framework using MrBayes v3.1 and implementing the GTR + I + G model.

Results

Hydatids were found in yaks (n = 129) when organs were inspected at the slaughterhouse in Maqu county, Gannan Tibetan Autonomous Prefecture, Gansu Province, China in October 2016. Of these, 33 (25.6%) harbored up to a dozen hydatid cysts. One cyst from each yak was characterized by sequencing its mitochondrial (mt) cox1 and nad1 genes. On the basis of these sequence data, 32 cysts were identified as Echinococcus granulosus (sensu stricto) (G1-G3) and the remaining one was identified as the G10 genotype of E. canadensis. Its mt genome was then fully sequenced and compared with that of the G10 genotype in GenBank (AB745463). Phylogenetic analysis using complete mt genomes confirmed the Chinese cyst as belonging to the G10 genotype.

Conclusions

To our knowledge, this is the first report globally of E. canadensis (G10) from yaks in China, which suggests that the G10 genotype has a wider geographical distribution and broader host range than previously believed. This genotype has therefore potential risks to human health and animal husbandry.