Abstract

Insects are one of the most successful animal groups on earth. Some insects, such as the silkworm and honeybee, are beneficial to humans, whereas others are notorious pests of crops. At present, the genomes of 38 insects have been sequenced and made publically available. In addition, the transcriptomes of dozens of insects have been sequenced. As gene data rapidly accumulate, constructing the pathway of molecular interactions becomes increasingly important for entomological research. Here, we developed an improved tool, iPathCons, for knowledge-based construction of pathways from the transcriptomes or the official gene sets of genomes. Considering the high evolution diversity in insects, iPathCons uses a voting system for Kyoto Encyclopedia of Genes and Genomes Orthology assignment. Both stand-alone software and a web server of iPathCons are provided. Using iPathCons, we constructed the pathways of molecular interactions of 52 insects, including 37 genome-sequenced and 15 transcriptome-sequenced ones. These pathways are available in the iPathDB, which provides searches, web server, data downloads, etc. This database will be highly useful for the insect research community.

Database URL:http://ento.njau.edu.cn/ipath/

Introduction

Insects are one of the most successful animal groups on earth. They comprise more than a million species, representing about half of all known living organisms. Some insects, such as silkworm (Bombyx mori) and honeybee (Apis mellifera) are beneficial to humans by producing valuable products and/or services (silk, honey, pollination). In contrast, other species damage crops by feeding on leaves or fruits, causing huge economic losses.

As the sequencing cost has dramatically declined in recent decade, gene sequences data have accumulated rapidly in insects. The genome sequences of 38 insect species have been reported, including 12 species of Drosophila (1, 2), seven kinds of ants (3–8), three wasps (9), three mosquitoes (10–12), two butterflies (13, 14), the human body louse Pediculus humanus humanus (15), the kissing bug Rhodnius prolixus (16), the tsetse fly Glossina morsitans (16), a tick Ixodes scapularis (16), the honeybee A mellifera (17), the silkworm B. mori (18), stick insect Timema cristinae (19) and several agricultural insect pests, such as the red flour beetle Tribolium castaneum (20), pea aphid Acyrthosiphon pisum (21), diamondback moth Plutella xylostella (22) and locust Locusta migratoria (23). In addition, dozens of insects have been sequenced for their transcriptome (the SRA database, September 2014).

Constructing the pathway of molecular interaction from the insect genomes or transcriptomes is important for gene function analysis. Large-scale gene expression analysis is an efficient and widely used technique in molecular biology experiment. However, selecting the right candidate genes for experiment validation is still a challenge. One solution is to find differently expressed genes in a related pathway. To construct pathways, several knowledge-based methods were developed, such as PANTHER (24), Gene Ontology (25), Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO) (26), Reactome (27) and PharmGKB (28). Although these methods can be applied to the insect gene data, insect pathway construction is still a difficult work. First, most insects have a heterozygous genome, reducing the quality of genome assembly and annotation and increasing the difficulty in pathway construction. Second, compared with mammals and other groups of animals, insects are species-rich and have high evolution diversity. Here, we developed an improved tool for insect pathway construction and built an insect pathway database, which should be helpful to the entomological community.

Data resources

Official gene sets of 37 insect species

The genome sequences of 37 insects were downloaded from the NCBI (National Center for Biotechnology Information) genome database or in species-specific databases, such as FlyBase (29), AphidBase (30), BeeBase, ButterflyBase (31), BeetleBase (32), DBM-DB (33), MonarchBase (34), NasoniaBase, HessianFlyBase, ManducaBase, Ant genome (35), VectorBase (36), SilkDB (37) and ChiloDB. When we prepared this article, the genome of locust was published (23). We did not include locust in this work. We downloaded the official gene sets (OGS) of 37 insects (Supplementary Table S1), including Aedes aegypti (v2.2), Anopheles gambiae (v3.8), Anopheles darlingi (v2.2), Anopheles stephensi (v2.2), Culex quinquefasciatus (v1.4), G. morsitans (v1.3), Lutzomyia longipalpis (v1.1), Pediculus humanus corporis (v1.3), Phlebotomus papatasi (v1.1), Rh. prolixus (v1.1), Acromyrmex echinatior (v3.8), Atta cephalotes (v1.2), Camponotus floridanus (v3.3), Harpegnathos saltator (v3.3), Linepithema humile (v1.2), Pogonomyrmex barbatus (v1.2), Solenopsis invicta (v2.2.3), A. mellifera (v3.2), Nasonia vitripennis (v1.2), Drosophila ananassae (v1.3), Drosophila erecta (v1.3), Drosophila grimshawi (v1.3), Drosophila melanogaster (v5.54), Drosophila mojavensis (v1.3), Drosophila persimilis (v1.3), Drosophila pseudoobscura (v3.1), Drosophila sechellia (v1.3), Drosophila simulans (v1.4), Drosophila virilis (v1.2), Drosophila willistoni (v1.3), Drosophila yakuba (v1.3), Danaus plexippus (v2.0), Heliconius melpomene (v1.1), B. mori (v2.0), T. castaneum (v3.0), Ac. pisum (v2.1b) and Ti. cristinae (v1.0).

Insect transcriptomes of 15 insects

