Abstract
The milk fat globule (MFG) is one of the most representative of mammary gland tissues and can be utilized to study gene expression of lactating cows during lactation. In this study, RNA-seq technology was employed to detect differential expression of genes in MFGs at day 10 and day 70 after calving between two groups of cows with extremely high (H group) and low (L group) 305-day milk yield, milk fat yield and milk protein yield. In total, 1232, 81, 429 and 178 significantly differentially expressed genes (false discovery rate q < 0.05) were detected between H10 and L10, H70 and L70, H10 and H70, and L10 and L70 respectively. Gene Ontology enrichment and pathway analysis revealed that these differentially expressed genes were enriched in biological processes involved in mammary gland development, protein and lipid metabolism process, signal transduction, cellular process, differentiation and immune function. Among these differentially expressed genes, 178 (H10 vs. L10), 4 (H70 vs. L70), 68 (H10 vs. H70) and 22 (L10 vs. L70) were found to be located within previously reported QTL regions for milk production traits. Based on these results, some promising candidate genes for milk production traits in dairy cattle were suggested.