Microbial Metagenomes and Host Transcriptomes Reveal the Dynamic Changes of Rumen Gene Expression, Microbial Colonization and Co-Regulation of Mineral Element Metabolism in Yaks from Birth to Adulthood

Yaks are the main pillar of plateau animal husbandry and the material basis of local herdsmen’s survival. The level of mineral elements in the body is closely related to the production performance of yaks. In this study, we performed a comprehensive analysis of rumen epithelial morphology, transcriptomics and metagenomics to explore the dynamics of rumen functions, microbial colonization and functional interactions in yaks from birth to adulthood. Bacteria, eukaryotes, archaea and viruses colonized the rumen of yaks from birth to adulthood, with bacteria being the majority. Bacteroidetes and Firmicutes were the dominant phyla in five developmental stages, and the abundance of genus Lactobacillus and Fusobacterium significantly decreased with age. Glycoside hydrolase (GH) genes were the most highly represented in five different developmental stages, followed by glycosyltransferases (GTs) and carbohydrate-binding modules (CBMs), where the proportion of genes coding for CBMs increased with age. Integrating host transcriptome and microbial metagenome revealed 30 gene modules related to age, muscle layer thickness, nipple length and width of yaks. Among these, the MEmagenta and MEturquoise were positively correlated with these phenotypic traits. Twenty-two host genes involved in transcriptional regulation related to metal ion binding (including potassium, sodium, calcium, zinc, iron) were positively correlated with a rumen bacterial cluster 1 composed of Alloprevotella, Paludibacter, Arcobacter, Lactobacillus, Bilophila, etc. Therefore, these studies help us to understand the interaction between rumen host and microorganisms in yaks at different ages, and further provide a reliable theoretical basis for the development of feed and mineral element supplementation for yaks at different ages.