Genome-wide differences in DNA methylation changes in caprine ovaries between oestrous and dioestrous phases

DNA methylation plays a vital role in reproduction. Entire genome DNA methylation changes during the oestrous phase (ES) and dioestrous phase (DS) in the ovaries of Guanzhong dairy goats were investigated using bisulphite sequencing to understand the molecular biological mechanisms of these goats' oestrous cycle.

We discovered distinct genome-wide DNA methylation patterns in ES and DS ovaries. A total of 26,910 differentially methylated regions were upregulated and 21,453 differentially methylated regions were downregulated in the ES samples compared with the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Differentially methylated region analysis showed hypomethylation in the gene body regions and hypermethylation in the joining region between upstream regions and gene bodies. The methylation ratios of the STAR, FGF2, FGF12, BMP5 and SMAD6 genes in the ES samples were lower than those of the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Conversely, the methylation ratios of the EGFR, TGFBR2, IGF2BP1 and MMD2 genes increased in the ES samples compared with the DS samples. In addition, 223 differentially methylated genes were found in the GnRH signalling pathway (KO04912), ovarian steroidogenesis pathway (KO04913), oestrogen signalling pathway (KO04915), oxytocin signalling pathway (KO04921), insulin secretion pathway (KO04911) and MAPK signalling pathway (KO04010).

This study is the first large-scale comparison of the high-resolution DNA methylation landscapes of oestrous and dioestrous ovaries from dairy goats. Previous studies and our investigations have shown that the NR5A2, STAR, FGF2 and BMP5 genes might have potential application value in regulating caprine oestrus.