We downloaded the RNA-seq raw data of 15 insects from Sequence Read Archive database (SRA, http://www.ncbi.nlm.nih.gov/sra/). The genomic data of these species are not available at present. They are A. cerana, Lucilia sericata, Rhagoletis pomonella, Aedes albopictus, Galleria mellonella, Chilo suppressalis, Spodoptera exigua, Manduca sexta, Melitaea cinxia, P. xylostella, Zygaena filipendulae, Dendroctonus ponderosae, Nilaparvata lugens, Oncopeltus fasciatus and Bemisia tabaci. The SRA accession numbers are given in Supplementary Table S2. The transcriptome of C. quinquefasciatus was also downloaded and used to estimate the precision and coverage of the iPathCons.

The SRA database only provides the sequences of raw reads. For most insects, the assembled transcriptomes are not available. Therefore, we assembled these 15 transcriptomes de novo. The statistic of assembled transcriptome was presented in the Supplementary Table S3. First, the raw data were cleaned by removing adaptor sequences, empty reads, and low-quality reads that contain N or whose average nucleotides quality is less than 15. Second, we merged the raw reads from different samples of same species to obtain contigs as many as possible. Third, Trinity was used to assemble Illumina Solexa raw reads with default parameters (38). The insects included A. cerana, L. sericata, Ch. suppressalis, P. xylostella, N. lugens, Be. tabaci and C. quinquefasciatus. The Newbler was used to assemble Roche/454 raw data with default parameters, including Sp. exigua, R. pomonella, Ae. albopictus, De. ponderosae, O. fasciatus and M. sexta. For the EST (Expressed Sequence Tag) data of Ga. mellonella, Me. cinxia, and Z. filipendulae, we assembled the transcripts using the Cap3 software (39). Finally, the assembled transcripts were annotated using BLASTX (Basic Local Alignment Search Tool) against the NCBI nr database.

KEGG data

The KEGG database provides the most commonly used resources for pathway analysis (40). KEGG contains the genes of 21 insects (Table 1), which were downloaded from the NCBI RefSeq database (41). Among these 21 insects, 14 were from Diptera, including 11 Drosophila species and three mosquitoes. The KEGG Markup Language (KGML) is a format of the KEGG pathway maps, which can be used to draw KEGG pathways and to model gene and chemical networks. Given that KEGG requires a subscription to its FTP server, we downloaded the KGML files from their individual web download pages.

Table 1.

Data statistic of iPathDB

OrderSpeciesNucleotide or protein sequencesAnnotated sequencesKO termPathway
ColeopteraDe. ponderosae¤25 36071661997261
T. castaneum*16 73320932093256
DipteraAedes aegypti*15 78620972097257
Ae. albopictus¤48 42215 7942013253
Anopheles darlingi11 43043091849187
Anopheles gambiae*12 81121302130258
Anopheles stephensi23 28749801867189
C. quinquefasciatus*18 95520682068259
D. ananassae*15 14320872087257
D. erecta*15 12320892088256
D. grimshawi*15 06620952094257
D. melanogaster*29 45321592156258
D. mojavensis*14 66520822082257
D. persimilis*16 94820372036254
D. pseudoobscura*16 92920722072256
D. sechellia*16 54420732073251
D. simulans*15 48619841983254
D. virilis*14 56720902089257
D. willistoni*15 59620792079257
D. yakuba*16 15821022101257
G. morsitans12 45747261658182
L. sericata¤146 25045 7682272263
Lu. longipalpis10 11039461483176
Ph. papatasi11 16443881506178
R. pomonella¤13 07131571032223
HemipteraA. pisum*36 19820462046252
Be. tabaci¤54 86071612078255
N. lugens¤36 74894752163258
O. fasciatus¤4863922520171
Rh. prolixus15 44145791580178
HymenopteraA. cerana cerana¤50 37311 1432150260
A. mellifera*15 31919811981253
At. cephalotes18 09354241957187
Acromyrmex echinatior17 27856511957187
C. floridanus17 06453651959187
H. saltator18 56458801964187
Li. humile16 11652571966187
Na. vitripennis*18 90019681968256
Po. barbatus17 18952481963187
So. invicta16 52251041793187
LepidopteraB. mori*14 62321382467296
Ch. suppressalis¤37 04070052029256
D. plexippus15 13051111648184
Ga. mellonella¤40 03516 3141675255
He. melpomene12 82944931580181
M. sexta¤67 9021213744189
Me. cinxia¤195 08676211598249
P. xylostella¤456 02693 9412330266
Sp. exigua¤153 61931 0711957260
Z. filipendulae¤42 73511 7091861254
PhthirapteraPe. humanus corporis*10 77320872087256
PhasmatodeaTi. cristinae17 093156 5362073260
OrderSpeciesNucleotide or protein sequencesAnnotated sequencesKO termPathway
ColeopteraDe. ponderosae¤25 36071661997261
T. castaneum*16 73320932093256
DipteraAedes aegypti*15 78620972097257
Ae. albopictus¤48 42215 7942013253
Anopheles darlingi11 43043091849187
Anopheles gambiae*12 81121302130258
Anopheles stephensi23 28749801867189
C. quinquefasciatus*18 95520682068259
D. ananassae*15 14320872087257
D. erecta*15 12320892088256
D. grimshawi*15 06620952094257
D. melanogaster*29 45321592156258
D. mojavensis*14 66520822082257
D. persimilis*16 94820372036254
D. pseudoobscura*16 92920722072256
D. sechellia*16 54420732073251
D. simulans*15 48619841983254
D. virilis*14 56720902089257
D. willistoni*15 59620792079257
D. yakuba*16 15821022101257
G. morsitans12 45747261658182
L. sericata¤146 25045 7682272263
Lu. longipalpis10 11039461483176
Ph. papatasi11 16443881506178
R. pomonella¤13 07131571032223
HemipteraA. pisum*36 19820462046252
Be. tabaci¤54 86071612078255
N. lugens¤36 74894752163258
O. fasciatus¤4863922520171
Rh. prolixus15 44145791580178
HymenopteraA. cerana cerana¤50 37311 1432150260
A. mellifera*15 31919811981253
At. cephalotes18 09354241957187
Acromyrmex echinatior17 27856511957187
C. floridanus17 06453651959187
H. saltator18 56458801964187
Li. humile16 11652571966187
Na. vitripennis*18 90019681968256
Po. barbatus17 18952481963187
So. invicta16 52251041793187
LepidopteraB. mori*14 62321382467296
Ch. suppressalis¤37 04070052029256
D. plexippus15 13051111648184
Ga. mellonella¤40 03516 3141675255
He. melpomene12 82944931580181
M. sexta¤67 9021213744189
Me. cinxia¤195 08676211598249
P. xylostella¤456 02693 9412330266
Sp. exigua¤153 61931 0711957260
Z. filipendulae¤42 73511 7091861254
PhthirapteraPe. humanus corporis*10 77320872087256
PhasmatodeaTi. cristinae17 093156 5362073260

*21 genome-sequenced and KEGG-annotated insects; 16 genome-sequenced insects; ¤15 transcriptome-sequenced insects.

Table 1.

Data statistic of iPathDB

OrderSpeciesNucleotide or protein sequencesAnnotated sequencesKO termPathway
ColeopteraDe. ponderosae¤25 36071661997261
T. castaneum*16 73320932093256
DipteraAedes aegypti*15 78620972097257
Ae. albopictus¤48 42215 7942013253
Anopheles darlingi11 43043091849187
Anopheles gambiae*12 81121302130258
Anopheles stephensi23 28749801867189
C. quinquefasciatus*18 95520682068259
D. ananassae*15 14320872087257
D. erecta*15 12320892088256
D. grimshawi*15 06620952094257
D. melanogaster*29 45321592156258
D. mojavensis*14 66520822082257
D. persimilis*16 94820372036254
D. pseudoobscura*16 92920722072256
D. sechellia*16 54420732073251
D. simulans*15 48619841983254
D. virilis*14 56720902089257
D. willistoni*15 59620792079257
D. yakuba*16 15821022101257
G. morsitans12 45747261658182
L. sericata¤146 25045 7682272263
Lu. longipalpis10 11039461483176
Ph. papatasi11 16443881506178
R. pomonella¤13 07131571032223
HemipteraA. pisum*36 19820462046252
Be. tabaci¤54 86071612078255
N. lugens¤36 74894752163258
O. fasciatus¤4863922520171
Rh. prolixus15 44145791580178
HymenopteraA. cerana cerana¤50 37311 1432150260
A. mellifera*15 31919811981253
At. cephalotes18 09354241957187
Acromyrmex echinatior17 27856511957187
C. floridanus17 06453651959187
H. saltator18 56458801964187
Li. humile16 11652571966187
Na. vitripennis*18 90019681968256
Po. barbatus17 18952481963187
So. invicta16 52251041793187
LepidopteraB. mori*14 62321382467296
Ch. suppressalis¤37 04070052029256
D. plexippus15 13051111648184
Ga. mellonella¤40 03516 3141675255
He. melpomene12 82944931580181
M. sexta¤67 9021213744189
Me. cinxia¤195 08676211598249
P. xylostella¤456 02693 9412330266
Sp. exigua¤153 61931 0711957260
Z. filipendulae¤42 73511 7091861254
PhthirapteraPe. humanus corporis*10 77320872087256
PhasmatodeaTi. cristinae17 093156 5362073260
OrderSpeciesNucleotide or protein sequencesAnnotated sequencesKO termPathway
ColeopteraDe. ponderosae¤25 36071661997261
T. castaneum*16 73320932093256
DipteraAedes aegypti*15 78620972097257
Ae. albopictus¤48 42215 7942013253
Anopheles darlingi11 43043091849187
Anopheles gambiae*12 81121302130258
Anopheles stephensi23 28749801867189
C. quinquefasciatus*18 95520682068259
D. ananassae*15 14320872087257
D. erecta*15 12320892088256
D. grimshawi*15 06620952094257
D. melanogaster*29 45321592156258
D. mojavensis*14 66520822082257
D. persimilis*16 94820372036254
D. pseudoobscura*16 92920722072256
D. sechellia*16 54420732073251
D. simulans*15 48619841983254
D. virilis*14 56720902089257
D. willistoni*15 59620792079257
D. yakuba*16 15821022101257
G. morsitans12 45747261658182
L. sericata¤146 25045 7682272263
Lu. longipalpis10 11039461483176
Ph. papatasi11 16443881506178
R. pomonella¤13 07131571032223
HemipteraA. pisum*36 19820462046252
Be. tabaci¤54 86071612078255
N. lugens¤36 74894752163258
O. fasciatus¤4863922520171
Rh. prolixus15 44145791580178
HymenopteraA. cerana cerana¤50 37311 1432150260
A. mellifera*15 31919811981253
At. cephalotes18 09354241957187
Acromyrmex echinatior17 27856511957187
C. floridanus17 06453651959187
H. saltator18 56458801964187
Li. humile16 11652571966187
Na. vitripennis*18 90019681968256
Po. barbatus17 18952481963187
So. invicta16 52251041793187
LepidopteraB. mori*14 62321382467296
Ch. suppressalis¤37 04070052029256
D. plexippus15 13051111648184
Ga. mellonella¤40 03516 3141675255
He. melpomene12 82944931580181
M. sexta¤67 9021213744189
Me. cinxia¤195 08676211598249
P. xylostella¤456 02693 9412330266
Sp. exigua¤153 61931 0711957260
Z. filipendulae¤42 73511 7091861254
PhthirapteraPe. humanus corporis*10 77320872087256
PhasmatodeaTi. cristinae17 093156 5362073260

*21 genome-sequenced and KEGG-annotated insects; 16 genome-sequenced insects; ¤15 transcriptome-sequenced insects.

iPathCons

We developed a pipeline, iPathCons, to construct insect pathways using either OGSs or transcriptomes (Figure 1).

The pipeline of iPathCons.
Figure 1.

The pipeline of iPathCons.

Data preparation

We downloaded the KEGG genes of 21 genome-sequenced insects, the OGS of 37 genome-sequenced insects and the transcriptome data of 15 species. For the 21 KEGG-annotated insects with OGS data, we compared the sequences in KEGG gene data and OGS data sets, (i) if there is a length difference, we used the long transcript; (ii) we kept those genes even they appear in only one gene data set.

KO assignment

Assigning KO terms is a crucial step in pathway construction. We used a voting system for KO assignment in the iPathCons. The protein data sets of 21 KEGG-annotated insects were divided by species. The protein sequences of each insect were formatted to build the local BLAST database, respectively, which were used as the template for KO assignment of other 16 genome-sequenced insects. For each insect need be annotated, its protein sequences were used to BLASTP against the template of every KEGG-annotated insect. The best BLASTP hit was used to assign KO terms (E-value ≤ 10−5), which has been widely used in the KOBAS (42), KAAS (43) and Blast2GO (44, 45). In this way, every protein was assigned KO terms for 21 times. The term that appears at the highest frequency (the minimum cutoff is ≥ 2) was used as the final KO assignment for the protein sequences.

A similar procedure was used to deduce pathway from the transcriptomes of 15 insects. All 37 genome-sequenced insects, including 21 KEGG-annotated and 16 iPathCons-annotated ones, were used as the template for KO assignment. The protein sequences of each genome-sequenced insect were used as the local BLAST database, respectively. The transcriptome sequences were used to BLASTP against the local BLAST database (E-value ≤ 10−5). The KO term that appeared at the highest frequency (the minimum cutoff is ≥ 2) was used as the final KO annotation.

Validation of iPathCons

We used C. quinquefasciatus gene data to validate the iPathCons. The protein sequences of C. quinquefasciatus have been annotated by the KEGG database. We removed all protein sequences of C. quinquefasciatus from the KEGG template and used them to deduce pathways using the iPathCons. The results indicated that the precision reached 95% and the coverage was 94% (E-value ≤ 10−5). We also used the transcriptome data of C. quinquefasciatus for pathway construction and obtained a similar result.

We compared the results of the iPathCons with that of other relate tool KAAS, which is a widely used pathway annotation tool provided by the KEGG database. The transcriptome of A. cerana cerana and Ga. mellonella were used to deduce pathways by both iPathCons and KAAS. The results indicated that similar number of KO terms and pathways were annotated in A. cerana cerana by two tools, whereas the iPathCons found 1675 KO terms and 255 pathways, much more than 1511 KO terms and 239 pathways annotated by the KAAS (Table 2). The iPathCons annotated significantly more contigs than the KAAS, possibly because much more templates were used in the iPathCons. However, it should be noticed that both iPathCons and KAAS relied on homology analysis to deduce the pathway. So, the results should have some false positive and need to be confirmed by molecular experiments.

Table 2.

Comparison of iPathCons and KAAS

A. cerana cerana
Ga. mellonella
iPathConsKAASiPathConsKAAS
Total nucleotide sequences50 37340 035
Annotated sequences11 143346916 3142083
KO term2150214616751511
Pathway260254255239
A. cerana cerana
Ga. mellonella
iPathConsKAASiPathConsKAAS
Total nucleotide sequences50 37340 035
Annotated sequences11 143346916 3142083
KO term2150214616751511
Pathway260254255239
Table 2.

Comparison of iPathCons and KAAS

A. cerana cerana
Ga. mellonella
iPathConsKAASiPathConsKAAS
Total nucleotide sequences50 37340 035
Annotated sequences11 143346916 3142083
KO term2150214616751511
Pathway260254255239
A. cerana cerana
Ga. mellonella
iPathConsKAASiPathConsKAAS
Total nucleotide sequences50 37340 035
Annotated sequences11 143346916 3142083
KO term2150214616751511
Pathway260254255239

Availability of iPathCons software

Both stand-alone software and the web server were provided. The stand-alone, command-line program was written using Perl language. The program consists of three parts: the main program, the ‘doc’ folder containing the index of K number and KO terms, the ‘db’ folder containing the local BLAST database. iPathCons can complete the following tasks: (i) annotating an insect transcriptome or gene sets for pathway construction; (ii) generating KGML files that can be opened by VANTED (22) and KEGG-ED (23); and (iii) generating links for each pathway showing the KEGG pathways.

Database construction by iPathDB

Database system implementation

We constructed an insect pathway database named as the iPathDB, which was developed on a Linux operating system (Redhat 5.6, Raleigh, NC, USA). The Apache HTTP server was used to handle queries from web clients through PHP scripts to perform searches. The web pages were written using html, PHP, CSS and JavaScript. The architecture of iPathDB is presented in Figure 2.

Overview of iPathDB web pages.
Figure 2.

Overview of iPathDB web pages.

Search

Users can search insect pathways using keywords for species, pathway ID and pathway name. When using species name as the search keyword, all pathways for that species will be presented. When using pathway ID or pathway name as the search keyword, the pathway will be given for all species in the database. Search results provide gene sequences, annotations and a pathway map.

Online server

An online iPathCons server was provided. The KEGG- and iPathCons-annotated gene sets from different insect orders, including Diptera, Lepidoptera, Coleoptera, Hymenoptera, Phthiraptera and Hemiptera, are used as the template for constructing insect pathways. Users can select a template according to their requirements. When the queried sequences are less than 10, the results are displayed in the Webpage directly. If the queried sequences are more than 10, a URL link of the iPathCons results will be sent to the user via e-mail.

Download page

Both FTP and HTTP download options are provided. The iPathDB FTP site is ftp://ftp.insect-genome.com/pub/iPathDB/. On the download page, insect pathway information classified by species and the stand-alone version of iPathCons can be downloaded. A phylogenetic tree is provided to show evolutionary relationships. The insect pathway files can be directly downloaded by clicking the links in the tree. On the FTP server, insect pathway information is provided by species. The zipped archive contains a readme file, sequences in FASTA format, KEGG map links, a KGML file, a pathway summary and a pathway list. In total, iPathDB contains 11 581 pathways from 52 species and 387 478 annotated sequences (Table 1).

Insect pathways

Disease-associated pathways

Insects have been studied to model human diseases. Insect disease models can provide an efficient way to study mechanisms and screen drugs. Interestingly, the results showed that 72% of human disease pathways could be found in insects (Figure 3). In total, 17 human disease-associated pathways were found in insects, including bacterial, viral and parasitic infectious disease. In contrast, only two ‘immune disease’ pathways were found, suggesting that the immune systems of insects and humans are quite different. These results suggested that insects are good candidates for modeling human infectious diseases. A successful example is that of D. melanogaster, which has been used to model cholera (46).

Human-disease and insecticide-resistance pathways. Most pathways in the ‘immune disease’ category were not found. Five pathways involved in insecticide resistance were found in the insect transcriptomes.
Figure 3.

Human-disease and insecticide-resistance pathways. Most pathways in the ‘immune disease’ category were not found. Five pathways involved in insecticide resistance were found in the insect transcriptomes.

Xenobiotic metabolism pathways

Most insects feed on plants. To protect themselves, plants produce many kinds of secondary metabolites. Insect herbivores have evolved many of xenobiotic degradation and metabolism pathways in response. Almost all insects have the pathways belonging to the category ‘xenobiotics biodegradation and metabolism’. We found that all insects contained the ‘caprolactam degradation’ pathway. Caprolactam is a pesticide intermediate (Figure 3).

Signaling pathways

Signaling pathways are important signal transduction pathways related with proteins that pass signals from outside of a cell to the inside of the cell. In total, 29 signal pathways were found. The well-studied important signal pathways exist in almost all 52 insects, including Toll-like receptor signal pathway, MAPK signal pathway, NFKB signal pathway, Notch signal pathway, etc. This suggests that these pathways are highly conserved and also play important functions in insects.

Insect hormone biosynthesis

Almost all insects undergo incomplete metamorphosis from immature nymphs, which resemble the adults, or complete metamorphosis from immature larvae, which are significantly different from the adult. Both molting and juvenile hormone control the insect metamorphosis. Hormone biosynthesis pathways were identified in all 52 insects. In the genome-sequenced insects, almost all genes in the insect hormone biosynthesis pathway were found, suggesting that this pathway is highly conserved in insects. All insects with transcriptome data had juvenile hormone epoxide hydrolase, juvenile-hormone esterase and ecdysone oxidase. Ecdysteroid 25-hydroxylase, CYP306A1 (Phm), ecdysteroid 22-hydroxylase and CYP302A1 (Dib) were found in almost all insects (Figure 4). We compared the pathway members between holometabous and hemimetabolous insects, finding no apparent difference from present data. A detail analysis of the pathway differences is worthy of further investigation. Because of the low quality of the insect transcriptome data, some genes in the insect hormone biosynthesis pathway were missing. The completeness of this pathway can be used as a parameter to estimate the quality of genome annotation or transcriptome assembly.

Insect hormone biosynthesis pathway. Most of the enzymes in the pathway could be found in all 15 insect transcriptomes; JHEH, JHE, and EO were found in all 15. Abbreviations: JHAMT, juvenile hormone acid methyltransferase; JHE, juvenile-hormone esterase; JHEH, juvenile hormone epoxide hydrolase; Nvd, cholesterol 7-dehydrogenase; Spo/Spok, CYP307A; Phm, CYP306A1; ecdysteroid 25-hydroxylase; Dib, CYP302A1, ecdysteroid 22-hydroxylase; Sad, ecdysteroid 2-hydroxylase; EO, ecdysone oxidase; SHD, ecdysone 20-monooxygenase.
Figure 4.

Insect hormone biosynthesis pathway. Most of the enzymes in the pathway could be found in all 15 insect transcriptomes; JHEH, JHE, and EO were found in all 15. Abbreviations: JHAMT, juvenile hormone acid methyltransferase; JHE, juvenile-hormone esterase; JHEH, juvenile hormone epoxide hydrolase; Nvd, cholesterol 7-dehydrogenase; Spo/Spok, CYP307A; Phm, CYP306A1; ecdysteroid 25-hydroxylase; Dib, CYP302A1, ecdysteroid 22-hydroxylase; Sad, ecdysteroid 2-hydroxylase; EO, ecdysone oxidase; SHD, ecdysone 20-monooxygenase.

Wing development pathway

Insects are characterized by having six legs and four wings, which enable diverse mobile abilities. Insect wing development is an important research topic. However, no wing development pathway is available in the KEGG or other gene network databases. Therefore, we constructed a wing development pathway after reference mining research on wing development in D. melanogaster (47–57) and Ac. pisum (58). KGML files of wing development pathways in these two species were produced. Then, those files were used as templates to construct the pathways in other species (Figure 5). To best of our knowledge, this is the first report of an insect wing development pathway. The results indicated that almost all insects have genes in this pathway. However, major parts of genes associated with wing development were missing in the flightless silkworm, B. mori. Since the silkworm has been domesticated for thousands of years, the impact of domestication on the evolution of wing development requires further investigation.

Insect wing development pathway. Arrowheads and bars indicate activation and repression, respectively. Dashed lines indicate Ubx-related regulation specific to fly halteres. Abbreviations: en, engrailed; hh, hedgehog; dpp, decapentaplegic; sal, spalt major; Ubx, Ultrabithorax; vg, vestigial; hth, homothorax; ap, apterous; Ser, Serrate; wg, wingless; Dll, Distalless; ac/sc, achaete/scute.
Figure 5.

Insect wing development pathway. Arrowheads and bars indicate activation and repression, respectively. Dashed lines indicate Ubx-related regulation specific to fly halteres. Abbreviations: en, engrailed; hh, hedgehog; dpp, decapentaplegic; sal, spalt major; Ubx, Ultrabithorax; vg, vestigial; hth, homothorax; ap, apterous; Ser, Serrate; wg, wingless; Dll, Distalless; ac/sc, achaete/scute.

Conclusion

We developed an improved analysis tool for constructing insect pathways. Both stand-alone software and web servers are provided. Users can construct insect pathways from a list of genes. An insect pathway database was also built that contains well-annotated insect pathways from 52 species.

Future study

  1. Knowledge-based construction of insect pathways relies on sequence data. Therefore, we will continually update iPathDB by adding more insect genomes once they are sequenced and published. We will also reconstruct the pathway when new versions of OGSs are released.

  2. Evolutionary analysis of insect pathway is an interesting topic that is worthy of further investigation. In the future, as more reliable insect pathways are added to iPathDB, we will carry out insect pathway conservation analysis. IPathDB will display conserved insect pathways in various insect species.

Funding

This work was supported by the National High Technology Research and Development Program (‘863’Program) of China (2012AA101505), the National Science Foundation of China (31171843, 31301691) and the Jiangsu Science Foundation for Distinguished Young Scholars (BK2012028). Funding for open access charge: BK2012028.

Conflict of interest. None declared.

References

1

Drosophila 12 Genomes Consortium
(
2007
)
Evolution of genes and genomes on the Drosophila phylogeny
.
Nature
,
450
,
203
218
.

2

Adams
M.D.
Celniker
S.E.
Holt
R.A.
et al. . (
2000
)
The genome sequence of Drosophila melanogaster
.
Science
,
287
,
2185
2195
.

3

Smith
C.D.
Zimin
A.
Holt
C.
et al. . (
2011
)
Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile)
.
Proc. Natl Acad. Sci. USA
,
108
,
5673
5678
.

