78 Effects of invitro production on the epigenome and transcriptome of bovine embryos determined through a multi-omics data integration approach

In cattle, ovarian superovulation followed by invivo embryo collection and transfer (MOET), and the invitro production (IVP) of embryos are used all over the world to improve animal genetics. Application of MOET has resulted in the production of billions of healthy animals during the past 40 years, and IVP has evolved and given rise to significant numbers of calves during the past 10 years. Nevertheless, the use of MOET and IVP can affect the embryo epigenome, and therefore its transcriptome, before and after elongation, as shown by different studies. The integration of publicly available epigenome-transcriptome datasets generated by these studies could lead to a robust characterisation of the impacts of the application of MOET and IVP. The goal of this study was to integrate all publicly available data about MOET and IVP embryos to determine temporally differentially methylated regions (DMRs) and differentially expressed genes (DEGs) from blastocyst to elongation between IVP and MOET embryos. Datasets were downloaded from the Gene Expression Omnibus (GEO) database. Accession numbers were (1) for epigenomics: GSE69173, GSE97517, and GSE101895, plus one provided dataset from O’Doherty et al. (2018 BMC Genomics, 19, 438; https://doi.org/10.1186/s12864-018-4818-3), all hybridized to the EDMA platform GPL18384; (2) for transcriptomics: GSE12327, GSE21030, GSE24596, GSE24936, GSE27817, and GSE40101, all hybridized to the Affymetrix platform GPL2112. Both types of data were analysed with the limma package for R software, and functional enrichment analysis was done with the DAVID database. For DMRs, comparisons between IVP and MOET were made from spherical blastocysts (n=16 per group) on Day 7, to embryos on Day 15, specifically in the trophectoderm (TE) or embryonic disc (ED) regions (n=4 per region and per group). For DEGs, comparisons between IVP and MOET were made from spherical blastocysts (n=9 per group) to elongated blastocysts on Day 13 and embryos undergoing gastrulation on Day 16 (n=6 per group). Considering a P-value <0.05 and fold-change >2, there were 16 672 (TE) and 26 264 (ED) DMRs and 2236 DEGs that temporally differed between IVP and MOET. Most of the identified DMRs were found in intronic regions (around 36%) rather than exonic regions (8%). However, DMRs that were more methylated at IVP compared with MOET contained exons encoding for genes that enriched the Wnt signalling Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway in the ED, and focal adhesion and ECM-receptor interaction KEGG pathways (P<0.05) in the TE. Accordingly, DEGs with lower expression in elongated embryos (Day 13 and Day 16) at IVP as opposed to MOET were mainly associated with these three pathways. In conclusion, this multi-omics analysis demonstrates that even when embryos are produced under different conditions and experiments, the main changes imposed by IVP affected genes involved in embryonic development and adhesion to the endometrium, which could explain the lower survival rates at IVP compared with MOET.