Weighted Correlation Network Analysis (WGCNA) of Japanese Flounder ( Paralichthys olivaceus ) Embryo Transcriptome Provides Crucial Gene Sets for Understanding Haploid Syndrome and Rescue by Diploidization

Artificial gynogenesis is of great research value in fish genetics and breeding technology. However, existing studies did not explain the mechanism of some interesting phenomena. Severe developmental defects in gynogenetic haploids can lead to death during hatching. After diploidization of chromosomes, gynogenetic diploids may dispense from the remarkable malformation and restore the viability, although the development time is longer and the survival rate is lower compared with normal diploids. The aim of this study was to reveal key mechanism in haploid syndrome of Japanese flounder, a commercially important marine teleost in East Asia. We measured genome-scale gene expression of flounder haploid, gynogenetic diploid and normal diploid embryos using RNA-Seq, constructed a module-centric co-expression network based on weighted correlation network analysis (WGCNA) and analyzed the biological functions of correlated modules. Module gene content analysis revealed that the formation of gynogenetic haploids was closely related to the abnormality of plasma proteins, and the up-regulation of p53 signaling pathway might rescue gynogenetic embryos from haploid syndrome via regulating cell cycle arrest, apoptosis and DNA repair. Moreover, normal diploid has more robust nervous system. This work provides novel insights into molecular mechanisms in haploid syndrome and the rescue process by gynogenetic diploidization.