4

Smith
C.R.
Smith
C.D.
Robertson
H.M.
et al. . (
2011
)
Draft genome of the red harvester ant Pogonomyrmex barbatus
.
Proc. Natl Acad. Sci. USA
,
108
,
5667
5672
.

5

Wurm
Y.
Wang
J.
Riba-Grognuz
O.
et al. . (
2011
)
The genome of the fire ant Solenopsis invicta
.
Proc. Natl Acad. Sci. USA
,
108
,
5679
5684
.

6

Nygaard
S.
Zhang
G.
Schiott
M.
et al. . (
2011
)
The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming
.
Genome Res.
,
21
,
1339
1348
.

7

Suen
G.
Teiling
C.
Li
L.
et al. . (
2011
)
The genome sequence of the leaf-cutter ant Atta cephalotes reveals insights into its obligate symbiotic lifestyle
.
PLoS Genet.
,
7
,
e1002007
.

8

Bonasio
R.
Zhang
G.
Ye
C.
et al. . (
2010
)
Genomic comparison of the ants Camponotus floridanus and Harpegnathos saltator
.
Science
,
329
,
1068
1071
.

9

Werren
J.H.
Richards
S.
Desjardins
C.A.
et al. . (
2010
)
Functional and evolutionary insights from the genomes of three parasitoid Nasonia species
.
Science
,
327
,
343
348
.

