Identification of lncRNAs involved in maternal-to-zygotic transition of in vitro-produced porcine embryos by single-cell RNA-seq.

Long non-coding RNAs (lncRNAs) function through multiple tiers of molecular circuits and are vital to gamete maturation and early embryo development. However, in pig early embryos, identification and expression dynamics of lncRNAs remain less studied. Here, we systematically analysed the expression dynamics of lncRNAs based on our previously published single-cell RNA-seq data from pig mature oocytes (GSE160334), and single blastomeres biopsied from pig in vitro fertilized (IVF) and early parthenogenetically activated (PA) embryos (1- to 8-cell stages; GSE164812). With the progression of embryo development, the total number of expressed lncRNAs gradually decreased and showed great variation at each developmental stage for both IVF and PA groups. Consecutive stage pairwise comparison of MII oocytes, 1-cell zygotes, 2-cell, 4-cell and 8-cell IVF embryos identified 151, 245, 1119 and 188 differentially expressed (DE) lncRNAs, including 119, 80, 867, 77 up-regulated and 32, 165, 252, 111 down-regulated, while 289, 437, 895 and 495 DE lncRNAs (141, 89, 768, 97 up-regulated and 148, 348, 127, 398 down-regulated) were identified in PA embryos at the same stages. The DE lncRNAs identified within IVF embryos were much different from that identified within PA embryos, showing embryo type-specific manner. Further cross-comparison between PA and IVF embryos identified 184, 656, 2502 and 266 DE lncRNAs for the 1- to 8-cell embryo stages, respectively. Further GO and KEGG enrichment analysis of DE mRNAs targeted by DELs indicated that different signalling pathways were involved in maternal-only and bi-parental embryo development. Collectively, comparative profiling of lncRNA expression dynamics between pig IVF and PA embryos provides a valuable resource, to investigate further regulatory mechanisms of lncRNAs associated with ZGA and maternal RNA decay during early embryo development.