1 Whole-genome bisulfite sequencing of bovine gametes and invivo-produced pre-implantation embryos

Abstract

Dynamic changes in DNA methylation are crucial in the epigenetic regulation of mammalian embryogenesis. Global DNA methylation studies in the bovine, however, remain mostly at the immunostaining level. We adopted the single-cell whole-genome bisulfite sequencing method to characterise stage-specific genome-wide DNA methylation in bovine sperm, individual oocytes derived invivo and invitro, and invivo-developed embryos at the 2-, 4-, 8-, and 16-cell stages. This method allowed us to theoretically cover all CpG sites in the genome using a limited number of cells from single embryos. Pools of 20 sperm were selected from a bull with proven fertility. Single oocytes (n=6) and embryos (n=4 per stage) were collected from Holstein cows (n=10). Single-cell whole-genome bisulfite sequencing libraries were prepared and sequenced using the Illumina HiSEqn 4000 platform (Illumina, San Diego, CA, USA). Sequencing reads were filtered and aligned to the bovine reference genome (UMD 3.1.1) using Bismark (Krueger and Andrews 2011Bioinformatics27, 1571-1572, DOI: 10.1093/bioinformatics/btr167).A 300-bp tile-based method was applied to bin the genome into consecutive windows to facilitate comparison across samples. The DNA methylation level was calculated as the fraction of read counts of the total number of cytosines (methylated) in the total read counts of reported cytosines and thymines (methylated and unmethylated), only if more than 3 CpG sites were covered in this tile. Gamete-specific differentially methylated regions were identified when DNA methylation levels were greater than 75% in one type of gamete and less than 25% in the other with false discovery rate-corrected Fisher’s exact test P-values of less than 0.05. The major wave of genome-wide DNA demethylation was complete at the 8-cell stage when de novo methylation became prominent. Sperm and oocytes had numerous differentially methylated regions that were enriched in intergenic regions. Differentially methylated regions were also identified between invivo- and invitro-matured oocytes. Moreover, X chromosome methylation followed the global dynamic patterns. Virtually no (less than 1.5%) DNA methylation was found in mitochondrial DNA. Finally, using our RNA sequencing data generated from the same developmental stages (Jiang et al. 2014 BMC Genomics 15, 756; DOI: 10.1186/1471-2164-15-756), we revealed an inverse correlation between gene expression and promoter methylation. Our study provides the first fully comprehensive analysis of the global dynamics of DNA methylation in bovine gametes and single early embryos using single-cell whole-genome bisulfite sequencing. These data provide insights into the critical features of the methylome of bovine embryos and serve as an important reference for embryos produced by assisted reproduction, such as IVF and cloning, and a model for human early embryo epigenetic regulation.