Characterization of Anaerobic Rumen Fungal Community Composition in Yak, Tibetan Sheep and Small Tail Han Sheep Grazing on the Qinghai-Tibetan Plateau

Abstract

The anaerobic rumen fungal community play a critical role in fibrous material degradation. However, there is a lack of data describing the composition of anaerobic rumen fungal community of full grazing ruminants in the Qinghai-Tibetan Plateau. For this reason, we employed the next-generation sequencing technique to elucidate the rumen fungal structure composition and evaluate the effects of host species on fungal communities. Community comparisons (Bray–Curtis index) between yak and Tibetan sheep revealed that the rumen fungal community was affected by host species (p < 0.05). The alpha diversity indices in the yak were significantly higher than in the Tibetan sheep and Small Tail Han sheep. Neocallimastigomycota was predominant regardless of host species. Within this phylum, unidentified genus of Neocallimastigaceae was the most dominant in all samples, followed by Piromyces and Orpinomyces. Moreover, the shared and unique OTUs in the rumen were identified and most of them belonged to the Orpinomyces. Co-occurrence network analysis identified that each animal species had their own keystone species and most of them were non-dominant flora. Our data indicate that host breeds override living environment as the key factor that determines fungal community in the rumen of grazing ruminants in the Qinghai-Tibetan Plateau.