BACKGROUND
cotton bollworm  and oriental tobacco budworm  are two sibling noctuid species of lepidoptera.
they are distributed in almost same region from 50°s to 50°n and from 45°s to 45°n, respectively.
h. armigera are slightly broader than oriental tobacco budworms  <cit> .
in the field, similar external morphology makes them easily confused.
interestingly, these two species can mate under artificial manipulation and produce offspring.
however, when female h. armigera mated with male h. assulta, the first filial generations are all males  <cit> .
this result further confirms that they are two distinct species  <cit> .
under natural conditions, because of differences in sex pheromone composition, the two species seldom mate.
their sex pheromones comprise cis-11-hexadecenal  and cis-9-hexadecenal , but the compositions are reversed in the two species.
the ratio of these two components is 97: <dig> in h. armigera, but it is 7: <dig> in h. assulta <cit> .
in addition, the host ranges of these two species are significantly different.
h. armigera is a euryphagous insect whose host range includes  <dig> families of over  <dig> different plants.
however, h. assulta is an oligophagous insect, they are mainly feeding on the plants of the solanaceae, for example, tobacco and hot pepper  <cit> .
although each species have their own preferred host plants, both of them love feed on tobacco and hot pepper  <cit> .
these similarity and difference may be depend on the host plant selection by adult, or depend on the food digestion or detoxification enzymes from larvae.
host plant selection is a complicated and continuous process.
the color, odor and shape of the plants will affect the insects choice on host plants, among which odors are a critically important element for the lifestyle and reproduction of an insect species  <cit> .
accordingly, insects with different feeding habits possess their own specific odor identification and odor-binding proteins  <cit> .
at the same time, enzymes for food digestion and detoxification are also important factors for insect growth and development.
these enzymes probably effect on the survival of insects, consequently affecting on the host range of an insect species.
therefore, identifying enzymes related to insect development and feeding habits will benefit to research on insect host range and on insect pest control.
the two sibling species are non-model insects.
their genome sequences are not available till now.
to identify differentially expressed genes from h. armigera and h. assulta, digital gene expression tag  profile libraries were constructed and sequenced using high throughput second-generation sequencing technology  <cit> .
a dge-tag profile is according to the theory and method of sage  combined with high-throughput sequencing technology  <cit> .
in this study, eight dge-tag libraries were constructed and sequenced for four developmental stages  of the two species.
differentially expressed transcripts or genes between h. armigera and h. assulta were analyzed by bioinformatics.
most growth and development related genes have similar expression modes.
the differentially expressed genes are mainly focus on olfactory-related genes and enzymes for food detoxification or digestion.
therefore, the two sibling species represent a good model for host plant selection and adaptability.
these results also provide valuable data for insect pest control.
RESULTS
the main identifying characteristics of h. armigera and h. assulta
the two insect species in genus helicoverpa, cotton bollworm  and oriental tobacco budworm , are important insect pests of crops in china.
they have similar external morphology.
figure  <dig> shows some taxonomic characteristics to distinguish these two species.
their eggs, larvae and pupae look like nearly same .
only under the microscope, according to some taxonomic characteristics they can be distinguished .
the adult is the easiest stage to be distinguished by entomologist, because there are some special speckles and stripes on the wings .
to show the relationship of h. armigera and h. assulta, multiple sequence alignments of spanning the 18s rrna across  <dig> species from  <dig> orders  and the expansion segment of the coi gene across  <dig> species of lepidopteran moths  were constructed and supplied as additional files.
the results provide clues about the evolutionary origin of the phytophagous noctuidae.
the sister group of h. assulta and h. armigera is clustered on a clade.
some hypotheses on the sister group relatedness based on morphology are concordant with our molecular results.
sequencing of dge-tag libraries and unique tag annotation
dge-tag profile libraries were constructed from total rna of h. armigera and h. assulta for four development stages .
the summary sequence results are shown in table  <dig>  low frequency tags were discounted under the assumption that many could have arisen through sequencing errors such as base substitution, deletion or addition at a single position  <cit> .
therefore, after eliminating low quality tags , copy numbers less than two and adaptor sequences, the remaining reads were called clean tags, of which more than 50% were singletons , which is typically observed in sage experiments  <cit> .
we obtained approximately  <dig> million 17nt clean tags from both insect species.
their total unique clean tag  numbers were also similar at approximately  <dig>  .
unique tag-to-gene assignments were conducted for the four development stages of h. armigera and h. assulta using soapdenovo program just permitting 1 bp mismatch  <cit> .
on average, more than 75% of the uni-tags of h. armigera were mapped on transcripts; however, only  <dig> % uni-tags of h. assulta mapped on transcripts.
the total numbers of transcripts or genes were  <dig>  and  <dig>  for h. armigera and h. assulta, respectively, among which the annotated transcripts or genes were  <dig>  and  <dig> , respectively .
the illumina short-reads sequence of h. armigera and h. assulta were submitted to ncbi sequence read archive under the accession number of srr <dig> and srr <dig>  respectively.
although we obtained similar amounts of total unique clean tags from h. armigera and h. assulta, the numbers of uni-tags obtained from the different developmental stages in each insect were quite different .
in theory, each uni-tag should be derived from one transcript  <cit> .
according to this theory, the embryo stage has the highest number of transcripts compared with any other stages in the two species .
unexpectedly, the larval stage of h. armigera has the lowest number of transcripts,  <dig>  .
the gene annotation result indicated that  <dig> % and  <dig> % uni-tags of h. armigera and h. assulta corresponded to est sequences in the h. armigera transcriptome library.
the number of identified genes in h. armigera is higher than in h. assulta; however, the number of identified genes has no significantly difference at each developmental stage of the same species .
global analysis of differentially expressed genes between the two species
the unique clean tags provide transcripts information for one species.
using a venn diagram, developmental stage-specific transcripts and coexpressed transcripts between two to four developmental stages were shown between of h. armigera and h. assulta .
the analysis results revealed that the minimum numbers of coexpressed transcripts existed between the larval and adult stage in both h. armigera and h. assulta .
this indicated that the biggest differences exist between these two stages among the four developmental stages.
however, in h. armigera, the embryo and adult stage are probably “the closest neighbors” and have the most amount of  <dig> coexpressed transcripts or genes.
in h. assulta, the largest number of coexpressed transcripts was  <dig> between the embryo and larvae stage.
the uni-tag annotation results were also analyzed for differential expression and coexpressed transcripts or genes using a venn diagram .
the copy number of each unique tag provides quantitative information for the abundance of the transcripts or genes detected by the tags.
using the tag copy number, we can roughly estimate the expression level of each transcripts or gene.
the dynamics of gene expression can be reflected by up- or down-regulation among the four development stages by pairwise comparisons .
overall, the changes in gene expression levels between two developmental stages in h. assulta are more extreme than in h. armigera.
coexpressed transcripts or genes between two species
comparing the unique clean tags between h. armigera and h. assulta, approximately 30% transcripts or genes are coexpressed in each developmental stage .
this reflects the real situation of coexpression transcripts in these two insects.
because there are no whole genome annotation information for the two insect species, most of the differentially or coexpressed genes are unknown proteins, hypothetical proteins, enzymes or cytoskeletal proteins, which are annotated according to the h. armigera transcriptome results.
in terms of annotated genes, about 67% to 85% genes in h. armigera and h. assulta are coexpression during the four developmental stages .
go analysis also confirmed that only in the larval stage there were more functional genes in h. assulta than in h. armigera .
development and host range related transcripts or genes between two insect species
to explain why the two species have these similarities and differences, we focused on comparing the transcripts or genes that are related to growth and development, food digestion or detoxification enzymes, and host plant recognition.
we identified  <dig> and  <dig> growth and development related transcripts,  <dig> and  <dig> transcripts for food digestion or detoxification related enzymes, and  <dig> and  <dig> olfaction related transcripts from h. armigera and h. assulta, respectively.
the relative expression levels  for these transcripts in each developmental stage of the two species and the annotation results are listed in additional files  <dig>   <dig> and  <dig>  in summary, the amounts of each type of transcript are similar between the two species .
the biggest difference in the number is odorant binding proteins .
there are  <dig> obp transcripts in h. armigera, but only  <dig> in h. assulta .
expression patterns of possible host-range-related transcripts or genes
the expression patterns of transcripts can be divided into two categories by tag copy number among the four developmental stages between two species .
the first type is the similar expression pattern.
for example, in table  <dig>  the trypsin- <dig> gene just express at the embryo stage in the two species .
the cytochrome c oxidase polypeptide iii gene has nearly same high expression levels at all four development stages .
more information is shown in additional files  <dig>   <dig>  and  <dig>  genes and transcripts with similar expression patterns were not further analyzed in this study.
the other expression pattern is showing a significant difference between two species.
for example, the obp <dig> gene is a carrier of odor molecules, which can protect odor molecules from enzymatic degradation  <cit> .
thus, obp <dig> is likely to be an important gene in host plant selection.
in this study, the expression pattern of obp <dig> was significantly different among the four development stages between the two insect species .
these kinds of transcripts or genes are probably the main reasons underlying the differences between the two species.
the expressions of these types of transcripts or genes were confirmed by reverse transcription polymerase chain reaction  .
most of expression patterns are nearly consistent with the digital expression tag copy number.
these transcripts or genes should be further studied in host selection and adaptation.
* s = similarity; d = difference.
CONCLUSIONS
in this paper, we systematically analyzed the differences and similarities between the two sibling insect species, h. armigera and h. assulta.
these characteristics make the two sibling species fitting for research on host range, evolution and pest control.
by comparing the tag copy number of each developmental stage, we identified some differentially expressed transcripts or genes that are probably associated with host plant recognition and food digestion or detoxification.
these genes provide important clues for further study.
sage is a method of large-scale gene expression analysis  <cit> .
it is an ‘open’ system that permits the relative expression levels of almost all transcripts in an organism.
dge-tag profiles are a development of this technology using second-generation high throughput sequencing technology  <cit> .
the dge-tag profile results indicated that even between two sibling species, only 30% of transcripts are coexpressed in each developmental stage.
the annotation results indicated that 67–85% genes are coexpressed in each developmental stage between the two species .
these results further confirm that they are distinct two species at the genome level.
the unique clean tags number can provide quantitative information for the number of transcripts detected by tags.
using traditional sage technology,  <dig>  to  <dig>  tags could be collected, which represent  <dig>  to  <dig>  unique tags  <cit> .
in this study, the high-throughput approach was adopted to implement the tag sequencing protocol on the illumina platform  <cit> .
using this technology, we obtained more than  <dig> million clean tags from h. armigera and h. assulta, respectively.
these data far more exceed the saturation requirements of sequencing  <cit> .
the total unique clean tags were  <dig>  and  <dig>  for each species, which probably represent the total transcripts in the whole life cycle of the two insects.
a total of  <dig>  to  <dig>  transcripts were identified for each developmental stage .
these results indicated that the data sets were suitable for analyzing the differential expression of transcripts or genes.
in this study, we found that embryo stage expressed the most transcripts or genes in both species.
unexpectedly, the larval stage expressed the lowest numbers of transcripts in h. armigera .
because the total clean tag number far more exceeded the saturation requirements, the difference is unlikely to be caused by sequencing bias.
the most significant difference between the two insects is in their host ranges, which should be reflected at the larval stage.
actually, the selection and adaptation of the host plant is decided by the insect’s internal factors and external environmental stimuli.
when insects choose a plant as a host, they will be spawning and feeding, and then growth and reproduction on the plant.
during this process, as external environmental factors, plant volatile odors are the most important cues for host plant selection.
at the same time, as internal factors, insect obps not only can selectively bind certain types of odor molecules, but also can remove toxic substances and protect the odor molecules from enzymatic degradation  <cit> .
then, the downstream ors  will be activated, the chemical odor molecule information will be converted into an electrical signal, spread to the central nervous system, triggering an insect behavioral response  <cit> .
in this research, we identified  <dig> and  <dig> obp-related transcripts or genes from h. armigera and h. assulta, respectively .
this is the biggest difference among all the transcripts types.
this is probably the main reason why these two sibling species have different host ranges.
we also identified two and four or-related transcripts or genes from h. armigera and h. assulta, respectively.
these transcripts should be further studied to increase our understanding of the host range difference between the two species.
in addition to host selection, the other important aspect is host adaptation.
during feeding, insect will inevitably swallow some poisonous secondary metabolites from plants.
therefore, insects have to develop an adaptation mechanism involving a series of detoxification enzymes  <cit> .
these detoxification enzymes include cytochrome p450-dependent monooxygenases , glutathione-s-transferases and carboxylesterases   <cit> .
in this study, we found that there are more transcripts or genes for p450s, coes and atpases in h. assulta than in h. armigera .
go analysis also confirmed that only in the larval stage there are more functional genes in h. assulta than in h. armigera .
therefore, we suspected that because oriental tobacco budworm has a narrower host range, more detoxification enzymes would help them increase food detoxification and utilization efficiency.
this should be further investigated in a future study.
the two species are also important agricultural insect pests.
considering all the similarities and differences between the two sibling species, we think they are a good model insect-pair for developing species-selective rnai technology.
many studies have shown that rnai is feasible technology in insect pest control  <cit> .
the appeal of rnai technology in pest control is that it is possible to design the pesticide to target only a single species or a group of related species, with minimal threat to other organisms  <cit> .
to this end, it is necessary to identify species-specific target genes.
the present study represents an effective strategy for identifying differentially expressed genes from related species that do not have genome sequences.
using dge-tag profile technology, we identified many differentially expression transcripts or genes from the two sibling insect species.
these genes not only provide clues for host range difference studies, but may also represent important targets for species-specific control by rnai technology.

 cotton bollworm and oriental tobacco budworm are two sibling noctuid species of lepidoptera. they are distributed in almost same region from 50°s to 50°n and from 45°s to 45°n, respectively. h. armigera are slightly broader than oriental tobacco budworms. in the field, similar external morphology makes them easily confused. interestingly, these two species can mate under artificial manipulation and produce offspring. however, when female h. armigera mated with male h. assulta, the first filial generations are all males. this result further confirms that they are two distinct species. under natural conditions, because of differences in sex pheromone composition, the two species seldom mate. their sex pheromones comprise cis-11-hexadecenal and cis-9-hexadecenal, but the compositions are reversed in the two species. the ratio of these two components is 97: in h. armigera, but it is 7: in h. assulta. in addition, the host ranges of these two species are significantly different. h. armigera is a euryphagous insect whose host range includes families of over different plants. however, h. assulta is an oligophagous insect, they are mainly feeding on the plants of the solanaceae, for example, tobacco and hot pepper. although each species have their own preferred host plants, both of them love feed on tobacco and hot pepper. these similarity and difference may be depend on the host plant selection by adult, or depend on the food digestion or detoxification enzymes from larvae. host plant selection is a complicated and continuous process. the color, odor and shape of the plants will affect the insects choice on host plants, among which odors are a critically important element for the lifestyle and reproduction of an insect species. accordingly, insects with different feeding habits possess their own specific odor identification and odor-binding proteins. at the same time, enzymes for food digestion and detoxification are also important factors for insect growth and development. these enzymes probably effect on the survival of insects, consequently affecting on the host range of an insect species. therefore, identifying enzymes related to insect development and feeding habits will benefit to research on insect host range and on insect pest control. the two sibling species are non-model insects. their genome sequences are not available till now. to identify differentially expressed genes from h. armigera and h. assulta, digital gene expression tag profile libraries were constructed and sequenced using high throughput second-generation sequencing technology. a dge-tag profile is according to the theory and method of sage combined with high-throughput sequencing technology. in this study, eight dge-tag libraries were constructed and sequenced for four developmental stages of the two species. differentially expressed transcripts or genes between h. armigera and h. assulta were analyzed by bioinformatics. most growth and development related genes have similar expression modes. the differentially expressed genes are mainly focus on olfactory-related genes and enzymes for food detoxification or digestion. therefore, the two sibling species represent a good model for host plant selection and adaptability. these results also provide valuable data for insect pest control.  the main identifying characteristics of h. armigera and h. assulta the two insect species in genus helicoverpa, cotton bollworm and oriental tobacco budworm, are important insect pests of crops in china. they have similar external morphology. figure  shows some taxonomic characteristics to distinguish these two species. their eggs, larvae and pupae look like nearly same. only under the microscope, according to some taxonomic characteristics they can be distinguished. the adult is the easiest stage to be distinguished by entomologist, because there are some special speckles and stripes on the wings. to show the relationship of h. armigera and h. assulta, multiple sequence alignments of spanning the 18s rrna across species from orders and the expansion segment of the coi gene across species of lepidopteran moths were constructed and supplied as additional files. the results provide clues about the evolutionary origin of the phytophagous noctuidae. the sister group of h. assulta and h. armigera is clustered on a clade. some hypotheses on the sister group relatedness based on morphology are concordant with our molecular results. sequencing of dge-tag libraries and unique tag annotation dge-tag profile libraries were constructed from total rna of h. armigera and h. assulta for four development stages. the summary sequence results are shown in table  low frequency tags were discounted under the assumption that many could have arisen through sequencing errors such as base substitution, deletion or addition at a single position. therefore, after eliminating low quality tags, copy numbers less than two and adaptor sequences, the remaining reads were called clean tags, of which more than 50% were singletons, which is typically observed in sage experiments. we obtained approximately million 17nt clean tags from both insect species. their total unique clean tag numbers were also similar at approximately. unique tag-to-gene assignments were conducted for the four development stages of h. armigera and h. assulta using soapdenovo program just permitting 1 bp mismatch. on average, more than 75% of the uni-tags of h. armigera were mapped on transcripts; however, only % uni-tags of h. assulta mapped on transcripts. the total numbers of transcripts or genes were and for h. armigera and h. assulta, respectively, among which the annotated transcripts or genes were and, respectively. the illumina short-reads sequence of h. armigera and h. assulta were submitted to ncbi sequence read archive under the accession number of srr and srr respectively. although we obtained similar amounts of total unique clean tags from h. armigera and h. assulta, the numbers of uni-tags obtained from the different developmental stages in each insect were quite different. in theory, each uni-tag should be derived from one transcript. according to this theory, the embryo stage has the highest number of transcripts compared with any other stages in the two species. unexpectedly, the larval stage of h. armigera has the lowest number of transcripts,. the gene annotation result indicated that % and % uni-tags of h. armigera and h. assulta corresponded to est sequences in the h. armigera transcriptome library. the number of identified genes in h. armigera is higher than in h. assulta; however, the number of identified genes has no significantly difference at each developmental stage of the same species. global analysis of differentially expressed genes between the two species the unique clean tags provide transcripts information for one species. using a venn diagram, developmental stage-specific transcripts and coexpressed transcripts between two to four developmental stages were shown between of h. armigera and h. assulta. the analysis results revealed that the minimum numbers of coexpressed transcripts existed between the larval and adult stage in both h. armigera and h. assulta. this indicated that the biggest differences exist between these two stages among the four developmental stages. however, in h. armigera, the embryo and adult stage are probably “the closest neighbors” and have the most amount of coexpressed transcripts or genes. in h. assulta, the largest number of coexpressed transcripts was between the embryo and larvae stage. the uni-tag annotation results were also analyzed for differential expression and coexpressed transcripts or genes using a venn diagram. the copy number of each unique tag provides quantitative information for the abundance of the transcripts or genes detected by the tags. using the tag copy number, we can roughly estimate the expression level of each transcripts or gene. the dynamics of gene expression can be reflected by up- or down-regulation among the four development stages by pairwise comparisons. overall, the changes in gene expression levels between two developmental stages in h. assulta are more extreme than in h. armigera. coexpressed transcripts or genes between two species comparing the unique clean tags between h. armigera and h. assulta, approximately 30% transcripts or genes are coexpressed in each developmental stage. this reflects the real situation of coexpression transcripts in these two insects. because there are no whole genome annotation information for the two insect species, most of the differentially or coexpressed genes are unknown proteins, hypothetical proteins, enzymes or cytoskeletal proteins, which are annotated according to the h. armigera transcriptome results. in terms of annotated genes, about 67% to 85% genes in h. armigera and h. assulta are coexpression during the four developmental stages. go analysis also confirmed that only in the larval stage there were more functional genes in h. assulta than in h. armigera. development and host range related transcripts or genes between two insect species to explain why the two species have these similarities and differences, we focused on comparing the transcripts or genes that are related to growth and development, food digestion or detoxification enzymes, and host plant recognition. we identified and growth and development related transcripts, and transcripts for food digestion or detoxification related enzymes, and and olfaction related transcripts from h. armigera and h. assulta, respectively. the relative expression levels for these transcripts in each developmental stage of the two species and the annotation results are listed in additional files and in summary, the amounts of each type of transcript are similar between the two species. the biggest difference in the number is odorant binding proteins. there are obp transcripts in h. armigera, but only in h. assulta. expression patterns of possible host-range-related transcripts or genes the expression patterns of transcripts can be divided into two categories by tag copy number among the four developmental stages between two species. the first type is the similar expression pattern. for example, in table  the trypsin- gene just express at the embryo stage in the two species. the cytochrome c oxidase polypeptide iii gene has nearly same high expression levels at all four development stages. more information is shown in additional files and genes and transcripts with similar expression patterns were not further analyzed in this study. the other expression pattern is showing a significant difference between two species. for example, the obp gene is a carrier of odor molecules, which can protect odor molecules from enzymatic degradation. thus, obp is likely to be an important gene in host plant selection. in this study, the expression pattern of obp was significantly different among the four development stages between the two insect species. these kinds of transcripts or genes are probably the main reasons underlying the differences between the two species. the expressions of these types of transcripts or genes were confirmed by reverse transcription polymerase chain reaction. most of expression patterns are nearly consistent with the digital expression tag copy number. these transcripts or genes should be further studied in host selection and adaptation. * s = similarity; d = difference.  in this paper, we systematically analyzed the differences and similarities between the two sibling insect species, h. armigera and h. assulta. these characteristics make the two sibling species fitting for research on host range, evolution and pest control. by comparing the tag copy number of each developmental stage, we identified some differentially expressed transcripts or genes that are probably associated with host plant recognition and food digestion or detoxification. these genes provide important clues for further study. sage is a method of large-scale gene expression analysis. it is an ‘open’ system that permits the relative expression levels of almost all transcripts in an organism. dge-tag profiles are a development of this technology using second-generation high throughput sequencing technology. the dge-tag profile results indicated that even between two sibling species, only 30% of transcripts are coexpressed in each developmental stage. the annotation results indicated that 67–85% genes are coexpressed in each developmental stage between the two species. these results further confirm that they are distinct two species at the genome level. the unique clean tags number can provide quantitative information for the number of transcripts detected by tags. using traditional sage technology, to tags could be collected, which represent to unique tags. in this study, the high-throughput approach was adopted to implement the tag sequencing protocol on the illumina platform. using this technology, we obtained more than million clean tags from h. armigera and h. assulta, respectively. these data far more exceed the saturation requirements of sequencing. the total unique clean tags were and for each species, which probably represent the total transcripts in the whole life cycle of the two insects. a total of to transcripts were identified for each developmental stage. these results indicated that the data sets were suitable for analyzing the differential expression of transcripts or genes. in this study, we found that embryo stage expressed the most transcripts or genes in both species. unexpectedly, the larval stage expressed the lowest numbers of transcripts in h. armigera. because the total clean tag number far more exceeded the saturation requirements, the difference is unlikely to be caused by sequencing bias. the most significant difference between the two insects is in their host ranges, which should be reflected at the larval stage. actually, the selection and adaptation of the host plant is decided by the insect’s internal factors and external environmental stimuli. when insects choose a plant as a host, they will be spawning and feeding, and then growth and reproduction on the plant. during this process, as external environmental factors, plant volatile odors are the most important cues for host plant selection. at the same time, as internal factors, insect obps not only can selectively bind certain types of odor molecules, but also can remove toxic substances and protect the odor molecules from enzymatic degradation. then, the downstream ors will be activated, the chemical odor molecule information will be converted into an electrical signal, spread to the central nervous system, triggering an insect behavioral response. in this research, we identified and obp-related transcripts or genes from h. armigera and h. assulta, respectively. this is the biggest difference among all the transcripts types. this is probably the main reason why these two sibling species have different host ranges. we also identified two and four or-related transcripts or genes from h. armigera and h. assulta, respectively. these transcripts should be further studied to increase our understanding of the host range difference between the two species. in addition to host selection, the other important aspect is host adaptation. during feeding, insect will inevitably swallow some poisonous secondary metabolites from plants. therefore, insects have to develop an adaptation mechanism involving a series of detoxification enzymes. these detoxification enzymes include cytochrome p450-dependent monooxygenases, glutathione-s-transferases and carboxylesterases. in this study, we found that there are more transcripts or genes for p450s, coes and atpases in h. assulta than in h. armigera. go analysis also confirmed that only in the larval stage there are more functional genes in h. assulta than in h. armigera. therefore, we suspected that because oriental tobacco budworm has a narrower host range, more detoxification enzymes would help them increase food detoxification and utilization efficiency. this should be further investigated in a future study. the two species are also important agricultural insect pests. considering all the similarities and differences between the two sibling species, we think they are a good model insect-pair for developing species-selective rnai technology. many studies have shown that rnai is feasible technology in insect pest control. the appeal of rnai technology in pest control is that it is possible to design the pesticide to target only a single species or a group of related species, with minimal threat to other organisms. to this end, it is necessary to identify species-specific target genes. the present study represents an effective strategy for identifying differentially expressed genes from related species that do not have genome sequences. using dge-tag profile technology, we identified many differentially expression transcripts or genes from the two sibling insect species. these genes not only provide clues for host range difference studies, but may also represent important targets for species-specific control by rnai technology. 
