Whole Genome DNA Methylation Variations in Mammary Gland Tissues from Holstein Cattle Producing Milk with Various Fat and Protein Contents

Milk fat and protein contents are among key elements of milk quality, and they are attracting more attention in response to consumers′ demand for high-quality dairy products. To investigate the potential regulatory roles of DNA methylation underlying milk component yield, whole genome bisulfite sequencing was employed to profile the global DNA methylation patterns of mammary gland tissues from 17 Canada Holstein cows with various milk fat and protein contents. A total of 706, 2420 and 1645 differentially methylated CpG sites (DMCs) were found between high vs. low milk fat (HMF vs. LMF), high vs. low milk protein (HMP vs. LMP), and high vs. low milk fat and protein (HMFP vs. LMFP) groups, respectively (q value < 0.1). Twenty-seven, 56 and 67 genes harboring DMCs in gene regions (denoted DMC genes) were identified for HMF vs. LMF, HMP vs. LMP and HMFP vs. LMFP, respectively. DMC genes from HMP vs. LMP and HMFP vs. LMFP comparisons were significantly overrepresented in GO terms related to aerobic electron transport chain and/or mitochondrial ATP (adenosine triphosphate) synthesis coupled electron transport. A total of 83 (HMF vs. LMF), 708 (HMP vs. LMP) and 408 (HMFP vs. LMFP) DMCs were co-located with 87, 147 and 158 quantitative trait loci (QTL) for milk component and yield traits, respectively. In conclusion, the identified methylation changes are potentially involved in the regulation of milk fat and protein yields, as well as the variation in reported co-located QTLs.