Comparative Transcriptome Analysis of Wing Discs from Bombyx Mori and Bombyx Mandarina

The insect wing is developed from the wing imaginal disc which is designed from the embryonic ectoderm. To get insight into gene expression profiles in wing discs of Bombyx mori during metamorphosis, we compared the gene expression in the wing between B. mori and B. mandarina moth through RNA-seq. Out of total valid reads identified from the 5th day of 5th instar larvae of silkworm (L5), 7th day of pupae (P7), 1st day of moth (M1) and 1st day of wild silkworm moth (WM1), 20,092,004, 29,251,647, 24,654,695 and 19,753,089 reads were mapped to the mRNA reference sequences of silkworm, respectively. 9229, 7048, 9268 and 6701 differentially expressed genes (DEGs) were respectively recorded in P7 vs L5, Ml vs P7, Ml vs L5 and WM1 vs M1. Further, the peroxisome, ribosome, endocytosis and oxidative phosphorylation pathways were significantly regulated in the metamorphosis of the silkworm. Our study identified 16 orthologous genes with a positive selection from M1, which might be subjected to artificial selection in the domestication of B. mori and would play vital roles in the flight of B. mandarina.