high-throughput sequencing results indicated that we have got more than million 17nt clean tags from both species, respectively. the number of unique clean tags was nearly same in both species. according to the gene annotation results, we identified and olfaction related transcripts from h. armigera and h. assulta, respectively. among the olfaction related transcripts, more odorant binding protein and g protein-coupled receptor were identified in h. armigera than in h. assulta. p carboxypeptidase and atpase in h. assulta than in h. armigera. these differences partially explain that because of the narrow host plant range of h. assulta, more detoxification enzymes would help them increase the food detoxification and utilization efficiency. cotton bollworm and oriental tobacco budworm are noctuid sibling species. as serious agricultural insect pests, cotton bollworms are euryphagous insects, but oriental tobacco budworms are oligophagous insects. these genes will be useful information for studying on the evolution of host plant selection. under artificial manipulation, they can mate and produce fertile offspring. to identify the differentially expressed genes that affect host recognition and host adaptation between the two species, we constructed digital gene expression tag profiles for four developmental stages of the two species. it also provides some important target genes for insect species-specific control by rnai technology. this study supplied some differentially expressed genes affecting host selection and adaptation between the two sibling species.
 cotton bollworm and oriental tobacco budworm are noctuid sibling species. under artificial manipulation, they can mate and produce fertile offspring. as serious agricultural insect pests, cotton bollworms are euryphagous insects, but oriental tobacco budworms are oligophagous insects. to identify the differentially expressed genes that affect host recognition and host adaptation between the two species, we constructed digital gene expression tag profiles for four developmental stages of the two species. high-throughput sequencing results indicated that we have got more than million 17nt clean tags from both species, respectively. the number of unique clean tags was nearly same in both species.  according to the gene annotation results, we identified and olfaction related transcripts from h. armigera and h. assulta, respectively. at the same time, and transcripts of digestion enzymes were identified from the two species. among the olfaction related transcripts, more odorant binding protein and g protein-coupled receptor were identified in h. armigera than in h. assulta. among the digestion enzymes, there are more detoxification enzyme, e.g. p carboxypeptidase and atpase in h. assulta than in h. armigera. these differences partially explain that because of the narrow host plant range of h. assulta, more detoxification enzymes would help them increase the food detoxification and utilization efficiency.  this study supplied some differentially expressed genes affecting host selection and adaptation between the two sibling species. these genes will be useful information for studying on the evolution of host plant selection. it also provides some important target genes for insect species-specific control by rnai technology. developmenthost plant rangetranscriptsdigital gene expression tag profile sibling speciesdifferential expression genehelicoverpa armigerahelicoverpa assulta 