Transcriptomic analysis of the housefly (Musca domestica) larva using massively parallel pyrosequencing

To explore the transcriptome of Musca domestica larvae and to identify unique sequences, we used massively parallel pyrosequencing on the Roche 454-FLX platform to generate a substantial EST dataset of this fly. As a result, we obtained a total of 249,555 ESTs with an average read length of 373 bp. These reads were assembled into 13,206 contigs and 20,556 singletons. Using BlastX searches of the Swissprot and Nr databases, we were able to identify 4,814 contigs and 8,166 singletons as unique sequences. Subsequently, the annotated sequences were subjected to GO analysis and the search results showed a majority of the query sequences were assignable to certain gene ontology terms. In addition, functional classification and pathway assignment were performed by KEGG and 2,164 unique sequences were mapped into 184 KEGG pathways in total. As the first attempt on large-scale RNA sequencing of M. domestica, this general picture of the transcriptome can establish a fundamental resource for further research on functional genomics.