10

Nene
V.
Wortman
J.R.
Lawson
D.
et al. . (
2007
)
Genome sequence of Aedes aegypti, a major arbovirus vector
.
Science
,
316
,
1718
1723
.

11

Holt
R.A.
Subramanian
G.M.
Halpern
A.
et al. . (
2002
)
The genome sequence of the malaria mosquito Anopheles gambiae
.
Science
,
298
,
129
149
.

12

Arensburger
P.
Megy
K.
Waterhouse
R.M.
et al. . (
2010
)
Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics
.
Science
,
330
,
86
88
.

13

Heliconius Genome Consortium
(
2012
)
Butterfly genome reveals promiscuous exchange of mimicry adaptations among species
.
Nature
,
487
,
94
98
.

14

Zhan
S.
Merlin
C.
Boore
J.L.
et al. . (
2011
)
The monarch butterfly genome yields insights into long-distance migration
.
Cell
,
147
,
1171
1185
.

15

Kirkness
E.F.
Haas
B.J.
Sun
W.
et al. . (
2010
)
Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle
.
Proc. Natl Acad. Sci. USA
,
107
,
12168
12173
.

16

Megy
K.
Emrich
S.J.
Lawson
D.
et al. . (
2012
)
VectorBase: improvements to a bioinformatics resource for invertebrate vector genomics
.
Nucleic Acids Res.
,
40
,
D729
D734
.

