- Research article
- Open Access
Transcriptomic analysis of developmental features of Bombyx mori wing disc during metamorphosis
BMC Genomics volume 15, Article number: 820 (2014)
Wing discs of B. mori are transformed to pupal wings during the larva-to-pupa metamorphosis with dramatic morphological and structural changes. To understand these changes at a transcriptional level, RNA-seq of the wing discs from 6-day-old fifth instar larvae (L5D6), prepupae (PP) and pupae (P0) was performed.
In total, 12,254 transcripts were obtained from the wing disc, out of which 5,287 were identified to be differentially expressed from L5D6 to PP and from PP to P0. The results of comprehensive analysis of RNA-seq data showed that during larvae-to-pupae metamorphosis, many genes of 20E signaling pathway were up-regulated and those of JH signaling pathway were down-regulated. Seventeen transcription factors were significantly up-regulated. Cuticle protein genes (especially wing cuticle protein genes), were most abundant and significantly up-regulated at P0 stage. Genes responsible for the degradation and de novo synthesis of chitin were significantly up-regulated. There were A and B two types of chitin synthases in B. mori, whereas only chitin synthase A was up-regulated. Both trehalose and D-fructose, which are precursors of chitin synthesis, were detected in the hemolymph of L5D6, PP and P0, suggesting de novo synthesis of chitin. However, most of the genes that are related to early wing disc differentiation were down-regulated.
Extensive transcriptome and DGE profiling data of wing disc during metamorphosis of silkworm have been generated, which provided comprehensive gene expression information at the transcriptional level. These results implied that during the larva-to-pupa metamorphosis, pupal wing development and transition might be mainly controlled by 20E signaling in B. mori. The 17 up-regulated transcription factors might be involved in wing development. Chitin required for pupal wing development might be generated from both degradation of componential chitin and de novo synthesis. Chitin synthase A might be responsible for the chitin synthesis in the pupal wing, while both trehalose and D-fructose might contribute to the de novo synthesis of chitin during the formation of pupal wing.
The domestic silkworm, B. mori, is not only an important economic insect, but also a model insect of Lepidoptera. The genome of silkworm has been reported [1, 2], and the genomic resources are available in Silkworm Genome Database: http://silkworm.genomics.org.cn/ and KAIKObase: http://sgp.dna.affrc.go.jp/KAIKObase/[3, 4].
Wing disc is an imaginal disc of insect wing during larval stage. It consists of undifferentiated and proliferating cells . In B. mori, as larvae develop to pupae, the wing discs undergo dramatic morphological changes and evaginate from the body to form the pupal wing. These changes and processes are regulated by insect hormones and hundreds of genes .
20-hydroxyecdysone (20E) and juvenile hormone (JH) are two major hormones in insects and regulate different biological processes, including growth, molting, and reproduction [7, 8]. Both 20E and JH are present at larval stages and they are antagonistic in many aspects. For example, 20E induces periodic ecdysis while JH maintains the "status quo" status of insects at larval stage [9, 10]. Insect metamorphosis is initiated by a surge of 20E at the end of last larval stage when JH titer is low . In the 20E signaling pathway, ecdysone receptor (EcR) binds with ultraspiracle (USP) to form a heterodimer, which binds to the response element of early response genes, such as Broad complex (BR-C), E74 and E75, to initiate the expression. Early-late response genes, including HR3, HR39 and E78, are also involved in the 20E signal transduction pathway . In the JH signal transduction pathway, methoprene-tolerant (Met) was regarded as a potential candidate receptor for JH [14–17]. Knock-out Met led to a down-regulation of JH-response gene Krüppel homolog 1 (Kr-h1), and the insect showed a precocious phenotype in Tribolium castaneum. In B. mori, Met has two isoforms (Met1 and Met2). Met2 combines with steroid receptor coactivator (SRC, also known as "FISC" or "Taiman") to form a heterodimer and induces the expression of Kr-h1. In B. mori, 20E and JH signaling pathways have been found to be involved in normal wing development [19–21], but no detailed analyses of the two signaling pathways was performed at the transcriptional level.
Genes related to wing disc development at the embryonic stage and early stages were well elucidated in Drosophila. Transcription factors, such as Engrailed (En), Apterous (Ap), Scalloped (Sd) and Vestigial (Vg), secreted proteins, such as Hedgehog (Hh), Decapentaplegic (Dpp), Wingless (Wg) and Fringe (Fng), receptors such as Epidermal growth factor receptor (Egfr) and Notch (N), and their ligands such as Serrate (Ser) and Delta (Dl) were found to be involved in the early differentiation and embryonic development of Drosophila wing discs [5, 22–27]. However their roles in wing transition development from wing disc to pupal wing were not reported.
Cuticle protein and chitin are main structural components of cuticle in insects . Thirteen wing cuticle protein genes (WCP) have been identified in B. mori wing disc, including WCP1a, WCP1b, WCP2, WCP3, WCP4, WCP5, WCP6, WCP71a, WCP7b, WCP8, WCP9, WCP10 and WCP11[29–31]. Chitin degradation and synthesis in insects were well-studied . Chitinase, β-N-acetylglucosaminidase and chitin deacetylase are three main enzymes that involved in chitin degradation, while four critical enzymes, trehalase, hexokinase, glucose-6-phosphate isomerase and glutamine and fructose-6-phosphate aminotransferase, are needed for the de novo synthesis of chitin (KEGG pathway map: ko00520). There are two kinds of chitin synthases in insects: chitin synthase 1 (class A) is responsible for the chitin synthesis in the epidermis and trachea, while chitin synthase 2 (class B) is critical for the chitin synthesis in the gut peritrophic membrane .
For the past decade, a variety of methods such as expressed sequence tags (ESTs) [32, 33], serial analysis of gene expression (SAGE) [34, 35], microarrays [36, 37], and RNA-seq technology  have been applied to gene expression analysis in the silkworm, but no information was available on detailed analysis of wing disc development during metamorphosis at the transcriptional level. In this study, RNA-seq method [39, 40] was used to get a comprehensive view on the gene expression of B. mori wing disc during metamorphosis.
Transcriptome sequencing and statistics of gene expression
In total, 4,842,606, 4,741,356 and 4,852,047 clean reads were obtained from the wing disc at 6-day-old fifth instar larvae (L5D6), prepupae (PP) and pupae (P0) stages, respectively. These clean reads were mapped to reference sequences from Silkworm Genome Database using SOAPaligner/soap2 software  and their nucleotide sequences were assembled. A total of 12,254 assembled transcripts were obtained and the average length of the transcripts was 1,226 base pairs.At L5D6, PP and P0 stages 11,316, 11,236, 10,938 assembled transcripts were obtained, respectively (Figure 1). Of the 12,254 assembled transcripts, 9,822 transcripts are commonly shared by the three stages; 607, 373, 312 transcripts are shared by two stages (L5D6 and PP, PP and P0, L5D6 and P0), respectively; 575, 434, 431 transcripts are specifically expressed at L5D6, PP and P0 stage, respectively (Figure 1). These results showed that most transcripts were commonly expressed in each stage, while the stage-specific transcripts might be important for the wing disc development at the specific stage.
Differentially expressed transcripts and gene ontology (GO) classification
In total, 5,287 differentially expressed transcripts were identified based on the standard process (see Methods). From L5D6 to PP (L5D6 vs PP), 2,778 transcripts were identified as differentially expressed, among which 926 were up-regulated and 1,852 were down-regulated (Figure 2). From PP to P0 (PP vs P0), 3,864 transcripts were identified as differentially expressed, of which 672 were up-regulated and 3,192 were down-regulated (Figure 2). It is surprising that there were more down-regulated transcripts than up-regulated transcripts during the transition from L5D6 to PP and from PP to P0 (Additional files 1 and 2), implying that with the transitional development of the wing disc during metamorphosis, expression of more and more genes was turned off.
It is reasonable to hypothesize that developmental features of the wing disc depend on the activities of differentially expressed transcripts. Therefore, GO classification analysis of differentially expressed transcripts may explain the process of wing disc development. The GO ID numbers of differentially expressed transcripts were obtained from Blast2go software and the GO classification histogram was obtained by uploading GO ID number to BGI WEGO online plot software (see Methods). The results showed that binding and catalytic functions were the most abundant GO function items (>10%), both from L5D6 to PP and PP to P0 developmental stages (Figure 3A and 3B). Specifically, in the cellular component ontology, no dominant GO items of the up-regulated transcripts were found for L5D6 vs PP, while several dominant GO items of the down-regulated transcripts were found, including cell, cell part, macromolecular complex, membrane-enclosed lumen, organelle and organelle part (Figure 3A). The dominant GO items of the significantly up-regulated transcripts over down-regulated transcripts for L5D6 vs PP are catalytic, molecular transducer, structural molecular and transporter in molecular function ontology (Figure 3A). The significantly dominant GO item of the down-regulated transcripts over the up-regulated transcripts for L5D6 vs PP in molecular function ontology was transcription regulator, while the down-regulated transcripts of other functional GO items were not significantly more than the up-regulated transcripts in the same catalogs, including translation regulator, in which no up-regulated transcripts were detected for L5D6 vs PP (Figure 3A). For biological process ontology for L5D6 vs PP, the dominant GO items of the up-regulated transcripts included establishment of localization, localization and post-embryonic development, while the dominant GO items of the down-regulated transcripts included cellular component biogenesis, cellular component organization, cellular process and pigmentation (Figure 3A). Although some GO items in the three ontologies, such as synapse, synapse part, antioxidant, electron carrier and translation regulator, showed obviously distinctive differences in the numbers of transcript expression, they were not statistically significant (p < 0.05) by Pearson test  (Figure 3A).On the other hand, for PP vs P0, the dominant GO items of significantly down-regulated transcripts in cellular component ontology included most of catalogues except extracellular (part) and synapse (part), while no dominant GO items of significantly up-regulated transcripts were found (Figure 3B). In molecular function ontology, the only dominant GO item of the significantly up-regulated transcripts was structural molecule; the dominant GO items of the significantly down-regulated transcripts included binding, catalytic, transcription regulator and translation regulator, while the other GO items were not significantly different between up- and down-regulated transcripts (Figure 3B). For cellular process, no dominant GO items of the significantly up-regulated transcripts were detected, while most catalogues were significantly dominant for the down-regulated transcripts over the up-regulated transcripts (Figure 3B).
The 20E signal transduction pathway was more active than the JH signal transduction pathway during metamorphosis
To investigate activity of the 20E and JH signal transduction pathways in the wing disc during metamorphosis, the expression patterns of the genes that are involved in the 20E and JH signal transduction pathways were identified and analyzed.
The results revealed that most of the genes that are associated with the 20E signaling pathway, such as E74, E75 and HR3, were significantly up-regulated during metamorphosis, whereas most of the genes that are involved in the JH signaling pathway, such as SRC and Kr-h1, were significantly down-regulated (Tables 1 and 2). The expression levels of the 20E-related genes appeared to be higher than those of the JH-related genes. This implied that 20E signal is actively involved in the wing disc development during the larva-to-pupa metamorphosis. The assembled nucleotide sequences of transcripts (RNA-seq ID) in Tables 1 and 2 were provided in Additional files 3 and 4, respectively.
Expression of transcription factors in the wing disc during metamorphosis
Transcription factors play critical roles in regulation of expression of their target genes. To identify transcription factors that may regulate the development of B. mori wing discs during metamorphosis, transcription factors were identified by KEGG analysis. In total, 147 out of 12,254 assembled transcripts were identified as transcription factors, 60 of which were differentially expressed during metamorphosis. Among the 60 differentially expressed transcription factors, 17 were up-regulated at PP and/or P0 stage (Table 3). Some of them were selected for confirmation by qRT-PCR analysis. The results showed that the transcription factors Bm_nscaf2589_261 (hormone receptor HR38), Bm_nscaf1898_501 (D-ETS-4-like isoform X2), Bm_nscaf2876_46 (bHLH protein 15-like), Bm_nscaf2888_249 (nuclear receptor GRF), Bm_nscaf2902_040 (protein delilah-like), Bm_nscaf3090_4 (chorion specific C/EBP) and Bm_scaffold316_2 (protein embryonic gonad-like) were significantly up-regulated at PP stage, which were well consistent with the RNA-Seq results (Table 3 and Figure 4A-C, E-H). The transcription factors Bm_nscaf1898_212 (protein pokkuri-like), Bm_nscaf2770_58 (SOX-4-like), Bm_nscaf2847_110 (Krueppel-like factor 10-like) and Bm_nscaf2930_060 (enhancer of split mbeta) were significantly up-regulated at P0 stage, which were also consistent with its RNA-Seq results (Table 3 and Figure 4I-L). The transcription factor Bm_nscaf2847_349 (protein escargot-like) was highly expressed at W1 stage, which was not consistent with the RNA-Seq result (Table 3 and Figure 4D). The qRT-PCR results showed that the RNA-Seq results were generally consistent. Bm_nscaf2859_001 (E75 isoform B), Bm_nscaf2964_066 (hormone receptor 3C) and Bm_nscaf2888_433 (transcription factor E74) were significantly up-regulated at PP or P0 stage (Table 3), which was consisted with the reported results [43, 44]. ASH2 was essential for scales formation in B. mori wing during pupal stage  and we found this transcription factor (Bm_nscaf1898_167) was significantly up-regulated at P0 stage. The assembled nucleotide sequences of transcripts (RNA-seq ID) in Table 3 were provided in Additional file 5, and primers for qRT-PCR were provided in Additional file 6.
Changes in expression level of cuticle protein genes in wing disc during metamorphosis
Cuticle proteins are important component of wing disc in the post-larva development. To study the relationship between expression of cuticle proteins and wing development, cuticle protein genes were analyzed. Totally 205 cuticle protein transcripts were identified using Blast search. Ninety-seven transcripts were expressed at all the three stages, while 7, 9 and 58 transcripts were specifically expressed at L5D6, PP and P0 stage respectively (Figure 5A). One hundred forty eight out of 205 transcripts were differentially expressed, with most of them being up-regulated during the larva-to-pupa metamorphosis, particularly from PP to P0 stage (Figure 5B). This result was consistent with the results by microarray analysis . By similarity analysis of sequences, 10 WCPs genes were identified in our RNA-Seq dataset. Expression of WCP1a, WCP1b, WCP3, WCP4, WCP5, WCP8 and WCP9 was slightly up-regulated at the PP stage and then greatly increased at the P0 stage. Expression of WCP10 was slightly up-regulated at the PP stage and then greatly decreased at the P0 stage, whereas the expression of WCP11 was decreased at the PP stage and then greatly increased at the P0 stage (Table 4).
Change in expression level of chitin degradation and synthesis related genes
To study the dynamic of chitin in wing disc during the pupal wing transition, the genes that are involved in chitin degradation and synthesis pathways were investigated. Most of the genes that are related to chitin degradation were significantly up-regulated at the PP stage and the genes that are related to chitin synthesis were also up-regulated at the PP or P0 stage (Tables 5 and 6). Two critical chitin degradation related genes [Bm_nscaf2993_299 (chitin deacetylase 4) and Bm_nscaf3031_201 (β-N-acetylglucosaminidase 1 precursor)] were significantly up-regulated at the PP stage, then down-regulated at the P0 stage (Figure 6A, B), which were consistent with the results of RNA-Seq (Table 5). Chitin synthesis related genes including Bm_nscaf3027_063 (hexokinase), Bm_nscaf2887_059 (glutamine:F-6-P aminotransferase), Bm_nscaf2823_132 (glucosamine-6-P-N-acetyltransferase), Bm_nscaf2589_266 (chitin synthase A) were significantly up-regulated at the PP or P0 stage, and Bm_nscaf2589_269 (chitin synthase B) was significantly down-regulated from L5D3 to P3 stage (Figure 7A-E), which were generally consistent with the results of RNA-Seq (Table 6). The assembled nucleotide sequences of transcripts (RNA-seq ID) in Tables 5 and 6 were provided in Additional files 7 and 8, respectively. Primers for qRT-PCR were provided in Additional file 6.
De novo synthesis of chitin in insects usually utilizes trehalose as the starting material. However, base on the pathway map of chitin synthesis (KEGG pathway map: ko00520), D-fructose can also be used for chitin de novo synthesis. In this study, high performance liquid chromatography (HPLC) was used to determine the levels of trehalose, D-glucose and D-fructose in the hemolymph of the silkworm at L5D6, PP and P0 stage. The results showed that trehalose, D-glucose and D-fructose were presented with high concentrations at L5D6 stage, and then their concentrations rapid decreased at the PP stage (Table 7). The concentrations of the three sugars were slightly decreased from the PP to P0 stage (Table 7). These results imply that the rapid use of the three sugars in the hemolymph might be taking place from L5D6 to PP stage.
Changes in the expression levels of wing disc development related genes
Expression patterns of the genes homologous to Drosophila wing early development were also analyzed in this study. The results showed that most of the genes, such as En, Ap, Hh, reported to be important for embryonic and early development of Drosophila wing disc were significantly decreased during metamorphosis (Table 8). However, Fng that was necessary for B. mori wing development during larvae to pupae metamorphosis , Ash2 that was vital for scales formation of B. mori wing during pupal stage , and urbain, an unknown gene that was up-regulated at PP stage in B. mori wing , were also significantly up-regulated at the PP or P0 stage in this study. The assembled nucleotide sequences of transcripts (RNA-seq ID) in Table 8 were provided in Additional file 9.
Insect metamorphosis is mediated by 20E and JH signals. During larval stage, JH exhibits its "status quo" function probably through two transcription factors, E75A and Kr-h1, to prevent the activation of BR-C, a 20E signal transduction gene for larvae-pupae transition [7, 47]. Though a relative high titer of JH was present in hemolymph of B. mori at the end of the final instar larvae , the JH does not inhibit the wing disc transformation. Some studies implied that the success of this transition was attributed to the loss of sensitivity of JH signal in B. mori wing disc at the end of the instar larvae [19, 20]. However, the mechanism for the JH insensitiveness in the wing disc is not revealed yet. In this study, during the larva-to-pupa metamorphosis, 20E signal transduction pathway related genes, such as E74, E75 and HR3, were significantly up-regulated in B. mori wing disc. However JH signal transduction pathway related genes, such as Met, SRC and Kr-h1, were significantly down-regulated and their expression levels were relatively low. High expression of BR-C was detected in B. mori wing disc, while the expression of E75A was not detectible before pupation (Tables 1 and 2). These results implied that the loss of sensitivity to JH is probably due to the low expression of JH-related genes during larvae-to-pupae metamorphosis, resulting that the B. mori wing disc can be transformed to pupal wing despite high JH titer in the hemolymph before pupation. However what causes the inactive expression of transmitters of JH signaling during metamorphosis has not been revealed yet.
Out of the 17 up-regulated transcription factors, 10 were significantly up-regulated at the PP stage (Table 3). These genes may be involved in positive regulation of the transition of wing disc to pupal wing. In Drosophila, HR38 was involved in regulating the expression of cuticle protein genes during pupal stage . Cuticle protein genes were actively transcribed in B. mori wing at P0 stage (Figure 5B), so HR38 might be involved in modulating their expression. Previous study indicated that GRF might be a cofactor of βFTZ-F1 and HR39 . βFTZ-F1 was reported to have an essential role in regulating the expression of cuticle protein [51, 52]. In this way, the transcription factor GRF may be positively involved in regulating the expression of cuticle protein, which induces the growth of the wing disc during the larva-to-pupa transformation. Transcription factors HR3 and E75 were involved in 20E signal transduction and they were found to be capable of binding with the response element of many cuticle protein genes [12, 53]. Therefore, they may also be involved in regulating the expression of cuticle proteins in the wing disc during metamorphosis. In Drosophila, gonad, also known as Eagle, is important for the development of nervous system, and egMZ360/egMZ360 (a partial loss-of-function allele) Drosophila showed a held-out wing phenotype and lost the ability of flying [54, 55], implying that Eagle may be involved in the normal development of wing. C/EBP transcription factor was reported to have a role in regulating chorion gene expression in silkmoth, however, its role in wing disc development has not been studied yet . Delilah is a marker gene for epidermal muscle attachment site in Drosophila and may play a role in epidermal muscle attachment site development [57, 58], but its role for B. mori wing development has not been clarified yet.
Out of the 17 transcription factors, except for the PP up-regulated transcription factors, six transcription factors were up-regulated only at the P0 stage (Table 3) which may play a role in pupal wing development in B. mori. For example, ASH2 is required for normal formation of wing scales in B. mori during pupal stage . Transcription factor Pokkuri is also known as Yan and Aop in Drosophila. Yan is regulated by epidermal growth factor (EGF). Ectopic expression of Yan in the A-P boundary of Drosophila wing disc can lead to the mis-expression of Vein, a gene involved in normal development of wing vein and the small wing, and result in the loss of bristles phenotype [59, 60]. This implies that Pokkuri is involved in wing development, but the exact role of Pokkuri in B. mori wing development is still uncovered. SOX-4, a HMG-box transcription factor, is an activator in regulating the expression of neuronal gene and necessary for neuronal maturation . Nervous system can be found in wing vein and the up-regulated SOX-4 at P0 stage may be involved in regulating the neuronal development in B. mori wing vein. Pupal wing of B. mori has reached to the proper size of its adult wing. Krueppel-like factor 10 is an anti-proliferative factor in mammals . It is possible that when B. mori wing disc transforms to pupal wing, cell proliferation is negatively regulated by Krueppel-like factor 10, resulting in the arrest of wing growth when it reaches the proper size. E74 is an early response gene in 20E signaling pathway and is involved in the expression regulation of many cuticle protein genes [12, 52]. Its high level of expression in the wing disc suggests its potential role in regulating the expression of cuticle protein during the wing development. Study in Drosophila suggests that Enhancer of split mbeta [E(spl)mβ] has a function to limit the expression of Rho (Rhomboid), a vein marker gene, in the Drosophila wing vein [59, 63]. In B. mori, wing veins are optimized during pupal stage. Therefore, E(spl)mβ may play similar role as in Drosophila. Functions of the other up-regulated transcription factors remain unknown.
GO classification analyses indicated that the significantly up-regulated transcripts are dominant in structural molecule GO item (Figure 3A and 3B). Actually, cuticle protein genes have high expression levels during the transition of from larval wing disc to pupal wing (Figure 5B and Table 4). Cuticle protein genes were highly expressed at the P0 stage and their RPKM value accounted for about 45.5% of RPKM from all genes (see Additional file 10). These results implied that the large amount of cuticle proteins are synthesized, which is one of the main events related to morphological changes in wing, during larva-to-pupa metamorphosis. In addition, among the 13 WCPs genes identified in B. mori, WCP 1–9 were highly expressed around pupation, while WCP10 and WCP11 have high expression levels before pupation and then significantly down-regulated during pupation [29–31]. In this study, 10 WCPs can be found and their expression patterns are similar to previous studies (Table 4). These results implied that WCP1-9 are the main cuticle proteins for pupal wing, while WCP10 and WCP11 are main cuticle proteins for wing disc in B. mori.
As insects metamorphosize from larvae to pupae, old integument of larvae is degraded and replaced by pupal integument. Wing disc undergoes rapid growth and differentiation during metamorphosis, evaginating to form pupal wing. No cuticle is degraded in B. mori wing disc during the transition from the wing disc to pupal wing. In this study, relevant genes for both chitin degradation and de novo synthesis were significantly up-regulated (Tables 5 and 6 and Figures 6 and 7). This implied that the pupal wing chitin may come from both chitin degradation and de novo synthesis. Chitin deacetylases can catalyze chitin to form chitosan and had a function in organizing tracheal tube length in Drosophila[64, 65]. Chitin deacetylases had higher expression in wing disc at PP stage in B. mori. Whether it involves in extending of tracheal tube in wing is still unknown. In B. mori, there are two kinds of chitin synthase, A and B. While synthase A is responsible for the chitin synthesis in the epidermis, synthase B is responsible for the chitin synthesis in the gut . However it is not clear which synthase is involved in chitin synthesis in insect wing. In this study, synthase A was significantly up-regulated in wing disc during metamorphosis, while the expression of synthase B was almost undetectable (Table 6 and Figure 7D, 7E). This result suggested that the chitin synthesis in pupal wing is majorly performed by synthase A. In Drosophila and Tribolium, Knickkopf (Knk) and Retroactive (Rtv) are essential for chitin synthesis and always along with the expression of chitin synthase A [66, 67]. In Lepidopterans, these two factors had not been studied so far. But Bm_nscaf3003_024 (protein Skeletor, isoforms D/E-like), the homolog of Knk, was found to be significantly up-regulated at PP stage in our data (Additional file 1), whose expression pattern was consistent with that of chitin synthase A (Figure 7D), implying that Protein Skeletor may have similar function as Knk did in Drosophila and Tribolium. The homolog of Rtv was not detected in our data.
In addition, D-fructose can be transformed to D-fructose 6-phosphate by hexokinase, which can be used by chitin synthesis pathway (KEGG pathway map: ko00520). Trehalose, D-glucose and D-fructose were detected in B. mori hemolymph, and their amounts rapidly decreased from L5D6 to PP (Table 7). This suggests that the rapid de novo synthesis of chitin might be taken place in the wing during transition from L5D6 to PP, and both trehalose and D-fructose may be used in de novo synthesis of chitin in the wing of B. mori.
Early gene Fng was also necessary for the B. mori wing development during larva-to-pupa metamorphosis [5, 21, 68]. To investigate other early function genes, expression patterns of these genes were analyzed during metamorphosis in B. mori (Table 8). Results showed that most of them, including Ap and En, which were reported to play important roles in cell fate determination in Drosophila wing disc during early developmental stages, were significantly down-regulated in B. mori wing during metamorphosis, implying that cell fate determination of B. mori wing disc may be decided before pupation.
In this work, total 12,254 assembled transcripts were obtained from the wing disc at L5D6, PP and P0 stages. Totally 5,287 differentially expressed assembled transcripts were identified. From L5D6 to PP (L5D6 vs PP) 2,778 transcripts were differentially expressed and from PP to P0 (PP vs P0) 3,864 transcripts were differentially expressed. More down-regulated transcripts than up-regulated transcripts were identified during the transition from L5D6 to PP and from PP to P0. The transition of wing disc to pupal wing appeared under the control of the 20E signal transduction pathway because high levels of the transcription factors related to 20E signal transduction pathway were found during larva-to-pupa metamorphosis. Seventeen up-regulated transcription factors were screened out, and their functions on B. mori wing development were predicated. High levels of the expression of cuticle proteins and chitin related genes reflect the active chitin metabolism related to the morphological and structural changes during the metamorphosis from the wing disc to pupal wing. Chitin source of pupal wing may come from both degradation of componential chitin and de novo synthesis. Chitin synthase A might be responsible for the chitin synthesis in B. mori pupal wing. De novo synthesis of chitin initiated from D-fructose may exist in B. mori wing.
Experimental animal and RNA isolation
Silkworm, B. mori, were reared with mulberry leaves at 27°C under a 12 h light/12 h dark cycle. Under this condition, the fifth instar stage lasted about 7 days; the wandering stage (W) lasted about 3 days and pupal stage lasted about 10 days. The stage between the end of larval spinning and the onset of pupation was defined as prepupae (PP), which lasted about 12 hours, equal to the period from 2.5-day-old wandering to 3-day-old wandering. P0 was defined as the stage when B. mori just finished pupation. B. mori wing discs were clearly dissected from larvae or pupae at stage of L5D6, PP and P0. Total RNA was isolated by using Rneasy Kit (QIAGEN) according to the manufacturer’s instructions.
RNA-seq sequencing and de novoassembly of sequences
Five μg of total RNA from L5D6, PP and P0 wing was prepared, respectively, for performing RNA-seq sequencing analysis by using Ilumina HiSeq™ 2000 platform. The clean reads of each gene detected in RNA-seq sequencing analysis was assembled by using Trinity software (de novo assembly)  and the assemble results were provided in supplemental files.
Calculation of gene expression level and annotation of predicted proteins
Expression level of genes is calculated by using RPKM (Reads Per kb per Million reads) method , and the formula is:
Here C means the number of reads that uniquely aligned to target gene; N means the total number of reads that uniquely aligned to all genes; L means the number of base on target gene.
Genes were annotated by using Blast2go software (http://www.blast2go.com/b2ghome) according to the translated nucleotide sequences in Blastp.
Isolation and analysis of differentially expressed genes
To compare the differentially expressed genes, the transcripts of PP were compared with those of L5D6 (L5D6 as control), and the transcripts of P0 were compared with those of PP (PP as control). Those genes that had two fold difference between two stages and False Discovery Rate (FDR) was less than 0.001  were identified as differentially expressed genes. The detailed information of their expression patterns can be obtained in the supplemental files (Additional files 1 and 2).
The Cluster 3.0 and Treeview software (http://www.eisenlab.org/eisen/) were used to analyze the expression pattern of differentially expressed genes. Kyoto Encyclopedia of Genes and Genomes (KEGG: http://www.kegg.jp/)  was used to identify transcription factors. BGI Web Gene Ontology Annotation Plotting (WEGO) (http://wego.genomics.org.cn/cgi-bin/wego/index.pl)  were used to classify genes to different GO catalogs and Pearson Chi-Square test was used to justify if the gene numbers of the two input dataset were significantly different.
Quantitative real time -PCR (qRT-PCR)
RNAiso Plus (Takara Dalian, China) was used to isolate total RNA from B. mori wing disc from 3-day-old fifth instar larvae (L5D3), L5D6, 1-day-old wandering (W1), PP, P0 and 3-day-old pupae (P3), respectively. Two μg of total RNAs from each stage was used to synthesize cDNA using reverse transcriptase (Takara) and oligo d (T) primer (Takara). The total volume of real-time quantitative PCR reactions was 20 μl, which contained 10 μl of 2 × SYBR® Premix EX Taq™ (Takara), 0.4 μl of 50 × ROX Reference Dye (Takara) and 0.4 μl of specific primers (10 μM). Detail information of qRT-PCR Primers were provided in Additional file 4. qRT-PCR was performed with an ABI7300 real-time PCR system (Applied Biosystems, Foster City, CA), using the following condition: 95°C for 30 s followed by 40 cycles in 95°C for 5 s, 60°C for 30 s and 72°C for 31 s. The mRNA quantity of each gene was calculated with the 2-△△Ct method and normalized to the abundance of a house-keeping gene, Bmβ-actin (GenBank NO: EU780706.1). The relative mRNA levels of each gene were represented as folds over the expression levels of Bmβ-actin.
Hemolymph collection and sugar detection
The hemolymph of B. mori was collected by cutting its proleg (L5D6 and PP stages) or pupal wing (P0 stage), and quickly transformed to centrifuge tube by pipette. All hemolymph samples were stored at -80°C, until performing the sugar detection. The sugar detection was performed by China National Analytical Center (Guangzhou) by HPLC method.
Availability of supporting data
The assembled unigenes (larger than 200 bp) were deposited in the Transcriptome Shotgun Assembly Sequence Database (TSA) at DDBJ/EMBL/GenBank under the accession GBJR00000000. The version described in this paper is the first version, GBJR01000000.
3-day-old fifth instar larvae
6-day-old fifth instar larvae
Krüppel homolog 1
Steroid receptor coactivator
Epidermal growth factor receptor
Wing cuticle protein
Reads Per kb per Million reads
False Discovery Rate
Kyoto Encyclopedia of Genes and Genomes
BGI Web Gene Ontology Annotation Plotting.
Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, Zhao P, Zha X, Cheng T, Chai C, Pan G, Xu J, Liu C, Lin Y, Qian J, Hou Y, Wu Z, Li G, Pan M, Li C, Shen Y, Lan X, Yuan L, Li T, Xu H, Yang G, Wan Y, Zhu Y, Yu M, Shen W, et al: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306 (5703): 1937-1940.
International Silkworm Genome C: The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem Mol Biol. 2008, 38 (12): 1036-1045. 10.1016/j.ibmb.2008.11.004.
Wang J, Xia Q, He X, Dai M, Ruan J, Chen J, Yu G, Yuan H, Hu Y, Li R, Feng T, Ye C, Lu C, Wang J, Li S, Wong GK, Yang H, Wang J, Xiang Z, Zhou Z, Yu J: SilkDB: a knowledgebase for silkworm biology and genomics. Nucleic Acids Res. 2005, 33 (Database issue): D399-D402.
Shimomura M, Minami H, Suetsugu Y, Ohyanagi H, Satoh C, Antonio B, Nagamura Y, Kadono-Okuda K, Kajiwara H, Sezutsu H, Nagaraju J, Goldsmith MR, Xia Q, Yamamoto K, Mita K: KAIKObase: an integrated silkworm genome database and data mining tool. BMC Genomics. 2009, 10: 486-10.1186/1471-2164-10-486.
Klein T: Wing disc development in the fly: the early stages. Curr Opinion Genetics Dev. 2001, 11 (4): 470-475. 10.1016/S0959-437X(00)00219-7.
Deng H, Zhang J, Li Y, Zheng S, Liu L, Huang L, Xu WH, Palli SR, Feng Q: Homeodomain POU and Abd-A proteins regulate the transcription of pupal genes during metamorphosis of the silkworm, Bombyx mori. Proc Natl Acad Sci U S A. 2012, 109 (31): 12598-12603. 10.1073/pnas.1203149109.
Dubrovsky EB: Hormonal cross talk in insect development. Trends in Endocrinol Metab TEM. 2005, 16 (1): 6-11. 10.1016/j.tem.2004.11.003.
Thummel CS: Molecular mechanisms of developmental timing in C. elegans and Drosophila. Dev Cell. 2001, 1 (4): 453-465. 10.1016/S1534-5807(01)00060-0.
Truman JW, Riddiford LM: The origins of insect metamorphosis. Nature. 1999, 401 (6752): 447-452. 10.1038/46737.
Truman JW, Riddiford LM: Endocrine insights into the evolution of metamorphosis in insects. Annu Rev Entomol. 2002, 47: 467-500. 10.1146/annurev.ento.47.091201.145230.
Riddiford LM, Hiruma K, Zhou X, Nelson CA: Insights into the molecular basis of the hormonal control of molting and metamorphosis from Manduca sexta and Drosophila melanogaster. Insect Biochem Mol Biol. 2003, 33 (12): 1327-1338. 10.1016/j.ibmb.2003.06.001.
Thummel CS: Ecdysone-regulated puff genes 2000. Insect Biochem Mol Biol. 2002, 32 (2): 113-120. 10.1016/S0965-1748(01)00112-6.
Huet F, Ruiz C, Richards G: Sequential gene activation by ecdysone in Drosophila melanogaster: the hierarchical equivalence of early and early late genes. Development. 1995, 121 (4): 1195-1204.
Ashok M, Turner C, Wilson TG: Insect juvenile hormone resistance gene homology with the bHLH-PAS family of transcriptional regulators. Proc Natl Acad Sci U S A. 1998, 95 (6): 2761-2766. 10.1073/pnas.95.6.2761.
Charles JP, Iwema T, Epa VC, Takaki K, Rynes J, Jindra M: Ligand-binding properties of a juvenile hormone receptor, Methoprene-tolerant. Proc Natl Acad Sci U S A. 2011, 108 (52): 21128-21133. 10.1073/pnas.1116123109.
Konopova B, Jindra M: Juvenile hormone resistance gene Methoprene-tolerant controls entry into metamorphosis in the beetle Tribolium castaneum. Proc Natl Acad Sci U S A. 2007, 104 (25): 10488-10493. 10.1073/pnas.0703719104.
Parthasarathy R, Tan A, Palli SR: bHLH-PAS family transcription factor methoprene-tolerant plays a key role in JH action in preventing the premature development of adult structures during larval-pupal metamorphosis. Mech Dev. 2008, 125 (7): 601-616. 10.1016/j.mod.2008.03.004.
Kayukawa T, Minakuchi C, Namiki T, Togawa T, Yoshiyama M, Kamimura M, Mita K, Imanishi S, Kiuchi M, Ishikawa Y, Shinoda T: Transcriptional regulation of juvenile hormone-mediated induction of Kruppel homolog 1, a repressor of insect metamorphosis. Proc Natl Acad Sci U S A. 2012, 109 (29): 11729-11734. 10.1073/pnas.1204951109.
Obara Y, Miyatani M, Ishiguro Y, Hirota K, Koyama T, Izumi S, Iwami M, Sakurai S: Pupal commitment and its hormonal control in wing imaginal discs. J Insect Physiol. 2002, 48 (10): 933-944. 10.1016/S0022-1910(02)00160-9.
Koyama T, Obara Y, Iwami M, Sakurai S: Commencement of pupal commitment in late penultimate instar and its hormonal control in wing imaginal discs of the silkworm, Bombyx mori. J Insect Physiol. 2004, 50 (2–3): 123-133.
Sato K, Matsunaga TM, Futahashi R, Kojima T, Mita K, Banno Y, Fujiwara H: Positional cloning of a Bombyx wingless locus flugellos (fl) reveals a crucial role for fringe that is specific for wing morphogenesis. Genetics. 2008, 179 (2): 875-885. 10.1534/genetics.107.082784.
Kubota K, Goto S, Eto K, Hayashi S: EGF receptor attenuates Dpp signaling and helps to distinguish the wing and leg cell fates in Drosophila. Development. 2000, 127 (17): 3769-3776.
Bray S: Drosophila development: Scalloped and Vestigial take wing. Current Biol CB. 1999, 9 (7): R245-R247. 10.1016/S0960-9822(99)80154-7.
Morata G: How Drosophila appendages develop. Nat Rev Mol Cell Biol. 2001, 2 (2): 89-97. 10.1038/35052047.
Schwank G, Basler K: Regulation of organ growth by morphogen gradients. Cold Spring Harbor Perspect Biol. 2010, 2 (1): a001669-
Swarup S, Verheyen EM: Wnt/Wingless signaling in Drosophila. Cold Spring Harbor Perspect Biol. 2012, 4: 6-
Neto-Silva RM, Wells BS, Johnston LA: Mechanisms of growth and homeostasis in the Drosophila wing. Annu Rev Cell Dev Biol. 2009, 25: 197-220. 10.1146/annurev.cellbio.24.110707.175242.
Merzendorfer H, Zimoch L: Chitin metabolism in insects: structure, function and regulation of chitin synthases and chitinases. J Exp Biol. 2003, 206 (Pt 24): 4393-4412.
Takeda M, Mita K, Quan GX, Shimada T, Okano K, Kanke E, Kawasaki H: Mass isolation of cuticle protein cDNAs from wing discs of Bombyx mori and their characterizations. Insect Biochem Mol Biol. 2001, 31 (10): 1019-1028. 10.1016/S0965-1748(01)00050-9.
Noji T, Ote M, Takeda M, Mita K, Shimada T, Kawasaki H: Isolation and comparison of different ecdysone-responsive cuticle protein genes in wing discs of Bombyx mori. Insect Biochem Mol Biol. 2003, 33 (7): 671-679. 10.1016/S0965-1748(03)00048-1.
Futahashi R, Okamoto S, Kawasaki H, Zhong YS, Iwanaga M, Mita K, Fujiwara H: Genome-wide identification of cuticular protein genes in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2008, 38 (12): 1138-1146. 10.1016/j.ibmb.2008.05.007.
Gan LP, Zhang WY, Niu YS, Xu L, Xi J, Ji MM, Xu SQ: Construction and application of an electronic spatiotemporal expression profile and gene ontology analysis platform based on the EST database of the silkworm, Bombyx mori. J Insect sci. 2010, 10: 114-
Mita K, Morimyo M, Okano K, Koike Y, Nohata J, Kawasaki H, Kadono-Okuda K, Yamamoto K, Suzuki MG, Shimada T, Goldsmith MR, Maeda S: The construction of an EST database for Bombyx mori and its application. Proc Natl Acad Sci U S A. 2003, 100 (24): 14121-14126. 10.1073/pnas.2234984100.
Funaguma S, Hashimoto S, Suzuki Y, Omuro N, Sugano S, Mita K, Katsuma S, Shimada T: SAGE analysis of early oogenesis in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2007, 37 (2): 147-154. 10.1016/j.ibmb.2006.11.001.
Zhang Y, Huang J, Jia S, Liu W, Li M, Wang S, Miao X, Xiao H, Huang Y: SAGE tag based cDNA microarray analysis during larval to pupal development and isolation of novel cDNAs in Bombyx mori. Genomics. 2007, 90 (3): 372-379. 10.1016/j.ygeno.2007.05.005.
Liang J, Zhang L, Xiang Z, He N: Expression profile of cuticular genes of silkworm, Bombyx mori. BMC Genomics. 2010, 11: 173-10.1186/1471-2164-11-173.
Xia Q, Cheng D, Duan J, Wang G, Cheng T, Zha X, Liu C, Zhao P, Dai F, Zhang Z, He N, Zhang L, Xiang Z: Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori. Genome Biol. 2007, 8 (8): R162-10.1186/gb-2007-8-8-r162.
Li Y, Wang G, Tian J, Liu H, Yang H, Yi Y, Wang J, Shi X, Jiang F, Yao B, Zhang Z: Transcriptome analysis of the silkworm (Bombyx mori) by high-throughput RNA sequencing. PLoS One. 2012, 7 (8): e43713-10.1371/journal.pone.0043713.
Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.
Curtis K, Youngquist ST: Part 21: categoric analysis: Pearson chi-square test. Air Med J. 2013, 32 (4): 179-180. 10.1016/j.amj.2013.04.007.
Matsuoka T, Fujiwara H: Expression of ecdysteroid-regulated genes is reduced specifically in the wing discs of the wing-deficient mutant (fl) of Bombyx mori. Dev Genes Evol. 2000, 210 (3): 120-128. 10.1007/s004270050018.
Ote M, Mita K, Kawasaki H, Seki M, Nohata J, Kobayashi M, Shimada T: Microarray analysis of gene expression profiles in wing discs of Bombyx mori during pupal ecdysis. Insect Biochem Mol Biol. 2004, 34 (8): 775-784. 10.1016/j.ibmb.2004.04.002.
Zhou Q, Yu L, Shen X, Li Y, Xu W, Yi Y, Zhang Z: Homology of dipteran bristles and lepidopteran scales: requirement for the Bombyx mori achaete-scute homologue ASH2. Genetics. 2009, 183 (2): 619-627. 10.1534/genetics.109.102848. 611SI-613SI
Besson MT, Chareyre P, Deleage G, Demont J, Quennedey B, Courrent A, Fourche J, Bosquet G: Isolation and characterization of Urbain, a 20-hydroxyecdysone-inducible gene expressed during morphogenesis of Bombyx mori wing imaginal discs. Rouxs Arch Dev Biol. 1996, 205 (7–8): 333-343.
Abdou M, Peng C, Huang J, Zyaan O, Wang S, Li S, Wang J: Wnt signaling cross-talks with JH signaling by suppressing Met and gce expression. PLoS One. 2011, 6 (11): e26772-10.1371/journal.pone.0026772.
Muramatsu D, Kinjoh T, Shinoda T, Hiruma K: The role of 20-hydroxyecdysone and juvenile hormone in pupal commitment of the epidermis of the silkworm, Bombyx mori. Mech Dev. 2008, 125 (5–6): 411-420.
Kozlova T, Lam G, Thummel CS: Drosophila DHR38 nuclear receptor is required for adult cuticle integrity at eclosion. Dev Dyn Off Publ Am Assoc Anatomists. 2009, 238 (3): 701-707.
Charles JP, Shinoda T, Chinzei Y: Characterization and DNA-binding properties of GRF, a novel monomeric binding orphan receptor related to GCNF and betaFTZ-F1. Eur J Biochem/FEBS. 1999, 266 (1): 181-190. 10.1046/j.1432-1327.1999.00842.x.
Wang HB, Nita M, Iwanaga M, Kawasaki H: betaFTZ-F1 and Broad-Complex positively regulate the transcription of the wing cuticle protein gene, BMWCP5, in wing discs of Bombyx mori. Insect Biochem Mol Biol. 2009, 39 (9): 624-633. 10.1016/j.ibmb.2009.06.007.
Deng H, Zheng S, Yang X, Liu L, Feng Q: Transcription factors BmPOUM2 and BmbetaFTZ-F1 are involved in regulation of the expression of the wing cuticle protein gene BmWCP4 in the silkworm, Bombyx mori. Insect Mol Biol. 2011, 20 (1): 45-60. 10.1111/j.1365-2583.2010.01041.x.
Ali MS, Iwanaga M, Kawasaki H: Ecdysone-responsive transcriptional regulation determines the temporal expression of cuticular protein genes in wing discs of Bombyx mori. Gene. 2013, 512 (2): 337-347. 10.1016/j.gene.2012.09.126.
Lundell MJ, Hirsh J: eagle is required for the specification of serotonin neurons and other neuroblast 7–3 progeny in the Drosophila CNS. Development. 1998, 125 (3): 463-472.
Gritzan U, Weiss C, Brennecke J, Bohmann D: Transrepression of AP-1 by nuclear receptors in Drosophila. Mech Dev. 2002, 115 (1–2): 91-100.
Sourmeli S, Papantonis A, Lecanidou R: A novel role for the Bombyx Slbo homologue, BmC/EBP, in insect choriogenesis. Biochem Biophys Res Commun. 2005, 337 (2): 713-719. 10.1016/j.bbrc.2005.09.103.
Armand P, Knapp AC, Hirsch AJ, Wieschaus EF, Cole MD: A novel basic helix-loop-helix protein is expressed in muscle attachment sites of the Drosophila epidermis. Mol Cell Biol. 1994, 14 (6): 4145-4154.
Staudt N, Molitor A, Somogyi K, Mata J, Curado S, Eulenberg K, Meise M, Siegmund T, Hader T, Hilfiker A, Bronner G, Ephrussi A, Rorth P, Cohen SM, Fellert S, Chung HR, Piepenburg O, Schafer U, Jackle H, Vorbruggen G: Gain-of-function screen for genes that affect Drosophila muscle pattern formation. PLoS Genet. 2005, 1 (4): e55-10.1371/journal.pgen.0010055.
Biehs B, Sturtevant MA, Bier E: Boundaries in the Drosophila wing imaginal disc organize vein-specific genetic programs. Development. 1998, 125 (21): 4245-4257.
Wessells RJ, Grumbling G, Donaldson T, Wang SH, Simcox A: Tissue-specific regulation of vein/EGF receptor signaling in Drosophila. Dev Biol. 1999, 216 (1): 243-259. 10.1006/dbio.1999.9459.
Bergsland M, Werme M, Malewicz M, Perlmann T, Muhr J: The establishment of neuronal properties is controlled by Sox4 and Sox11. Genes Dev. 2006, 20 (24): 3475-3486. 10.1101/gad.403406.
Dang DT, Pevsner J, Yang VW: The biology of the mammalian Kruppel-like family of transcription factors. Int J Biochem Cell Biol. 2000, 32 (11–12): 1103-1121.
de Celis JF, Bray S, Garcia-Bellido A: Notch signalling regulates veinlet expression and establishes boundaries between veins and interveins in the Drosophila wing. Development. 1997, 124 (10): 1919-1928.
Wang S, Jayaram SA, Hemphala J, Senti KA, Tsarouhas V, Jin H, Samakovlis C: Septate-junction-dependent luminal deposition of chitin deacetylases restricts tube elongation in the Drosophila trachea. Curr Biol CB. 2006, 16 (2): 180-185. 10.1016/j.cub.2005.11.074.
Zhao Y, Park RD, Muzzarelli RA: Chitin deacetylases: properties and applications. Marine Drugs. 2010, 8 (1): 24-46. 10.3390/md8010024.
Moussian B, Tang E, Tonning A, Helms S, Schwarz H, Nusslein-Volhard C, Uv AE: Drosophila Knickkopf and Retroactive are needed for epithelial tube growth and cuticle differentiation through their specific requirement for chitin filament organization. Development. 2006, 133 (1): 163-171. 10.1242/dev.02177.
Chaudhari SS, Arakane Y, Specht CA, Moussian B, Boyle DL, Park Y, Kramer KJ, Beeman RW, Muthukrishnan S: Knickkopf protein protects and organizes chitin in the newly synthesized insect exoskeleton. Proc Natl Acad Sci U S A. 2011, 108 (41): 17028-17033. 10.1073/pnas.1112288108.
Milan M, Perez L, Cohen SM: Short-range cell interactions and cell survival in the Drosophila wing. Dev Cell. 2002, 2 (6): 797-805. 10.1016/S1534-5807(02)00169-7.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.
Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36 (Database issue): D480-D484.
Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (Web Server issue): W293-W297.
This research was supported by grants from National Basic Research Program of China (2012CB114602) and National Natural Science Foundation of China (31172265; 31330071).
The authors declare that they have no competing interests.
QLF and LL conceived and designed the study. HMD collected and prepared RNA samples for RNA-seq. JO performed RNA extraction and qRT-PCR, JO, LL and LHH contributed to bioinformatic analysis. JO, QLF, LL and SCZ contributed to writing the manuscript and data analysis. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Differentially expressed transcripts in L5D6-VS-PP. The criteria applied for significance difference were FDR ≤ 0.001, and estimated absolute |log2Ratio (PP/L5D6)| ≥ 1. L5D6- and PP-RPKM: reads per kb per million reads for each transcript in stages L5D6 and PP, respectively. Log2 Ratio (PP/L5D6): the ratio between the RPKM in PP and the RPKM in L5D6. Up-Down (PP/L5D6) and FDR of each gene are also shown. GO Component, GO Function and Go Process: ontology information of Cellular Components, Molecular Function and Biological Processes of Gene-corresponding GO terms. "-": no hit. Blast nr: identification of homologues in GenBank; E-values are also shown. (XLS 980 KB)
Additional file 2: Table S2: Differentially expressed transcripts in PP-VS-P0. Description of data is similar to Table S1. Log2 Ratio (P0/PP): the ratio between the RPKM in P0 and the RPKM in PP. (XLS 1 MB)
Additional file 10: Table S4: All transcripts at P0 stage. All transcripts at stage P0 were listed. Description of data is similar to Table S1. (XLS 3 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ou, J., Deng, HM., Zheng, SC. et al. Transcriptomic analysis of developmental features of Bombyx mori wing disc during metamorphosis. BMC Genomics 15, 820 (2014). https://doi.org/10.1186/1471-2164-15-820
- B. mori
- Wing disc