Deregulation of oxidative phosphorylation pathways in embryos derived in vitro from prepubertal and pubertal heifers based on whole-transcriptome sequencing

Abstract

Abstract Background Although, oocytes from prepubertal donors are known to be less developmentally competent than those from adult donors it does not restrain their ability to produce full-term pregnancies. The transcriptomic profile of embryos could be used as a predictor for embryo’s individual developmental competence. The aim of the study was to compare transcriptomic profile of blastocysts derived from prepubertal and pubertal heifers oocytes. Bovine cumulus-oocyte complexes (COCs) were obtained by ovum pick- up method from prepubertal and pubertal heifers. After in vitro maturation COCs were fertilized and cultured to the blastocyst stage. Total RNA was isolated from both groups of blastocysts and RNA-seq was performed. Gene ontology analysis was performed by DAVID (Database for Annotation, Visualization and Integrated Discovery). Results A higher average blastocyst rate was obtained in the pubertal than in the pre-pubertal group. There were no differences in the quality of blastocysts between the examined groups. We identified 436 differentially expressed genes (DEGs) between blastocysts derived from researched groups, of which 247 DEGs were downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes, and 189 DEGs were upregulated. The genes involved in mitochondrial function, including oxidative phosphorylation (OXPHOS) were found to be different in studied groups using Kyoto Encyclopedia of Genes (KEGG) pathway analysis and 8 of those DEGs were upregulated and 1 was downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes. DEGs associated with mitochondrial function were found: ATP synthases (ATP5MF-ATP synthase membrane subunit f, ATP5PD- ATP synthase peripheral stalk subunit d, ATP12A- ATPase H+/K + transporting non-gastric alpha2 subunit), NADH dehydrogenases (NDUFS3- NADH: ubiquinone oxidoreductase subunit core subunit S3, NDUFA13- NADH: ubiquinone oxidoreductase subunit A13, NDUFA3- NADH: ubiquinone oxidoreductase subunit A3), cytochrome c oxidase (COX17), cytochrome c somatic (CYCS) and ubiquinol cytochrome c reductase core protein 1 (UQCRC1). We found lower number of apoptotic cells in blastocysts derived from oocytes collected from prepubertal than those obtained from pubertal donors. Conclusions Despite decreased expression of genes associated with OXPHOS pathway in blastocysts from prepubertal heifers oocytes, the increased level of ATP12A together with the lower number of apoptotic cells in these blastocysts might support their survival after transfer.