17

Honeybee Genome Sequencing Consortium
(
2006
)
Insights into social insects from the genome of the honeybee Apis mellifera
.
Nature
,
443
,
931
949
.

18

Xia
Q.
Zhou
Z.
Lu
C.
et al. . (
2004
)
A draft sequence for the genome of the domesticated silkworm (Bombyx mori)
.
Science
,
306
,
1937
1940
.

19

Soria-Carrasco
V.
Gompert
Z.
Comeault
A.A.
et al. . (
2014
)
Stick insect genomes reveal natural selection's role in parallel speciation
.
Science
,
344
,
738
742
.

20

Tribolium Genome Sequencing Consortium
,
Richards
S.
Gibbs
R.A.
et al. 
. (
2008
)
The genome of the model beetle and pest Tribolium castaneum
.
Nature
,
452
,
949
955
.

21

International Aphid Genomics Consortium
(
2010
)
Genome sequence of the pea aphid Acyrthosiphon pisum
.
PLoS Biol.
,
8
,
e1000313
.

22

You
M.
Yue
Z.
He
W.
et al. . (
2013
)
A heterozygous moth genome provides insights into herbivory and detoxification
.
Nat. Genet.
,
45
,
220
225
.

23

Wang
X.
Fang
X.
Yang
P.
et al. . (
2014
)
The locust genome provides insight into swarm formation and long-distance flight
.
Nat. Commun.
,
5
,
2957
.

