Single-cell RNA sequencing reveals blastomere heterogeneity of 2-cell embryos in pigs

In mammals, single blastomeres from as early as 2-cell embryos demonstrate heterogeneous developmental capacity and fate decision into different cell lineages. However, mechanisms underlying blastomere heterogeneity of 2-cell embryos remain largely unresolved. Here, we analysed the molecular heterogeneity of full-length mRNAs and their 3'UTR regions, based on the single-cell RNA-seq data of pig 2-cell embryos generated from in vivo fertilization (in vivo), in vitro fertilization (in vitro) and parthenogenetic activation (PA), respectively. First, unsupervised clustering helped discover two different groups of blastomeres for 2-cell pig embryos. Between these two groups of blastomeres in pig 2-cell embryos, 35, 301 and 428 full-length mRNAs respectively in in vivo, in vitro and PA embryo types were identified to be differentially expressed (p_adj  ≤ .05 and |log₂ [fold change]| ≥1) (DE mRNAs), while 92, 89 and 42 mRNAs were shown to be with significantly different 3'UTR lengths (3'UTR DE) (p_adj  ≤ .05). Gene enrichment for both DE mRNAs and 3'UTR DE mRNAs found multiple signalling pathways, including cell cycle, RNA processing. Few numbers of common DE mRNAs and 3'UTR DE mRNAs existed between in vitro and in vivo blastomeres derived from 2-cell embryos, indicating the larger differences between in vitro and in vivo fertilized embryos. Integrative genomics viewer analysis further identified that 3'UTRs of HSDL2 and SGTA (in vivo), FAM204A and phosphoserine phosphatase (in vitro), PRPF40A and RPIA (PA) had >100 nt average length changes. Moreover, numbers and locations of regulatory elements (polyadenylation site, cytoplasmic polyadenylation element and microRNA binding sites) within 3'UTRs of these DE mRNAs were predicted. These results indicate that molecular heterogeneity existed among blastomeres from different types of pig 2-cell embryos, providing useful information and novel insights into future functional investigation on its relationship with the subsequent embryo development and differentiation.