Talytic activity (3,869 unigenes) was prominently represented (Figure 5). The Cluster of Orthologous Groups (COG) is a database where the orthologous gene products were classified. All unigenes were aligned to the COG database to predict and classify possible functions [26]. Out of 30,427 nr hits, 
9,009 sequences were assigned to the COG classifications (Figure 6). Among the 25 COG function categories, the cluster for General function prediction only (3,519, 20.90 ) represented the largest group, followed by replication, recombination and repair (1,359, 8.07 ) (Figure 6).Results and Discussion Illumina Paired-end Sequencing and de novo AssembleTotal RNA was extracted from the worker heads of the different colonies. Using Illumina paired-end sequencing technology, a total of 57,271,634 raw sequencing reads were generated from a 200 bp insert library. An assembler, Trinity was employed for de novo assembly [21]. After stringent quality check and data cleaning, approximately 54 million high-quality reads were obtained with 98.09 Q20 bases (base quality more than 20). Based on the high quality reads, a total of 221,728 contigs were assembled with an average length of 302 bp. The size distribution of these contigs is shown in Figure 1. Then the reads were mapped back to contigs, with paired-end reads we were able to detect contigs from the same transcript as well as the distances between these contigs. After clustering these unigenes using TGICL software [22], contigs can finally generate 116,885 unigenes with 9,040 distinct clusters and 107,845 distinct singletons (Table 1). The length of assembled unigenes ranged from 150 to 17,355 bp. There were 83,002 unigenes (71.01 ) with length varying from 150 to 500 bp, 26,916 unigenes (23.03 ) in the length range of 501 to 1500 bp, and 6967 unigenes (5.96 ) with 23977191 length more than 1500 bp. The size distribution of these unigenes is showed in Figure 2.Functional Classification by KEGGThe Kyoto Encyclopedia of Genes and BTZ043 Genomes (KEGG) Pathway database records the networks of molecular interactions in the cells, and their variants of them specific to particular organisms. In order to identify the biological Lixisenatide web pathways involved, the assembled unigenes were annotated with corresponding Enzyme commission (EC) numbers from BLASTX alignments against the KEGG database [27]. 23727046 Firstly, based on a comparison against the KEGG database using BLASTX with an E-value cutoff of ,1025, out of the 116,885 unigenes, 19,611 (16.78 ) had significant matches in the database and were assigned to 242 KEGG pathways. The pathways most represented by unique sequences were metabolic pathways (2,282 members), Huntington’s disease (683 members), purine metabolism (661 members), RNA transport (629 members), and regulation of actin cytoskeleton (306 members). Taken together, 30,643 unique sequence-based annotations had BLAST scores exceeding our threshold (#1e-5) in nr, Swiss-Prot and KEGG databases (Figure 7A). The Venn diagram (Figure 7B) shows that an additional 3 unigenes were annotated by domainbased alignments. Overall, 30,646 unique sequence-based or domain-based annotations using the four selected public databases were assigned to O. formosanus unigenes (26.2 ). Among them, 8,458 unigenes had hits in all four public databases with relatively defined functional annotations of the assembled unigenes (TableFunctional Annotation by Searching Against Public DatabasesFor validation and annotation of assembled unigenes, sequence simil.Talytic activity (3,869 unigenes) was prominently represented (Figure 5). The Cluster of Orthologous Groups (COG) is a database where the orthologous gene products were classified. All unigenes were aligned to the COG database to predict and classify possible functions [26]. Out of 30,427 nr hits, 9,009 sequences were assigned to the COG classifications (Figure 6). Among the 25 COG function categories, the cluster for General function prediction only (3,519, 20.90 ) represented the largest group, followed by replication, recombination and repair (1,359, 8.07 ) (Figure 6).Results and Discussion Illumina Paired-end Sequencing and de novo AssembleTotal RNA was extracted from the worker heads of the different colonies. Using Illumina paired-end sequencing technology, a total of 57,271,634 raw sequencing reads were generated from a 200 bp insert library. An assembler, Trinity was employed for de novo assembly [21]. After stringent quality check and data cleaning, approximately 54 million high-quality reads were obtained with 98.09 Q20 bases (base quality more than 20). Based on the high quality reads, a total of 221,728 contigs were assembled with an average length of 302 bp. The size distribution of these contigs is shown in Figure 1. Then the reads were mapped back to contigs, 
with paired-end reads we were able to detect contigs from the same transcript as well as the distances between these contigs. After clustering these unigenes using TGICL software [22], contigs can finally generate 116,885 unigenes with 9,040 distinct clusters and 107,845 distinct singletons (Table 1). The length of assembled unigenes ranged from 150 to 17,355 bp. There were 83,002 unigenes (71.01 ) with length varying from 150 to 500 bp, 26,916 unigenes (23.03 ) in the length range of 501 to 1500 bp, and 6967 unigenes (5.96 ) with 23977191 length more than 1500 bp. The size distribution of these unigenes is showed in Figure 2.Functional Classification by KEGGThe Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway database records the networks of molecular interactions in the cells, and their variants of them specific to particular organisms. In order to identify the biological pathways involved, the assembled unigenes were annotated with corresponding Enzyme commission (EC) numbers from BLASTX alignments against the KEGG database [27]. 23727046 Firstly, based on a comparison against the KEGG database using BLASTX with an E-value cutoff of ,1025, out of the 116,885 unigenes, 19,611 (16.78 ) had significant matches in the database and were assigned to 242 KEGG pathways. The pathways most represented by unique sequences were metabolic pathways (2,282 members), Huntington’s disease (683 members), purine metabolism (661 members), RNA transport (629 members), and regulation of actin cytoskeleton (306 members). Taken together, 30,643 unique sequence-based annotations had BLAST scores exceeding our threshold (#1e-5) in nr, Swiss-Prot and KEGG databases (Figure 7A). The Venn diagram (Figure 7B) shows that an additional 3 unigenes were annotated by domainbased alignments. Overall, 30,646 unique sequence-based or domain-based annotations using the four selected public databases were assigned to O. formosanus unigenes (26.2 ). Among them, 8,458 unigenes had hits in all four public databases with relatively defined functional annotations of the assembled unigenes (TableFunctional Annotation by Searching Against Public DatabasesFor validation and annotation of assembled unigenes, sequence simil.