24

Mi
H.
Dong
Q.
Muruganujan
A.
et al. . (
2010
)
PANTHER version 7: improved phylogenetic trees, orthologs and collaboration with the gene ontology consortium
.
Nucleic Acids Res.
,
38
,
D204
D210
.

25

Ashburner
M.
Ball
C.A.
Blake
J.A.
et al. . (
2000
)
Gene ontology: tool for the unification of biology. The gene ontology consortium
.
Nat. Genet.
,
25
,
25
29
.

26

Kanehisa
M.
Goto
S.
Furumichi
M.
et al. . (
2010
)
KEGG for representation and analysis of molecular networks involving diseases and drugs
.
Nucleic Acids Res.
,
38
,
D355
D360
.

27

Croft
D.
O'Kelly
G.
Wu
G.
et al. . (
2011
)
Reactome: a database of reactions, pathways and biological processes
.
Nucleic Acids Res.
,
39
,
D691
D697
.

28

Eichelbaum
M.
Altman
R.B.
Ratain
M.
et al. . (
2009
)
New feature: pathways and important genes from PharmGKB
.
Pharmacogenet. Genomics
,
19
,
403
.

29

Marygold
S.J.
Leyland
P.C.
Seal
R.L.
et al. . (
2013
)
FlyBase: improvements to the bibliography
.
Nucleic Acids Res.
,
41
,
D751
D757
.

