Characterization of the transcriptome of the Asian gypsy moth Lymantria dispar identifies numerous transcripts associated with insecticide resistance

•We identified 62,063 unigenes with N50 of 993 bp in the transcriptome of Lymantria dispar.•GO and COG analysis identified 7,309 and 8,079 sequences, respectively.•We obtained a large number of transcripts associated with known insecticide targets.•Twenty-seven transcripts showed potential roles in L. dispar development.

Although the Asian gypsy moth Lymantria dispar causes extensive forest damage worldwide, little is known regarding the genes involved in its development or response to insecticides. Accordingly, characterization of the transcriptome of L. dispar larvae would promote the development of toxicological methods for its control. RNA-seq analysis of L. dispar larvae messenger RNA (mRNA) generated 62,063 unigenes with N50 of 993 bp, from which 23,975 unique sequences (E-value < 10−5) were identified using a BLASTx search of the NCBI non-redundant (nr) database. Using functional classification in the Gene Ontology (GO) and Clusters of Orthologous Groups (COG) databases, 7,309 indentified sequences were categorized into 51 functional groups and 8,079 sequences were categorized into 25 functional groups, respectively. Moreover, we identified a large number of transcripts encoding known insecticide targets, or proteins involved in the metabolism of insecticides. Reads per kilobase of unigene length per million mapped reads (RPKM) analysis identified 39 high abundance transcripts, of which 27 exhibited significantly altered expression patterns across the egg, larvae, pupae, male and female adult stages. Our study provides the most comprehensive transcriptomic sequence resource for L. dispar, which will form the basis for future identification of candidate insecticide resistance genes in L. dispar.

Graphical AbstractFigure optionsDownload full-size imageDownload as PowerPoint slide