Weighted gene co-expression network analysis of embryos and first instar larvae of the horseshoe crab Tachypleus tridentatus uncovers development gene networks.

Horseshoe crabs are marine chelicerates that have existed on Earth for about 450 million years, and they are often used as an experimental model for studying marine invertebrate embryology. In this study, we performed transcriptome gene expression profiling of four continuous embryonic stages (Stages 18-21) and first instar larvae of Tachypleus tridentatus. A mean of 50,742,995 high-quality clean reads was obtained from each library. We then conducted weighted gene co-expression network analysis (WGCNA) for 13,698 genes with fragments per kilobase of exon per million mapped fragments values >5. We identified 17 modules, six of which likely play critical roles in development. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differentially expressed genes was performed on the biologically significant modules. We found that several pathways, such as hedgehog signaling pathway, VEGF signaling pathway, dorso-ventral axis formation, may be involved in the embryonic development process of T. tridentatus. We also identified hub genes that were highly connected in the six critical modules. This is the first study to apply WGCNA to horseshoe crabs to identify hub genes that may play critical roles in development, and our results provide new insight into the mechanisms underlying early development in horseshoe crabs.