30

Legeai
F.
Shigenobu
S.
Gauthier
J.P.
et al. . (
2010
)
AphidBase: a centralized bioinformatic resource for annotation of the pea aphid genome
.
Insect Mol. Biol.
19
5
12
.

31

Papanicolaou
A.
Gebauer-Jung
S.
Blaxter
M.L.
et al. . (
2008
)
ButterflyBase: a platform for lepidopteran genomics
.
Nucleic Acids Res.
,
36
,
D582
D587
.

32

Kim
H.S.
Murphy
T.
Xia
J.
et al. . (
2010
)
BeetleBase in 2010: revisions to provide comprehensive genomic information for Tribolium castaneum
.
Nucleic Acids Res.
,
38
,
D437
D442
.

33

Tang
W.
Yu
L.
He
W.
et al. . (
2014
)
DBM-DB: the diamondback moth genome database
.
Database
,
2014
,
bat087
.

34

Zhan
S.
Reppert
S.M.
(
2013
)
MonarchBase: the monarch butterfly genome database
.
Nucleic Acids Res.
,
41
,
D758
D763
.

35

Wurm
Y.
Uva
P.
Ricci
F.
et al. . (
2009
)
Fourmidable: a database for ant genomics
.
BMC Genomics
,
10
,
5
.

36

