Identification of three annexin IX isoforms generated by alternative splicing of the carboxyl-terminal exon in silkworm, Bombyx mori

Improving Silkworm Genome Annotation Using a Proteogenomics Approach

The silkworm genome has been deeply sequenced and assembled, but accurate genome annotation, which is important for modern biological research, remains far from complete. To improve silkworm genome annotation, we carried out a proteogenomics analysis using 9.8 million mass spectra collected from different tissues and developmental stages of the silkworm. The results confirmed the translational products of 4307 existing gene models and identified 1701 novel genome search-specific peptides (GSSPs). Using these GSSPs, 74 novel gene-coding sequences were identified, and 121 existing gene models were corrected. We also identified 1182 novel junction peptides based on an exon-skipping database that resulted in the identification of 973 alternative splicing sites. Furthermore, we performed RNA-seq analysis to improve silkworm genome annotation at the transcriptional level. A total of 1704 new transcripts and 1136 new exons were identified, 2581 untranslated regions (UTRs) were revised, and 1301 alternative splicing (AS) genes were identified. The transcriptomics results were integrated with the proteomics data to further complement and verify the new annotations. In addition, 14 incorrect genes and 10 skipped exons were verified using the two analysis methods. Altogether, we identified 1838 new transcripts and 1593 AS genes and revised 5074 existing genes using proteogenomics and transcriptome analyses. Data are available via ProteomeXchange with identifier PXD009672. The large-scale proteogenomics and transcriptome analyses in this study will greatly improve silkworm genome annotation and contribute to future studies.

Transcriptome analysis of the silkworm (Bombyx mori) by high-throughput RNA sequencing

The domestic silkworm, Bombyx mori, is a model insect with important economic value for silk production that also acts as a bioreactor for biomaterial production. The functional complexity of the silkworm transcriptome has not yet been fully elucidated, although genomic sequencing and other tools have been widely used in its study. We explored the transcriptome of silkworm at different developmental stages using high-throughput paired-end RNA sequencing. A total of about 3.3 gigabases (Gb) of sequence was obtained, representing about a 7-fold coverage of the B. mori genome. From the reads that were mapped to the genome sequence; 23,461 transcripts were obtained, 5,428 of them were novel. Of the 14,623 predicted protein-coding genes in the silkworm genome database, 11,884 of them were found to be expressed in the silkworm transcriptome, giving a coverage of 81.3%. A total of 13,195 new exons were detected, of which, 5,911 were found in the annotated genes in the Silkworm Genome Database (SilkDB). An analysis of alternative splicing in the transcriptome revealed that 3,247 genes had undergone alternative splicing. To help with the data analysis, a transcriptome database that integrates our transcriptome data with the silkworm genome data was constructed and is publicly available at http://124.17.27.136/gbrowse2/. To our knowledge, this is the first study to elucidate the silkworm transcriptome using high-throughput RNA sequencing technology. Our data indicate that the transcriptome of silkworm is much more complex than previously anticipated. This work provides tools and resources for the identification of new functional elements and paves the way for future functional genomics studies.

Detection and analysis of alternative splicing in the silkworm by aligning expressed sequence tags with the genomic sequence

We identified 277 alternative splice forms in silkworm genes based on aligning expressed sequence tags with genomic sequences, using a transcript assembly program. A large fraction (74%) of these alternative splices are located in protein-coding regions and alter protein products, whereas only 26% are in untranslated regions. From the alternative splices located in protein-coding regions, some (43%) affect protein domains that bind various biological molecules. The vast majority of the detected alternative forms in this study appear to be novel, and potentially affect biologically meaningful control of function in silkworm genes. Our results indicate that alternative splicing in silkworm largely produces protein diversity and functional diversity, and is a widely used mechanism for regulating gene expression.