Transcriptomic analysis of differentially expressed genes in the Ras1CA-overexpressed and wildtype posterior silk glands
BMC Genomics volume 15, Article number: 182 (2014)
Using the piggyBac-mediated GAL4/UAS transgenic system established in the silkworm, Bombyx mori, we have previously reported that overexpression of the Ras1CA oncogene specifically in the posterior silk gland (PSG) improved cell growth, fibroin synthesis, and thus silk yield. However, the detailed molecular mechanism remains to be fully elucidated. To achieve this goal, Illumina sequencing was used in the present study to compare the transcriptomes of the Ras1CA-overexpressed and wildtype PSGs.
The transcriptomic sequencing results in 56 million reads following filtering steps. Most of the reads (~70%) are successfully mapped to the Bombyx genome. The mapped reads are situated within at least 9,133 predicted genes, covering 62.46% genes of the Bombyx genome. GO annotation shows that 2512 of the 2,636 differentially expressed genes (DEGs) are mostly distributed in metabolic process, cell and cell part, and binding, and KEGG annotation shows that 1,941 DEGs are mapped into 277 pathways. Importantly, Ras1CA overexpression in the PSG upregulated many DEGs distributed in “pathways in cancer”, “insulin signaling pathway”, and “MAPK signaling pathway” as well as “purine metabolism” and “pyrimidine metabolism”. Transcriptional regulation of these DEGs was verified by quantitative real-time PCR. Moreover, injection of small-molecule chemical inhibitors of the Ras1 downstream effectors into the Ras1CA-overexpressed silkworms revealed that both Raf-MAPK and PI3K-TORC1 pathways are required for the Ras1-induced DEG expression.
The transcriptomic analysis illustrates that, apart from phosphorylational regulation, Ras1 activates its downstream Raf-MAPK and PI3K-TORC1 pathways at the transcriptional level. Meanwhile, Ras1 increases DNA content and induces endoreplication, at least in part, by upregulating genes in “nucleotide metabolism” and “cell cycle”. This study provides further insights into the molecular mechanism of how Ras1CA overexpression in the PSG improves silk yield.
As a traditional agricultural industry, sericulture is economically important to China and several other countries. The domesticated silkworm, Bombyx mori, is the most important insect species for sericulture. In the past, sericulture has been greatly advanced by applying the hybrid breeding technique to Bombyx. However, it has reached a plateau during the last decades, mostly due to the inherent threshold of this technique. To break through the bottleneck of silk production, new breeding techniques, such as the molecular breeding technique, should be developed. It has been hypothesized that fibroin production in the Bombyx posterior silk gland (PSG) is directly proportional to silk yield and determined by its gland size and protein synthesis , making it possible to improve silk yield by genetic manipulation of the PSG . On the other hand, Bombyx has been used as a model lepidopteran insect for a long time [3, 4]. In terms of protein synthesis, its PSG is one of the most efficient organs in animals. Therefore, studying the molecular mechanism controlling fibroin synthesis in the Bombyx PSG is of great value for both applied and basic research.
Using the piggyBac-mediated GAL4/UAS transgenic system established in Bombyx[5–7], we specifically overexpressed a constitutively active form of Ras1 (Ras1CA) in the PSG. In the transgenic silkworm, Fil-GAL4 > UAS-Ras1CA (Fil > Ras1CA), Ras1CA overexpression increases gland size and protein synthesis in the PSG, leading to silk yield improvement by 60% . This study not only provides an application prospect to silk yield improvement in sericulture, but also supports the previous hypothesis that fibroin production is determined by gland size and protein synthesis in the PSG . It is certain that Ras1 and its downstream Raf-MAPK and PI3K-TORC1 pathways play critical roles in regulating fibroin production , while the detailed molecular mechanism remains to be fully elucidated. The completed Bombyx genome sequence [8, 9] makes it possible to use functional genomics, such as proteomics and transcriptomics, to achieve the above goal. Using 2D-DIGE-MS/MS analysis, we previously compared the proteomic profiles of the Ras1CA-overexpressed and wild type (WT) PSGs. Further studies revealed that, via the downstream Raf-MAPK and PI3K-TORC1 pathways, Ras1CA upregulates bcpi, which inhibits cathepsin activity thus preventing PSG destruction during metamorphosis . Transcriptomics could be an alternative approach for analyzing how Ras1CA overexpression in the PSG improves fibroin production.
In terms of transcriptomic tools, a whole-genome microarray containing 22,987 oligonucleotides of 70-mers that cover the presently known and predicted genes in the Bombyx genome was designed on the basis of the whole-genome sequences . This microarray has been used to survey the silkworm gene expression patterns in multiple tissues, at different developmental stages, and under various conditions or treatments [11–15]. RNA-Seq (also known as Next Generation Sequencing), including Roche/454 pyrosequencing, Illumina-Solexa sequencing, and Applied Biosystems SOLiD sequencing, has led to a revolution in genomics and provided cheaper and faster delivery of sequencing information . The Illumina-Solexa sequencing strategy was adopted for the sequencing of 40 Bombyx genomes from 29 phenotypically and geographically diverse domesticated silkworm lines and 11 wild silkworms from various mulberry fields in China. This comprehensive study constructs a genome-wide genetic variation map which shed light on the history of silkworm domestication . RNA-Seq also led to the identification of new exons, novel genes, alternative splicing genes, and trans-splicing events in Bombyx[17, 18].
In this study, Illumina-Solexa sequencing revealed 2,636 differentially expressed genes (DEGs) in the Ras1CA-overexpressed and WT PSGs. Confirmed by quantitative real-time PCR (qPCR), the transcriptomic analysis shows that Ras1 increases gland size, protein synthesis, and DNA content in the PSG at the transcriptional level.
Identification of DEGs using RNA-seq
For better understanding the molecular mechanism how Ras1CA overexpression increases fibroin production in the Bombyx PSG, we compared the transcriptomes of the Ras1CA-overexpressed and WT PSGs at the early wandering stage. The RNA-seq raw data were deposited to NCBI SRA with the accession number SRP026709 (http://www.ncbi.nlm.nih.gov/sra/?term=SRP026709). The accession numbers for the two WT PSG RNA-seq libraries are SRX320122 and SRX320124, and those for the two Ras1CA-overexpressed PSG RNA-seq libraries are SRX320125 and SRX320126. The RNA-seq raw data of SRX320122 and SRX320124 were combined for raw data processing, so were SRX320125 and SRX320126. Using the pair-end Illumina-Solexa sequencing strategy, we obtained 69,662,027 raw reads, containing about 6.9 gigabases with an average read length of 101 nucleotides. The raw data was preprocessed, and 13,238,008 (19%) low quality reads were removed. The remaining 56,424,019 reads, with an average length of 90.71 nucleotides (Table 1), were used to map the Bombyx genome (release_2.0)  using TopHat .
In total, 39,967,028 reads were mapped to silkworm genome database (silkgenome.fa) and silkworm gene database (silkworm_glean.gff), with approximately 2.2 and 1.8 million reads from the Ras1CA-overexpressed and WT PSGs, respectively. The completed Bombyx genome contains 14,623 unigenes [8, 9]. Using the Fragments Per kb Per Million Reads (FPKM) method , we have found that 9,133 unigenes are expressed, with 6,962 and 6,429 unigenes in the Ras1CA-overexpressed and WT PSGs, respectively. The mapping coverage is 62.46% genes of the silkworm genome, showing a high confidence of RNA-seq in this study (Table 2).
Importantly, RNA-Seq revealed 2,636 DEGs, with 1708 upregulated and 938 downregulated genes in the Ras1CA-overexpressed PSG compared to the WT PSG (Figure 1). Next, the DEGs were subjected to functional annotation and qPCR verification.
Functional annotation of DEGs
Gene ontology (GO) assignments were used to classify the functions of DEGs revealed by transcriptomic analysis. The DEGs were termed by GO ontology in three categories, namely biological process, cellular component, and molecular function. In total, 2512 DEGs (95.3% of 2636) were annotated in 60 GO functional groups (Figure 2). In the category of biological process, over 43% of DEGs were distributed in “metabolic process”, nearly 40% DEGs were classified into “cellular process”, while very few DEGs were found in “growth”, “carbon utilization”, “viral reproduction”, “rhythmic process”, and “locomotion” (Figure 2A). In the category of cellular component, over 35% of DEGs were distributed in “cell” and “cell part”, however, very small numbers of DEGs were found in “virion part”, “synapse”, and “extracellular matrix part” (Figure 2B). In the category of molecular function, the terms “binding” and “catalytic activity” enriched 32% and 30% of DEGs, respectively. By contrast, few DEGs were distributed in “channel regulator activity”, “receptor regulator activity”, and “metallochaperone activity” (Figure 2C).
To identify the biological pathways that are active in the Ras1CA-overexpressed PSG, we mapped the DEGs to the reference canonical pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG) . Among the 2636 DEGs, 1,941 sequences predicted to encode enzymes with enzyme commission (EC) numbers were mapped into 277 KEGG pathways. The top 25 KEGG pathways are shown in Table 3. The top 5 KEGG pathways with most representations by the DEGs were “purine metabolism”, “spliceosome”, “pathways in cancer”, “RNA transport”, and “HTLV-1 infection”, implying that these diverse metabolic processes are active in the Ras1CA-overexpressed PSG. “Purine metabolism” is essential, providing components of the nucleotides, DNA and RNA, and the energy currency of the cell, ATP . In parallel, “Pyrimidine metabolism” also ranks on the top 6 KEGG pathway (Table 3). In the Ras1CA-overexpressed PSG, the DNA content is nearly doubled and endoreplication is enhanced in comparison with the WT PSG, indicating that more active nucleotide metabolism is required . “Spliceosome” is complex molecular machinery involved in removal of introns from mRNA precursors . A lot of DEGs in this pathway suggest occurrence of high frequency of mRNA splicing for translation in the Ras1CA-overexpressed PSG. Interestingly, a large number of DEGs are found distributed in “pathways in cancer”. As an important oncogene, Ras plays important roles in both normal development and aberrant biological processes, such as tumorigenesis and developmental disorders . It is likely that, apart from the phosphorylational regulation, Ras could also activate “pathways in cancer” at the transcriptional level. “RNA transport” from the nucleus to the cytoplasm is fundamental for gene expression regulation. Most eukaryotic RNAs are produced in the nucleus by RNA polymerase I, II, or III. The RNA molecules undergo a variety of posttranscriptional processing events, and a majority of them are localized to their functional sites in the cell . In accordance with this result, ribosome biogenesis for mRNA translation was stimulated in the Ras1CA-overexpressed PSG . Unexpectedly, some DEGs were annotated into “HTLV-1 infection pathway”, but the silkworm cannot be affected by this human virus. Our preliminary data shows that virus-resistance is enhanced when Ras1CA is globally overexpressed in the transgenic silkworm, Actin3-GAL4 > UAS-Ras1CA, suggesting that the DEGs in “HTLV-1 infection pathway” might play a role in antivirus in Bombyx.
In short, these annotations provide a valuable insight into the specific processes, functions, and pathways and facilitate the identification of DEGs resulted from the Ras1CA-overexpressed PSG. Next, the DEGs in several important and top KEGG pathways were chosen for qPCR verification (Table 3).
qPCR verification of DEGs in “pathways in cancer”, “insulin signaling pathway”, and “MAPK signaling pathway”
In our previous study, we have determined that Ras1CA overexpression in the PSG increases Ras activity, resulting in phosphorylation of the Ras downstream effector proteins, Raf and PI3K110, which in turn activate the Raf-MAPK and PI3K-TORC1 pathways, respectively . Surprisingly, among all the 277 KEGG pathways, a large number of DEGs are distributed in “pathways in cancer”, “insulin signaling”, and “MAPK signaling pathway” (Additional file 1: Figure S1, Additional file 2: Figure S2 and Additional file 3: Figure S3), which rank on the top 3, 9, and 15 KEGG pathways containing 49, 35, and 32 DEGs, respectively (Table 4).
Since these three KEGG pathways are highly related to each other, we analyzed their common DEGs and verified them by qPCR (Figure 3). First, 4 DEGs were found in all the three pathways, including crk, mek, jnk, and erk. qPCR data were in agreement with the transcriptomic results, showing that all these 4 genes are upregulated in the Ras1CA-overexpressed PSG (Figure 4A). Second, 5 DEGs were observed in both “pathways in cancer” and “insulin signaling pathway”, including pi3kl, pi3ks, cbl1, cbl2, and cbl3 (Figure 4B). The transcriptomic results and qPCR data showed that pi3ks, cbl1, cbl2, and cbl3 are all upregulated by Ras1CA, whereas the qPCR data of pi3kl did not match its transcriptomic result suggesting it might also be upregulated by Ras1CA. Third, there were 3 DEGs in both “pathways in cancer” and “MAPK signaling pathway”, including fgfr1, pkc, and fgfr2. The transcriptomic results and qPCR data showed the fgfr1 and pkc are upregulated and downregulated by Ras1CA, respectively, while fgfr2 might be downregulated by Ras1CA (Figure 4C). Finally, pka, a single DEG in both “insulin signaling pathway” and “MAPK signaling pathway”, was downregulated by Ras1CA (Figure 4D).
We randomly selected some DEGs in each individual pathway for qPCR verification (Figure 5). In “pathways in cancer”, both the transcriptomic results and qPCR data revealed that 7 DEGs, including cdk4/6, ptenl, fu, cdk2, ral, faks, and elongb, are all upregulated by Ras1CA, while fakl is downregulated. Although the transcriptomic results suggested that ptens and smo are upregulated by Ras1CA, qPCR data showed that they might be downregulated (Figure 5A). In “insulin signaling pathway”, 7 DEGs, including rheb, fasl, fass, apkc, apkc, eif4e, and ampk, are all upregulated by Ras1CA, while phk is downregulated, and shc might be upregulated (Figure 5B). In “MAPK signaling pathway”, p38 and mapkapk are upregulated by Ras1CA, while hsp70s and tak are downregulated, and daxx might be downregulated (Figure 5C).
qPCR verification of DEGs in “purine metabolism” and “pyrimidine metabolism”
Interestingly, abundant DEGs were annotated in two major nucleotide metabolism pathways, “purine metabolism” and “pyrimidine metabolism” (Additional file 4: Figure S4 and Additional file 5: Figure S5), ranking on the top 1 and 6 KEGG pathways which contain 67 and 42 DEGs, respectively (Table 3).
We thus verified the above-mentioned hypothesis that Ras might activate nucleotide metabolism by qPCR verification of some randomly selected DEGs in both “purine metabolism” and “pyrimidine metabolism” (Figure 6). We first analyzed 10 of the 28 common DEGs in both pathways, including pole4, pole2, rpb5, rpc10, rpb4, rpc37, apf, itpa, rpc25, and nt5e. Different from the transcriptomic results, qPCR data suggest that rpc25 and nt5e might be upregulated rather than downregulated by Ras1 (Figure 6A). We then analyzed 5 of the 39 DEGs only in “purine metabolism”, including adk, allc, prps, pde, and gart. Although gart expression was inconsistent between its transcriptomic and qPCR data, all the other DEGs are upregulated (Figure 6B). In addition, urh1, one of DEGs only in “pyrimidine metabolism” is also upregulated (Figure 6C). Taken together, Ras1CA overexpression in the PSG upregulates most, if not all, DEGs in “purine metabolism” and “pyrimine metabolism”.
Identification of Ras1 downstream signals in regulating DEGs by small-molecule inhibitor treatments
We next investigated which downstream pathway(s) that Ras1 utilizes to regulate DEGs in the Bombyx PSG by injecting small-molecule inhibitors of the Ras downstream effectors (Rafi, Raf inhibitor; LY294002, PI3K inhibitor; rapamycin, TORC1 inhibitor; 15 μg/larva) into the Ras1CA-overexpressed silkworm larvae . Some Ras1CA-upregulated DEGs, which are consistent in both transcriptomic results and qPCR data, were chosen for inhibitor treatment experiments by qPCR analysis to examine their expression levels.
First, we detected the common DEGs annotated in “pathways in cancer”, “insulin signaling pathway”, and “MAPK signaling pathway” by qPCR (Figure 7A-C). The mRNA levels of mek, erk, and jnk distributed in all the three pathways were decreased to 10-20% by Rafi and 20-40% by LY294002, whereas rapamycin treatment showed weaker inhibitory effects (40-60%) (Figure 7A). For pi3ks, cbl2, and cbl3, the three DEGs presented in both “pathways in cancer” and “insulin signaling pathway”, LY294002 and rapamycin showed the strongest and weakest inhibitory effects, respectively (Figure 7B). By contrast, rapamycin strongly inhibited expression of fgfr1, the DEGs distributed in “pathways in cancer” and “MAPK signaling” (Figure 7C).
Second, we detected the individual DEGs annotated in “pathways in cancer” (Figure 7D), “insulin signaling pathway” (Figure 7E), and “MAPK signaling pathway” (Figure 7F) by qPCR. For most of the DEGs, Raf inhibitor exhibited the strongest inhibitory effects, while rapamycin showed little to no inhibitory effects.
Third, we detected the DEGs annotated in “purine metabolism” and “pyrimidine metabolism” (Figure 8). For DEGs in both pathways, LY294002 exhibited the strongest inhibitory effects (Figure 8A). For the two DEGs only in “purine metabolism”, Raf inhibitor and LY294002 showed the strongest inhibitory effects on pde and allc, respectively (Figure 8B).
In summary, inhibitors of the Ras downstream effectors showed inhibitory efforts on different DEGs to varying degrees indicating that both Raf-MAPK and PI3K-TORC1 pathways are involved in the transcriptional regulation of those DEGs. Interestingly, similar results were observed in mammalian cells in which Ras is overexpressed or transformed [27, 28].
Ras1 transcriptionally activates its downstream Raf-MAPK and PI3K-TORC1 pathways
On a genome-wide scale, the identification of Ras-responsive genes has become feasible using different transcriptomic tools. For example, subtractive suppression hybridization was performed in immortalized, non-tumorigenic rat embryo fibroblasts and in Ras-transformed cells. The results have shown that many DEGs are involved in almost all aspects of cellular growth control and cell survival . A microarray was conducted in RasCA-transformed mouse embryonic fibroblasts, showing that many genes encoding cell growth-related proteins are upregulated . The results of Ras-induced gene expression profiling studies based on subtractive suppression hybridization and microarrays were extensively summarized , revealing that Kruppel-like factor 5, the CD44 antigen, and members of the epidermal growth factor (EGR)-family are common Ras downstream effectors . Using Illumina-Solexa sequencing to analyze DEGs in the PSG, here we found that many Ras1-induced genes are distributed in “pathways in cancer”, “insulin signaling”, and “MAPK signaling pathway”. The transcriptomic analysis illustrates that, apart from phosphorylational regulation, Ras1 can also activate its downstream Raf-MAPK and PI3K-TORC1 pathways at the transcriptional level (Figure 3, 4 and 5, Additional file 1: Figure S1, Additional file 2: Figure S2 and Additional 3: Figure S3, Table 4). To our knowledge, this is the first report that Ras1 can transcriptionally activate its downstream pathways at a global level.
Ras1 transcriptionally activates genes involved in nucleotide metabolism and cell cycle for increasing DNA content and inducing endoreplication
Earlier studies have shown that a lot of key enzymes of nucleotide metabolism and DNA biosynthesis, such as CTP synthetase, thymidylate synthase, dihydrofolate reductase, IMP dehydrogenase, ribonucleotide reductase, DNA polymerase, and DNA methyltransferase, are markedly upregulated in certain tumor cells, which supports the excessive proliferation of transformed cells . The microarray conducted in RasCA-transformed mouse embryonic fibroblasts revealed that many genes encoding DNA-associated proteins (involved in DNA replication and DNA repair) are upregulated as well . Interestingly, a microarray analysis of Ras-overexpressed hemocytes in the fruitfly, Drosophila melanogaster, showed that a large number of genes that are functionally important in cell cycle regulation and DNA replication were upregulated . We have previously shown that compared to the WT PSG, total DNA content is nearly doubled in the Ras1CA-overexpressed PSG . Moreover, in comparison with the WT PSG, BrdU incorporation in the Ras1CA-overexpressed PSG is much higher indicating enhanced endoreplicative cycles . In this study, we found many Ras1-induced genes are enriched in “purine metabolism”, “pyrimidine metabolism”, and “cell cycle”, which ranks top 1, 6, and 7, respectively (Figure 6, Table 3). Therefore, it is likely that Ras1 transcriptionally activates genes involved in nucleotide metabolism and cell cycle for increasing DNA content and inducing endoreplication in the PSG.
About 46 years before, it has been hypothesized that fibroin production in the Bombyx PSG is directly proportional to silk yield and determined by its gland size and protein synthesis . Based on this hypothesis, we have generated a transgenic silkworm, Fil > Ras1CA, for improving silk yield. Importantly, overexpression of Ras1CA increases gland size and protein synthesis in the PSG, improving fibroin production and silk yield by 60%. At the molecular level, we have determined that Ras activation enhances phosphorylation levels of Ras downstream effector proteins, Raf and PI3K110, and thus activates its downstream Raf-MAPK and PI3K-TORC1 pathways . To better understand the molecular mechanisms how Ras1CA overexpression in the PSG improves fibroin production and silk yield, we performed both proteomics and transcriptomics. Unfortunately, in spite of the discovery that Ras1CA upregulates bcpi to inhibit cathepsin activity and thus to prevent PSG destruction, we were not able to better understand how Ras1CA improves fibroin production and silk yield using proteomics .
The transcriptomic results of the Ras1CA-overexpressed PSG presented here underlie a wide array of DEGs in many KEGG pathways. Importantly, we have discovered that a large number of DEGs in “pathways in cancer”, “insulin signaling”, and “MAPK signaling pathway” are upregulated by Ras1CA overexpression in the PSG. Combined with our previous findings  and the present studies (Figure 9), we conclude that Ras1 activates its downstream Raf-MAPK and PI3K-TORC1 pathways at both phosphorylational and transcriptional levels. Moreover, we find that Ras1 upregulates genes in “nucleotide metabolism” and “cell cycle” for increasing DNA content and inducing endoreplication (Figure 9). This study has advanced our knowledge on how Ras1CA overexpression in the PSG improves fibroin production and silk production.
The Bombyx strain, Dazao, was reared on fresh mulberry leaves in the laboratory at 25°C under 14 h light/10 h dark cycles. The transgenic silkworm Fil > Ras1CA was obtained as described previously by crossing Fil-GAL4 with UAS-Ras1CA. The transgenic silkworms were reared under the same condition as the WT silkworms.
Dissection of the PSGs
The silkworm PSGs were dissected from the Ras1CA-overexpressed and WT silkworms at the early wandering stage, when silkworms just begin to spin [2, 10]. The PSGs were used for Illumina-Solexa sequencing and qPCR. Throughout the paper, all qPCR experiments were performed in 3 biological duplicates.
For Illumina sequencing, total RNA from the Ras1CA-overexpressed PSG or the WT PSG was isolated with TRIzol (Invitrogen, Carlsbad, CA, USA). To remove any residual DNA, samples were pretreated with RNase-free DNase I (New England BioLabs, Ipswich, MA, USA) for 30 minutes at 37°C. RNA quality was first verified using a 2100 Bioanalyzer RNA Nanochip (Agilent, Santa Clara, CA) with RNA Integrity Number (RIN) value over 8.5. RNA was then quantified using NanoDrop ND-1000 Spectrophotometer (NanoDrop, Wilmington, DE).
Library preparation and Illumina sequencing
The RNA-seq libraries were prepared using Illumina’skit following the manufacturer’s protocol (Illumina, San Diego, CA). Approximately 20 μg of total RNA from the Ras1CA-overexpressed PSG or the WT PSG was used to isolate mRNA using Sera-mag Magnetic Oligo (dT) Beads (Illumina). To avoid cDNA synthesizing bias by priming, the purified mRNA was fragmented into small pieces (100–400 bp) using divalent cations at 94°C for 5 minutes. The double-stranded cDNA was first synthesized using the SuperScript Double-Stranded cDNA Synthesis kit (Invitrogen, Camarillo, CA) with random hexamer (N6) primers (Illumina). Then the synthesized cDNA was subjected to end-repair and phosphorylation using T4 DNA polymerase, Klenow DNA polymerase, and T4 PNK. These repaired cDNA fragments were 3' adenylated using Klenow Exo- (3' to 5' exo minus, Illumina).Illumina paired-end adapters were ligated to the ends of these 3'-adenylated cDNA fragments.The ligated cDNA was then enriched with 15 rounds of PCR amplification using PCR Primer PE 1.0 and PCR Primer PE 2.0 (Illumina) with Phusion DNA Polymerase.The libraries were sequenced using Illumina Highseq 2000 platform according to the manufacturer’s instructions. Illumina sequencing was performed at Suzhou Encode Genomics Biotechnology Co-Ltd (Encode Genomics; Suzhou, China).
Raw data preprocess
Preprocessing was carried out with a stringent filtering process. First, we removed reads that do not pass the built-in Illumina's software Failed-Chastity filter according to the relation "failed-chastity < = 1", using a chastity threshold of 0.6, on the first 25 cycles. Second, we discarded all reads with adaptor contamination. Third, we ruled out low-quality reads containing ambiguous sequences "N". Finally, the reads with more than 10% Q < 20 bases were also removed.
Genome mapping and abundance analysis
Quality-filtered reads were then aligned to the Bombyx genome (release_2.0, ftp://silkdb.org/pub/current/Genome/ by TopHat (version 2.0.4)  with the parameters “--bowtie1 -r 0 --mate-std-dev 50 -N 3--solexa1.3-quals” (insert size is set as 0). The resulting alignment data from Tophat were then fed to an assembler Cufflinks (version 2.0.1) to assemble aligned RNA-Seq reads into silkworm genome database (silkworm_glean.gff) and silkworm gene database (silkgenome.fa). Unigene abundances were measured by Fragments per kb of exon per million fragments mapped (FPKM) using the formula FPKM = (1,000,000*C)/ (N*L*1,000).
The DEGs in Ras1CA-overexpressed and WT PSGs were functional annotated by GO annotation and KEGG annotation. For GO annotation, the DEGs were first blasted against uniprot knowledgebase (including Swiss-Prot and TrEMBL) (UniProtKB; http://www.uniprot.org) to get uniprot IDs. Then the uniprot IDs were assigned to GO terms at three basic catergories including molecular function, biological process, and cellular component. For KEGG annotation, DEGs were functionally annotated with KAAS (KEGG Automatic Annotation Server) by BLAST comparisons against the manually curated KEGG GENES database. The result contains KO (KEGG Orthology) assignments and automatically generated KEGG pathways.
Total RNA of the Ras1CA-overexpressed PSG or the WT PSG was extracted using TRIzol (Invitrogen). qPCR was performed as previously described [2, 33]. The primers used in this paper are listed in Additional file 6: Table S1 (Supporting Information).
Chemical inhibitor treatment
Small-molecule chemical inhibitors of the Ras downstream effectors (Rafi, Raf inhibitor, Santa Cruz; LY294002, PI3K inhibitor, Sigma; rapamycin, TORC1 inhibitor, Sigma; 15 μg/larva) were injected into the Fil > Ras1CA larvae at the EW stage. At 24 h after the first injection, the PSGs were dissected for qPCR as previously described .
Experimental data were analyzed with the Student’s t-test and ANOVA. t-test: *, p < 0.05; **, p < 0.01. ANOVA: the bars labeled with different lowercase letters are significantly different (p < 0.05). Throughout the paper, values are represented as the mean ± standard deviation of at least 3 independent experiments.
Availability of supporting data
The RNA-seq raw data were deposited to NCBI SRA with the accession number SRP026709 (http://www.ncbi.nlm.nih.gov/sra/?term=SRP026709).
Tashiro Y, Morimoto T, Matsuura S, Nagata S: Studies on the posterior silk gland of the silkworm Bombyx mori I Growth of posterior silk gland cells and biosynthesis of fibroin during the fifth larval instar. J Cell Biol. 1968, 38: 574-588. 10.1083/jcb.38.3.574.
Ma L, Xu H, Zhu J, Ma S, Liu Y, Jiang RJ, Xia Q, Li S: Ras1CA overexpression in the posterior silk gland improves silk yield. Cell Res. 2011, 21: 934-943. 10.1038/cr.2011.36.
Goldsmith MR, Shimada T, Abe H: The genetics and genomics of the silkworm Bombyx mori. Annu Rev Entomol. 2005, 50: 71-100. 10.1146/annurev.ento.50.071803.130456.
Xia QY, Li S, Feng QL: Advances in Silkworm Studies Accelerated by the Genome Sequencing of Bombyx mori. Annu Rev Entomol. 2014, 59: 513-536. 10.1146/annurev-ento-011613-161940.
Tamura T, Thibert C, Royer C, Kanda T, Abraham E, Kamba M, Komoto N, Thomas JL, Mauchamp B, Chavancy G, Shirk P, Fraser M, Prudhomme JC, Couble P: Germline transformation of the silkworm Bombyx mori L. using a piggyBac transposon-derived vector. Nat Biotechnol. 2000, 18: 81-84. 10.1038/71978.
Imamura M, Nakai J, Inoue S, Quan GX, Kanda T, Tamura T: Targeted gene expression using the GAL4/UAS system in the silkworm Bombyx mori. Genetics. 2003, 165: 1329-1340.
Tan AJ, Tanaka H, Tamura T, Shiotsuki T: Precocious metamorphosis in transgenic silkworms overexpressing juvenile hormone esterase. Proc Natl Acad Sci U S A. 2005, 102: 11751-11756. 10.1073/pnas.0500954102.
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: 1937-1940.
Xia Q, Guo Y, Zhang Z, Li D, Xuan Z, Li Z, Dai F, Li Y, Cheng D, Li R, Cheng T, Jiang T, Becquet C, Xu X, Liu C, Zha X, Fan W, Lin Y, Shen Y, Jiang L, Jensen J, Hellmann I, Tang S, Zhao P, Xu H, Yu C, Zhang G, Li J, Gao J, Liu S: Complete resequencing of 40 genomes reveals domestication events and genes in silkworm (Bombyx). Science. 2009, 326: 433-436. 10.1126/science.1176620.
Ma L, Liu S, Shi M, Chen XX, Li S: Ras1CA-Upregulated bcpi inhibits cathepsin activity to prevent tissue destruction of the bombyx posterior silk gland. J Proteome Res. 2013, 12: 1924-1934. 10.1021/pr400005g.
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: R162-10.1186/gb-2007-8-8-r162.
Tian L, Guo E, Diao Y, Zhou S, Peng Q, Cao Y, Ling E, Li S: Genome-wide regulation of innate immunity by juvenile hormone and 20-hydroxyecdysone in the Bombyx fat body. BMC Genomics. 2010, 11: 549-10.1186/1471-2164-11-549.
Tian L, Guo E, Wang S, Liu S, Jiang RJ, Cao Y, Ling E, Li S: Developmental regulation of glycolysis by 20-hydroxyecdysone and juvenile hormone in fat body tissues of the silkworm, Bombyx mori. J Mol Cell Biol. 2010, 2: 255-263. 10.1093/jmcb/mjq020.
Gan L, Liu X, Xiang Z, He N: Microarray-based gene expression profiles of silkworm brains. BMC Neurosci. 2011, 12: 8-10.1186/1471-2202-12-8.
Xu Q, Lu A, Xiao G, Yang B, Zhang J, Li X, Guan J, Shao Q, Beemtsen BT, Zhang P, Wang C, Ling E: Transcriptional profiling of midgut immunity response and degeneration in the wandering silkworm. Bombyx mori. PLoS One. 2012, 7: e43769-10.1371/journal.pone.0043769.
Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92: 255-264. 10.1016/j.ygeno.2008.07.001.
Shao W, Zhao QY, Wang XY, Xu XY, Tang Q, Li M, Li X, Xu YZ: Alternative splicing and trans-splicing events revealed by analysis of the Bombyx mori transcriptome. RNA. 2012, 18: 1395-1407. 10.1261/rna.029751.111.
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: e43713-10.1371/journal.pone.0043713.
Duan J, Li R, Cheng D, Fan W, Zha X, Cheng T, Wu Y, Wang J, Mita K, Xiang Z, Xia Q: SilkDB v2.0: a platform for silkworm (Bombyx mori) genome biology. Nucleic Acids Res. 2010, 38: D453-456. 10.1093/nar/gkp801.
Trapnell C, Pachter L, Salzberg SL: TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
Di Virgilio F: Purines, purinergic receptors, and cancer. Cancer Res. 2012, 72: 5441-7. 10.1158/0008-5472.CAN-12-1600.
Bonnal S, Vigevani L, Valcárcel J: The spliceosome as a target of novel antitumour drugs. Nat Rev Drug Discov. 2012, 11: 847-859. 10.1038/nrd3823.
Karnoub AE, Weinberg RA: Ras oncogenes: split personalities. Nat Rev Mol Cell Biol. 2008, 9: 517-531. 10.1038/nrm2438.
Lécuyer E, Yoshida H, Parthasarathy N, Alm C, Babak T, Cerovina T, Hughes TR, Tomancak P, Krause HM: Global analysis of mRNA localization reveals a prominent role in organizing cellular architecture and function. Cell. 2007, 131: 174-187. 10.1016/j.cell.2007.08.003.
Schäfer R, Tchernitsa OI, Sers C: Global effects of Ras signaling on the genetic program in mammalian cells. RAS family GTPases. Edited by: Der CJ. 2006, Dordrecht: Springer, 169-198.
Györffy B, Schäfer R: Biomarkers downstream of RAS: a search for robust transcriptional targets. Curr Cancer Drug Targets. 2010, 10: 858-868. 10.2174/156800910793357916.
Zuber J, Tchernitsa OI, Hinzmann B, Schmitz A-C, Grips M, Hellriegel M, Sers C, Rosenthal A, Schafer R: A genome-wide survey of Ras transformation targets. Nat Genet. 2000, 24: 144-152. 10.1038/72799.
Vasseur S, Malicet C, Calvo EL, Labrie C, Berthezene P, Dagorn JC, Iovanna JL: Gene expression profiling by DNA microarray analysis in mouse embryonic fibroblasts transformed by rasV12 mutated protein and the E1A oncogene. Mol Cancer. 2003, 19: 19-
Hatse S, De Clercq E, Balzarini : Role of antimetabolites of purine and pyrimidine nucleotide metabolism in tumor cell differentiation. J Biochem Pharmacol. 1999, 58: 539-555. 10.1016/S0006-2952(99)00035-0.
Asha H, Nagy I, Kovacs G, Stetson D, Ando I, Dearolf CR: Analysis of Ras-induced overproliferation in Drosophila hemocytes. Genetics. 2003, 163: 203-215.
Liu Y, Zhou S, Ma L, Tian L, Wang S, Sheng Z, Jiang R-J, Bendena WG, Li S: Transcriptional regulation of the insulin signaling pathway genes by starvation and 20-hydroxyecdysone in the Bombyx fat body. J Insect Physiol. 2010, 56: 1436-1444. 10.1016/j.jinsphys.2010.02.011.
This study was supported the by Natural Science Foundation of China (31201747 to LM, 31125025 to SL), the 973 program (2012CB114605 to SL, 2012CB316501 to XL), the Shanghai Municipal Education Commission (12ZR1452500 to LM; 13JC1409100 to SL), Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences (2012KIP305 to LM).
The authors declare that they have no competing interests.
LM (the first author) performed most of the experiments and analyzed the data. QM helped with qRT-PCR and inhibitor treatment. XL, LC and KL helped performing some experiments and writing the paper. SL designed the experiments, analyzed the data, wrote the paper, and coordinated the whole study. All authors approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: DEGs distributed in “pathways in cancer” (KEGG) signaling network. DEGs are marked with red square. (EPS 2 MB)
Additional file 2: Figure S2: DEGs distributed in “insulin signaling pathway” (KEGG) signaling network. DEGs are marked with red square. (EPS 1 MB)
Additional file 3: Figure S3: DEGs distributed in “MAPK signaling pathway” (KEGG) signaling network. DEGs are marked with red square. (EPS 2 MB)
Additional file 4: Figure S4: DEGs distributed in “purine metabolism” (KEGG) signaling network. DEGs are marked with red square. (EPS 7 MB)
Additional file 5: Figure S5: ADEGs distributed in “pyrimidine metabolism” (KEGG) signaling network. DEGs are marked with red square. (EPS 8 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ma, L., Ma, Q., Li, X. et al. Transcriptomic analysis of differentially expressed genes in the Ras1CA-overexpressed and wildtype posterior silk glands. BMC Genomics 15, 182 (2014). https://doi.org/10.1186/1471-2164-15-182
- MAPK Signaling Pathway
- Insulin Signaling Pathway
- Purine Metabolism
- Pyrimidine Metabolism
- Transgenic Silkworm