As one of the important ruminants of the Qinghai-Tibet Plateau, Tibetan sheep are able to reproduce and maintain their population in this harsh environment of extreme cold and low oxygen. However, the adaptive mechanism of Tibetan sheep when nutrients are scarce in the cold season of the Plateau environment is unclear. We conducted comparative analysis rumen fermentation parameters, rumen microbes, and expression of host genes related to nutrient absorption and rumen epithelial barrier function in cold and warm season Tibetan sheep. We found that concentrations of the volatile fatty acids (VFAs) acetate, propionate and butyrate of Tibetan sheep in the cold season were significantly higher than in the warm season (P < 0.05). Microbial 16S rRNA gene analysis revealed significant differences in rumen microbiota between the cold and warm seasons, and the abundance of microbial in the cold season was significantly higher than that in the warm season (P < 0.05), and the lack of nutrients in the cold season led to a significant reduction in the expression of SGLT1, Claudin-4, and ZO-1 genes in the rumen epithelium. Correlation analysis revealed significant associations of some rumen microorganisms with the fermentation product acetate and the rumen epithelial genes SGLT1, Claudin-4, and ZO-1.

Tibetan sheep can maintain a normal life and reproduce in harsh environments under extreme cold and lack of nutrition. However, the molecular and metabolic mechanisms underlying the adaptability of Tibetan sheep during the cold season are still unclear. Hence, we conducted a comprehensive analysis of rumen epithelial morphology, epithelial transcriptomics, microbiology and metabolomics in a Tibetan sheep model. The results showed that morphological structure of rumen epithelium of Tibetan sheep in cold season had adaptive changes. Transcriptomics analysis showed that the differential genes were primarily enriched in the PPAR signaling pathway (ko03320), legionellosis (ko05134), phagosome (ko04145), arginine and proline metabolism (ko00330), and metabolism of xenobiotics by cytochrome P450 (ko00980). Unique differential metabolites were identified in cold season, such as cynaroside A, sanguisorbin B and tryptophyl-valine, which were mainly enriched in arachidonic acid metabolism, arachidonic acid metabolism and linolenic acid metabolism pathways, and had certain correlation with microorganisms. Integrated transcriptome-metabolome-microbiome analysis showed that epithelial gene-GSTM3 expression was upregulated in the metabolism of xenobiotics by the cytochrome P450 pathway during the cold season, leading to the downregulation of some harmful metabolites; TLR5 gene expression was upregulated and CD14 gene expression was downregulated in the legionellosis pathway during the cold season. This study comprehensively described the interaction mechanism between the rumen host and microbes and their metabolites in grazing Tibetan sheep during the cold season. Rumen epithelial genes, microbiota and metabolites act together in some key pathways related to cold season adaptation.

Keratin-associated proteins (KAPs) are structural components of wool and hair fibres, and are believed to play a role in defining the physico-mechanical properties of the wool fibre. In this study, the putative ovine homologue of the human KAP26-1 gene (KRTAP26-1) was sequenced and four variants (named A–D) were identified. The sequences shared some identity with each other and with other KRTAPs, but they had the greatest similarity with the human KRTAP26-1 sequence. This suggests they represent different variants of ovine KRTAP26-1. The association of these KRTAP26-1 variants with wool traits was investigated in the 383 Merino-Southdown cross sheep. The presence of B was associated (p < 0.05) with an increase in mean fibre diameter (MFD), mean fibre curvature, and prickle factor (PF). The presence of C was found to be associated (p < 0.05) with an increase in wool yield (Yield) and mean staple length (MSL), and a decrease in MFD, fibre diameter standard deviation (FDSD), and PF. The results suggest that sheep with C have, on average, higher wool quality. These results may be useful in the future development of breeding programs based on decreasing wool MFD and FDSD, or on increasing wool MSL.

The gene encoding the high glycine/tyrosine keratin-associated protein 20-2 (KAP20-2) gene has been described in humans, but has not been identified in any livestock species. A search for similar sequences in the caprine genome using the human KAP20-2 gene (KRTAP20-2) revealed a homologous sequence on chromosome 1. Three different banding patterns representing distinct sequences (A–C) in Longdong cashmere goats were identified using polymerase chain reaction-single stranded conformational polymorphism (PCR-SSCP) analysis. These sequences shared high sequence similarity with the human and mouse KRTAP20-2 sequences, suggesting that A–C are caprine variants of the human and mouse genes. Four single nucleotide polymorphisms (SNPs) were identified, and three of them were non-synonymous. KRTAP20-2 was found to be expressed in secondary hair follicles, but not in heart, liver, lung, kidney, spleen, or longissimus dorsi muscle. The presence of A was associated with increased cashmere fibre weight, while the presence of B was associated with a decrease in cashmere fibre weight and curly fibre length. Goats with genotype AA had a higher cashmere fibre weight and a higher curly fibre length than those with genotypes AB or BB. These results indicate that caprine KRTAP20-2 variation may have value as a genetic marker for improving cashmere fibre weight.

As an important ruminant on the Qinghai-Tibet Plateau, Tibetan sheep can maintain their population reproduction rate in the harsh high-altitude environment of low temperature and low oxygen, which relies on their special plateau adaptations mechanism that they have formed for a long time. Microbiomes (known as “second genomes”) are closely related to the nutrient absorption, adaptability, and health of the host. In this study, rumen fermentation characteristics, the microbiota, and rumen epithelial gene expression of Tibetan sheep in various months were analyzed. The results show that the rumen fermentation characteristics of Tibetan sheep differed in different months. The total SCFAs (short-chain fatty acids), acetate, propionate, and butyrate concentrations were highest in October and lowest in June. The CL (cellulase) activity was highest in February, while the ACX (acid xylanase) activity was highest in April. In addition, the diversity and abundance of rumen microbes differed in different months. Bacteroidetes (53.4%) and Firmicutes (27.4%) were the dominant phyla. Prevotella_1 and Rikenellaceae_RC9_gut_group were the dominant genera. The abundance of Prevotella_1 was highest in June (27.8%) and lowest in December (17.8%). In addition, the expression of CLAUDIN4 (Claudin-4) and ZO1 (Zonula occludens 1) was significantly higher in April than in August and December, while the expression of SGLT1 (Sodium glucose linked transporter 1) was highest in August. Correlation analysis showed that there were interactions among rumen fermentation characteristics, the microbiota, and host gene expression, mainly by adjusting the amino acid metabolism pathway and energy metabolism pathway to improve energy utilization. At the same time, we adjusted the balance of the rumen “core microbiota” to promote the development of rumen and maintain the homeostasis of rumen environment, which makes Tibetan sheep better able to adapt to the harsh environment in different periods of the Qinghai-Tibet Plateau.