Lawson
D.
Arensburger
P.
Atkinson
P.
et al. . (
2007
)
VectorBase: a home for invertebrate vectors of human pathogens
.
Nucleic Acids Res.
,
35
,
D503
D505
.

37

Duan
J.
Li
R.
Cheng
D.
et al. . (
2010
)
SilkDB v2.0: a platform for silkworm (Bombyx mori) genome biology
.
Nucleic Acids Res.
,
38
,
D453
D456
.

38

Grabherr
M.G.
Haas
B.J.
Yassour
M.
et al. . (
2011
)
Full-length transcriptome assembly from RNA-Seq data without a reference genome
.
Nat. Biotechnol.
,
29
,
644
652
.

39

Huang
X.
Madan
A.
(
1999
)
CAP3: a DNA sequence assembly program
.
Genome. Res.
,
9
,
868
877
.

40

Kanehisa
M.
Goto
S.
Sato
Y.
et al. . (
2014
)
Data, information, knowledge and principle: back to metabolism in KEGG
.
Nucleic Acids Res.
,
42
,
D199
D205
.

41

Pruitt
K.D.
Tatusova
T.
Maglott
D.R.
(
2007
)
NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins
.
Nucleic Acids Res.
,
35
,
D61
D65
.

42

Xie
C.
Mao
X.
Huang
J.
et al. . (
2011
)
KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases
.
Nucleic Acids Res.
,
39
,
W316
W322
.

43

Moriya
Y.
Itoh
M.
Okuda
S.
et al. . (
2007
)
KAAS: an automatic genome annotation and pathway reconstruction server
.
Nucleic Acids Res.
,
35
,
W182
W185
.

44

Conesa
A.
Gotz
S.
(
2008
)
Blast2GO: a comprehensive suite for functional analysis in plant genomics
.
Int. J. Plant. Genomics
,
2008
,
619832
.

45

Altschul
S.F.
Gish
W.
Miller
W.
et al. . (
1990
)
Basic local alignment search tool
.
J. Mol. Biol.
,
215
,
403
410
.

46

Blow
N.S.
Salomon
R.N.
Garrity
K.
et al. . (
2005
)
Vibrio cholerae infection of Drosophila melanogaster mimics the human disease cholera
.
PLoS Pathog.
,
1
,
e8
.

47

Held
L.I.
Jr.
(
2002
)
Bristles induce bracts via the EGFR pathway on Drosophila legs
.
Mech. Dev.
,
117
,
225
234
.

48

Azpiazu
N.
Morata
G.
(
2000
)
Function and regulation of homothorax in the wing imaginal disc of Drosophila
.
Development
,
127
,
2685
2693
.

49

Panganiban
G.
Rubenstein
J.L.
(
2002
)
Developmental functions of the Distal-less/Dlx homeobox genes
.
Development
,
129
,
4371
4386
.

50

Lawrence
P.A.
Casal
J.
Struhl
G.
(
1999
)
Hedgehog and engrailed: pattern formation and polarity in the Drosophila abdomen
.
Development
,
126
,
2431
2439
.

51

Garcia-Bellido
A.
Santamaria
P.
(
1972
)
Developmental analysis of the wing disc in the mutant engrailed of Drosophila melanogaster
.
Genetics
,
72
,
87
104
.

52

Nfonsam
L.E.
Cano
C.
Mudge
J.
et al. . (
2012
)
Analysis of the transcriptomes downstream of Eyeless and the Hedgehog, Decapentaplegic and Notch signaling pathways in Drosophila melanogaster
.
PLoS One
,
7
,
e44583
.

53

Wilson
T.G.
(
1981
)
Expression of phenotypes in a temperature-sensitive allele of the apterous mutation in Drosophila melanogaster
.
Dev. Biol.
,
85
,
425
433
.

54

Seto
E.S.
Bellen
H.J.
(
2006
)
Internalization is required for proper wingless signaling in Drosophila melanogaster
.
J. Cell Biol.
,
173
,
95
106
.

55

McKay
D.J.
Estella
C.
Mann
R.S.
(
2009
)
The origins of the Drosophila leg revealed by the cis-regulatory architecture of the Distalless gene
.
Development
,
136
,
61
71
.

56

Slattery
M.
Ma
L.
Negre
N.
et al. . (
2011
)
Genome-wide tissue-specific occupancy of the Hox protein Ultrabithorax and Hox cofactor Homothorax in Drosophila
.
PLoS One
,
6
,
e14686
.

57

Glicksman
M.A.
Brower
D.L.
(
1988
)
Expression of the sex combs reduced protein in Drosophila larvae
.
Dev. Biol.
,
127
,
113
118
.

58

Brisson
J.A.
Ishikawa
A.
Miura
T.
(
2010
)
Wing development genes of the pea aphid and differential gene expression between winged and unwinged morphs
.
Insect Mol. Biol.
,
19
63
73
.

Author notes

Citation details: Zhang,Z., Yin, C., Liu,Y., et al. iPathCons and iPathDB: an improved insect pathway construction tool and the database. Database (2014) Vol. 2014: article ID bau105; doi:10.1093/database/bau105

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data