Network analysis reveals different hub genes and molecular pathways for pig in vitro fertilized early embryos and parthenogenotes

Maternal-to-zygotic transition (MZT) occurs when maternal transcripts decay and zygotic genome is activated gradually at the early stage of embryo development. Previously, single-cell RNA-seq (scRNA-seq) has helped us to uncover the MZT-associated mRNA dynamics of in vitro-produced pig early embryos. Here, to further investigate functional modules and hub genes associated with MZT process, the weighted gene co-expression network analysis (WGCNA) was performed on our previously generated 45 scRNA-seq datasets. For the in vitro fertilized embryo (IVF) group, 5 significant modules were identified (midnight blue/black/red and blue/brown modules, positively correlated with 1-cell (IVF1) and 8-cell (IVF8), respectively), containing genes mainly enriched in signalling pathways such as Wnt, regulation of RNA transcription, fatty acid metabolic process, poly(A) RNA binding and lysosome. For the parthenogenetically activated embryo (PA) group, 9 significant modules were identified (black/purple/red, brown/turquoise/yellow, and magenta/blue/green modules, positively correlated with MII oocytes, 1-cell (PA1) and 8-cell (PA8), respectively), mainly enriched in extracellular exosome, poly(A) RNA binding, mitochondrion and transcription factor activity. Moreover, some of identified hub genes within 3 IVF and 9 PA significant modules, including ADCY2, DHX34, KDM4A, GDF10, ABCC10, PAFAH2, HEXIM2, COQ9, DCAF11, SGK1 and ESRRB, have been reported to play vital roles in different biological processes. Our findings provide information and resources for subsequent in-depth study on the regulation and function of MZT in pig embryos.