Dynamic transcriptome analysis of Bombyx mori embryonic development

Bombyx mori has been extensively studied but the gene expression control of its embryonic development is unclear. In this study, we performed transcriptome profiling of six stages of B. mori embryonic development using RNA sequencing (RNA-seq). A total of 12 894 transcripts were obtained from the embryos. Of these, 12 456 transcripts were shared among the six stages, namely, fertilized egg, blastoderm, germ-band, organogenesis, reversal period, and youth period stages. There were 111, 48, 41, 54, 77, and 107 transcripts specifically expressed during the six stages, respectively. By analyzing weighted gene correlation networks and differently expressed genes, we found that during embryonic development, many genes related to DNA replication, transcription, protein synthesis, and epigenetic modifications were upregulated in the early embryos. Genes of cuticle proteins, chitin synthesis-related proteins, and neuropeptides were more abundant in the late embryos. Although pathways of juvenile hormone and the ecdysteroid 20-hydroxyecdysone, and transcription factors were expressed throughout the embryonic development stages, more regulatory pathways were highly expressed around the organogenesis stage, suggesting more gene expression for organogenesis. The results of RNA-seq were confirmed by quantitative real-time polymerase chain reaction of 16 genes of different pathways. Nucleic acid methylation and seven sites in histone H3 modifications were confirmed by dot blot and western blot. This study increases the understanding of the molecular mechanisms of the embryonic developmental process and information on the regulation of B. mori development.