Transcriptome profiling in rumen, reticulum, omasum, and abomasum tissues during the developmental transition of pre-ruminant to the ruminant in yaks

Abstract

The development of the four stomachs of yak is closely related to its health and performance, however the underlying molecular mechanisms are largely unknown. Here, we systematically analyzed mRNAs of four stomachs in five growth time points [0 day, 20 days, 60 days, 15 months and 3 years (adult)] of yaks. Overall, the expression patterns of DEmRNAs were unique at 0 d, similar at 20 d and 60 d, and similar at 15 m and adult in four stomachs. The expression pattern in abomasum was markedly different from that in rumen, reticulum and omasum. Short Time-series Expression Miner (STEM) analysis demonstrated that multi-model spectra are drastically enriched over time in four stomachs. All the identified mRNAs in rumen, reticulum, omasum and abomasum were classified into 6, 4, 7, and 5 cluster profiles, respectively. Modules 9, 38, and 41 were the most significant three colored modules. By weighted gene co-expression network analysis (WGCNA), a total of 5,486 genes were categorized into 10 modules. CCKBR, KCNQ1, FER1L6, and A4GNT were the hub genes of the turquoise module, and PAK6, TRIM29, ADGRF4, TGM1, and TMEM79 were the hub genes of the blue module. Furthermore, functional KEGG enrichment analysis suggested that the turquoise module was involved in gastric acid secretion, sphingolipid metabolism, ether lipid metabolism, etc., and the blue module was enriched in pancreatic secretion, pantothenate and CoA biosynthesis, and starch and sucrose metabolism, etc. Our study aims to lay a molecular basis for the study of the physiological functions of rumen, reticulum, omasum and abomasum in yaks. It can further elucidate the important roles of these mRNAs in regulation of growth, development and metabolism in yaks, and to provide a theoretical basis for age-appropriate weaning and supplementary feeding in yaks.