Single-cell RNA sequencing reveals developmental heterogeneity of blastomeres during major genome activation in bovine embryos

Embryonic development is initially controlled by maternal RNAs and proteins stored in the oocyte, until gene products gradually generated by the embryo itself take over. Major embryonic genome activation (EGA) in bovine embryos occurs at the eight- to 16-cell stage. Morphological observations, such as size of blastomeres and distribution of microvilli, suggested heterogeneity among individual cells already at this developmental stage. To address cell heterogeneity on the transcriptome level, we performed single-cell RNA sequencing of 161 blastomeres from 14 in vitro produced bovine embryos at Day 2 (n = 6) and Day 3 (n = 8) post fertilization. Complementary DNA libraries were prepared using the Single-Cell RNA-Barcoding and Sequencing protocol and sequenced. Non-supervised clustering of single-cell transcriptome profiles identified six clusters with specific sets of genes. Most embryos were comprised of cells from at least two different clusters. Sorting cells according to their transcriptome profiles resulted in a non-branched pseudo-time line, arguing against major lineage inclination events at this developmental stage. In summary, our study revealed heterogeneity of transcriptome profiles among single cells in bovine Day 2 and Day 3 embryos, suggesting asynchronous blastomere development during the phase of major EGA.