Varroa destructor is a devastating ectoparasitic mite which attacks Honeybee, Apis mellifera. V. destructor feeds on honeybee hemolymph, and often harbors small RNA viruses such as the deformed wing virus to transmit these viruses in the infested bee hive. To survey the genes of V. destructor, total RNA was subjected to high-throughput transcriptome sequencing to construct in silico cDNA library by using the Illumina HiSeq 2000 platform. Total of 2×107,748,792 paired-end short reads were obtained and quality filtered reads were subjected to Trinity de novo assembler followed by TransDecoder, and CD-HIT program to make a V. destructor reference cDNA library containing 28,023 of clustered contigs with protein coding capacity. These cDNA sequences will help us to understand the molecular biology of V. destructor.

The Asian honeybee, Apis cerana, is a native honeybee species in Korea which is important in agriculture for pollination and honey production. For better understanding of the physiology of A. cerana, high-throughput Illumina transcriptome sequencing was performed to analyze the gene expression profiles of queen, worker and larva. A total of 219,799,682 clean reads corresponding to 22.2Gb of nucleotide sequences was obtained from the whole body total RNA samples. The Apis mellifera reference mRNA sequence database was used to measure the gene expression level with Bowtie2 and eXpress software and the Illumina short reads were mapped to 11,459 out of 11,736 A. mellifera reference genes. Total of 9,221 genes with FPKM value greater than 5 of each sample group were subjected to evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) with BLASTX for gene ontology analysis. The differential gene expression between queen and worker, and worker and larva were analyzed to screen the overexpressed genes in each sample group. in the queen and worker sample group, total of 1,766 genes were differentially expressed with 887 and 879 genes overexpressed over two folds in queen and worker, respectively. In the worker and larva sample group, total of 1,410 genes were differentially expressed with 1,009 and 401 genes overexpressed over two folds in worker and larva, respectively.

Sacbrood virus (SBV) is one of the most fatal pathogens against Asian honeybee, Apis cerana. This virus cause failure of the insect larvae to pupate and death of the adult insects. This study has analyzed the host genes affected by viral infection, by comparing the expression level of host transcripts infected with or without SBV. As a first step, we sequenced the cDNA libraries of Asian honeybee by using illumina RNA sequencing. The sequences were de novo assembled to acquire honeybee transcriptome sequences. The transcriptome was annotated by the sequence comparison to known protein sequences by BLASTX and evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) database with functional categories and description. By mapping the RNA-seq data to de novo assembled transcripts, we characterized the differentially expressed transcripts between SBV-infected and non-infected Asian honeybee.