Age-specific composition of milk microbiota in Tibetan sheep and goats.

This study investigates the dynamic changes in milk nutritional composition and microbial communities in Tibetan sheep and goats during the first 56 days of lactation. Milk samples were systematically collected at five time points (D0, D7, D14, D28, D56) post-delivery. In Tibetan sheep, milk fat, protein, and casein contents were highest on D0, gradually decreased, and stabilized after D14, while lactose and galactose levels showed the opposite trend. Goat milk exhibited similar initial peaks, with significant changes particularly between D0, D7, D14, and D56. 16S rRNA gene sequencing revealed increasing microbial diversity in both species over the lactation period. Principal coordinates analysis identified distinct microbial clusters corresponding to early (D0-D7), transitional (D14-D28), and mature (D56) stages. Core phyla, including Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria, dominated the milk microbiota, with significant temporal shifts. Core microbes like Lactobacillus, Leuconostoc, and Streptococcus were common in both species, with species-specific taxa observed (e.g., Pediococcus in sheep, Shewanella in goats). Furthermore, we observed a highly shared core microbiota in sheep and goat milk, including Lactobacillus, Leuconostoc, and Streptococcus. Spearman correlation analysis highlighted significant relationships between specific microbial genera and milk nutrients. For instance, Lactobacillus positively correlated with total solids, non-fat milk solids, protein, and casein, while Mannheimia negatively correlated with protein content. This study underscores the complex interplay between milk composition and microbial dynamics in Tibetan sheep and goats, informing strategies for livestock management and nutritional enhancement. KEY POINTS: • The milk can be classified into three types based on the microbiota composition • The changes of milk microbiota are closely related to the variations in nutrition • Filter out microbiota with species specificity and age specificity in the milk.

Different feeding strategies can affect growth performance and rumen functions in Gangba sheep as revealed by integrated transcriptome and microbiome analyses.

Due to the harsh environment in the Tibetan Plateau, traditional grazing greatly limits the growth potential of local animals and causes severe ecosystem degradation. This is an urgent issue to be solved, which requires alternative strategies for grazing animals in the Tibetan alpine pastoral livestock systems. This study aimed to investigate the effects of different feeding strategies on growth performance and ruminal microbiota-host interactions in the local breed of sheep (Gangba sheep). Thirty 9-month old Gangba sheep (n = 10 per group) were assigned to natural grazing (G), semi-grazing with supplementation (T), and barn feeding (F) groups (supplementation of concentrate and oat hay) based on body weight. At the end of the experiment (75 d), all sheep were weighed, rumen fluid was obtained from six sheep per group, and ruminal epithelium was obtained from 3 sheep per group. The results showed that: (1) Compared with the G and T groups, the F group significantly increased dry matter intake, average daily gain, and feed conversion ratio of animals. Additionally, Gangba sheep in the F group had higher concentrations of ruminal short-chain volatile fatty acids (VFAs), especially propionate and butyrate (P <0.05) than sheep in the G and T groups. (2) The principal coordinates analysis indicated a significant difference in bacterial composition among different feed strategies. More specifically, the relative abundance of propionate (unidentified F082 and Succiniclasticum) and butyrate-producing (Eubacterium_coprostanoligenes_group) genera were also observed to be increased in the F group, in which unidentified F082 was identified as a differential biomarker among the three groups according to linear discriminant analysis effect size analysis. (3) The dynamics of the rumen epithelial transcriptome revealed that ECM-receptor interactions, focal adhesion, and PI3K-Akt signaling pathways, which are critical in mediating many aspects of cellular functions such as cell proliferation and motility, were upregulated in the F group. In conclusion, under harsh conditions in the Tibetan alpine meadow, barn feeding increased ruminal VFAs concentrations (especially propionate and butyrate), which stimulated gene expression related to cell proliferation in rumen epithelium, appearing to be superior to natural grazing and semi-grazing in gaining body weight of the local Gangba sheep.

Comparative analysis of long noncoding RNA and mRNA expression provides insights into adaptation to hypoxia in Tibetan sheep.

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS_00306477, TCONS_00306029, TCONS_00139593, TCONS_00293272, and TCONS_00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments.

Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep.

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of Fst and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding Fst significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5'UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying "C-INDEL" of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying "CA-INDEL" of HAG1. We concluded that "C-INDEL" variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

Complete mitochondrial genome analysis of Leptomias sp. (Coleoptera, Curculionidae) from Southeast Tibet of China.

Leptomias sp. belongs to Coleoptera, Curculionidae, Leptomias Faust. It is mainly distributed in the Himalayas and Yarlung Zangbo River valleys in Tibet. It is the first time to report the complete mitochondrial genome of Leptomias sp., which featured a typical circular molecule of 16,801 bp in length, including 13 PCGs, 22 tRNAs, two rRNAs, and one control region. The overall nucleotide composition was 38.7% of A, 10.0% of G, 34.3% of T, and 17.0% of C. The content of A + T was higher than G + C. In this study, we also determined that Leptomias sp. was the sister to the Sympiezomias velatus based the phylogenetic analysis of nucleotide sequence datasets. This study provides a useful resource for further studies on conservation and population genetics of Leptomias species.

Characteristics of size distributions and sources of water-soluble ions in Lhasa during monsoon and non-monsoon seasons.

To understand the physical and chemical characteristics, particle size distribution and sources of size-separated aerosols in Lhasa, which is located on the Tibetan Plateau (TP), six sizes of aerosol samples were collected in Lhasa in 2014. Ca2+, NH4+, NO3-, SO42- and Cl- were the dominant ions. The ratio of cation equivalents (CE) to anion equivalents (AE) for each particle size segment indicated that the atmospheric aerosols in Lhasa were alkaline. SO42- and NO3- could be neutralized by Ca2+, but could not be neutralized by NH4+, according to the [NH4+]/[NO3- + SO42-] and [Ca2+]/[NO3- + SO42-] ratios. Mobile sources were dominant in PM0.95-1.5, PM1.5-3 and PM3-7.2, while stationary sources were dominant in the other three size fractions according to the [NO3-]/[SO42-] ratios. The particle size distribution of all water-soluble ions during monsoon and non-monsoon periods was characterized by a bimodal distribution due to the different sources and formation mechanisms, and it was revealed that different ions had different sources in different seasons and different particle size segments by combining particle size distribution with correlation analysis. Source analysis of aerosols in Lhasa was performed using the Principal component analysis (PCA) for the first time, which revealed that combustion sources, motor vehicle exhaust, photochemical reaction sources and various types of dust were the main sources of Lhasa aerosols. Furthermore, Lhasa's air quality was also affected by long-distance transmission, expressed as pollutants from South Asia and West Asia, which were transmitted to Lhasa according to backward trajectory analysis.