BACKGROUND
insects comprise the largest group of animals on earth and most species are heavily dependent on chemoreception for survival and reproduction.
due to ecological adaptations, the gene families involved in chemoreception have evolved, both in terms of gene numbers and sequences, to become highly divergent across different insect taxa
<cit> .
furthermore, in a particular species, expression of the genes is regulated differently depending on, for instance, sex, tissue, life stage, or physiological state.
the sequencing of insect genomes and more recently also of transcriptomes from chemosensory appendages, have paved the way for identification and characterization of the gene families that are important for chemoreception e.g.
<cit> .
several multi-gene families encode proteins with crucial roles in olfaction and taste, including both receptors and non-receptor proteins.
receptors encoded by three large gene families are involved in chemoreception and expressed in sensory neurons, located mainly on the antennae and to a lesser extent on other sensory appendages
<cit> .
these families encode the odorant receptors , gustatory receptors , and ionotropic receptors
<cit> .
the ors and grs belong to the same receptor superfamily
<cit> .
while the ors detect volatile chemicals , the grs are responsible for contact chemoreception and also for detection of carbon dioxide
<cit> .
the irs constitute a more ancient family of chemoreceptors
<cit> , including both conserved "antennal" irs involved in olfaction, as well as "divergent" irs that seem to play a role in taste
<cit> .
the sensory neuron membrane proteins  are expressed in pheromone-responding neurons in drosophila and lepidoptera, and have in some cases been shown to be important for proper pheromone responses
<cit> .
finally, the odorant binding proteins  are small soluble proteins, highly abundant in the sensillum lymph
<cit> .
obps bind odor molecules, most of which are hydrophobic, and transport them through the hydrophilic environment in the sensillum to the membrane-bound receptor.
obps have also been shown to improve the specificity of pheromone-detecting neurons
<cit> , but see
<cit> .
chemosensory gene families of insects have been most extensively studied in two model organisms that both belong to the order diptera , namely drosophila melanogaster, which feeds on yeast and sugar both as larva and adult, and anopheles gambiae, which feeds on organic material as an aquatic larvae and on blood as an adult female or nectar as a male e.g.
<cit> .
the herbivores are a third major group within the diptera in which the chemosensory genes so far had not been studied.
here we performed the first analysis of the chemosensory gene families in a herbivorous dipteran species, namely the hessian fly, mayetiola destructor say .
the large gall midge family cecidomyiidae  is an example of a group of flies that uses plants as food and contains a large number of important pests of agriculture and forestry
<cit> .
cecidomyiids belong to the suborder nematocera  and are phylogenetically interesting because they are positioned between mosquitoes, also in the nematocera, and drosophila, in the other suborder of the diptera, the brachycera
<cit> .
plant-feeding cecidomyiid larvae have highly specialized antagonistic relationships with host plants and have been called ‘reprogrammers of plant genomes’
<cit>  because of their ability to manipulate the plant to create a specialized nutritive tissue that provides a superior diet for the larva
<cit> , a trait that is thought to have contributed to adaptive radiation
<cit> .
most gall midges have a very narrow host range
<cit> .
other characteristics of cecidomyiids are a very short adult lifespan of 1– <dig> days or less, during which all reproductive activities must be completed, and mate location by the male is mediated by species-specific female-produced sex pheromones
<cit> .
cecidomyiids are also special in that adults either do not feed or only rarely feed on water or perhaps nectar
<cit> .
the hessian fly is a serious pest on wheat  in many of the world’s major wheat-growing regions
<cit> .
it is one of the most thoroughly studied plant-feeding insect species, with particular attention on behavioral and sensory ecology, genetics contributing to plant interactions, pheromone communication, pest management, and plant resistance
<cit> .
the hessian fly larva has an interaction with plants that is common for plant pathogens but rare for insects, being the target of a ‘gene-for-gene’ plant defense, mediated by resistance genes to which it can adapt via mutations in a matching avirulence gene
<cit> .
like other gall midges, adults of m. destructor have a maximum lifespan of 1– <dig> days, which limits their behavioral repertoire
<cit> .
essentially, males emerge from pupation sites in or near the soil and take flight to find virgin females, whereas females emerge with a full complement of mature eggs, release a pheromone to attract a mate, and then continuously search for and lay eggs on a large number of plants until death occurs
<cit> .
chemical cues are crucial for mate finding, host finding, and oviposition behavior and are, thus, clearly essential for the reproductive success of the hessian fly.
by having an interesting phylogenetic position between drosophila and mosquitoes, as well as very different behavioral and ecological traits, the hessian fly is an attractive model for studies of the chemosensory gene families that underlie pheromone detection and host plant finding and assessment.
the hessian fly genome sequence has recently become available .
within the framework of the hessian fly genomics consortium, we annotated a total of  <dig> genes involved in chemoreception , corresponding to the first identification of the chemosensory gene families in the cecidomyiidae.
in the present study, we used illumina rnaseq to sequence transcriptomes of the male and female antennae and terminal abdominal segments, the latter being involved in mating and in the female also in sex pheromone production and host selection behavior.
the generated reads were then mapped to the genome sequence to examine the expression profiles of the chemosensory genes in these four tissues.
we also analyzed global expression profiles of predicted transcripts and performed gene ontology  annotation.
in particular, our results reveal that most of the chemosensory genes have up-regulated expression in the antennae, while a few have up-regulated expression in the terminal abdominal tissues.
in addition, a strikingly large number of genes, especially or genes, have sex-specific antennal expression, which is likely to reflect the sex-divergent behaviors of adult hessian flies.
RESULTS
mapping and transcript prediction
we performed illumina sequencing and subsequently mapped the sequenced reads to the current hessian fly genome assembly to predict transcripts and analyze gene expression profiles.
the sequencing yielded ca.
314 m read-pairs, of which ca.
82 m were derived from the male antennal library, 75 m from female antennae, 77 m from female terminal abdomens, and 79 m from male terminal abdomens.
the illumina reads have been submitted to the sequence read archive  at ncbi .
the proportion of reads that successfully mapped to the hessian fly genome assembly was high in all four transcriptomes, ranging from 75% to 80%, despite strict mapping parameters that allowed for only 2 bp mismatches to account for the high sequence similarity among recently duplicated chemosensory genes.
the mapping results were used by cufflinks to predict  <dig>  <dig> transcripts expressed in female antennae ,  <dig>  <dig> transcripts in male antennae ,  <dig>  <dig> in female terminal abdomen , and  <dig>  <dig> in male terminal abdomen .
based on these transcripts,  <dig>  <dig> genes and  <dig>  <dig> isoforms were predicted.
the gc content was 34-35% in all transcriptomes.
global transcriptome profiling
gene ontology annotation
go annotation
<cit> , using blast2go
<cit> , was performed to compare the four transcriptomes with regards to the frequency of transcript-associated go terms.
the proportion of transcript sequences with significant blast hits was similar in the four transcriptomes, i.e.
57%  for female antennae, 59%  for male antennae, 65%  for female terminal abdomen, and 60%  for male terminal abdomen.
the top-hit species distribution  showed that most hits were represented by sequences from the yellow-fever mosquito , followed by the southern house mosquito  and the malaria mosquito .
the vinegar fly  came fourth despite the fact that ca.
five times as many protein sequences from d. melanogaster have been deposited in the nr database as compared to sequences from a. aegypti .
the top blast hits for all predicted transcripts as well as their fpkm  values and genome locations are given in additional file
<dig>
for the female antennae, 41% of the total transcripts had assigned go terms in addition to blast hits.
corresponding numbers in the other transcriptomes were 43% for male antennae, 50% for female terminal abdomen, and 44% for male terminal abdomen.
the four transcriptomes turned out to be very similar, both with respect to "molecular function"  and "biological process"  go annotation.
the only notable differences were found in the molecular function annotation where the "odorant binding" term only was present in the two antennal tissues, and the "signal receptor activity" term was almost twice as abundant in the two antennal tissues as compared to the two terminal abdominal tissues .
this result is likely to reflect the -sensory role of the main sensory organ, the antenna.
however, it also indicates that genes coding for proteins with functions  common to all tissue types dominate in this type of analysis that provides no information on the expression levels of the individual genes themselves.
pair-wise transcriptome comparisons
the expression levels of the  <dig>  <dig> genes predicted by cufflinks were used for global transcriptome profiling.
transcriptomes were compared pair-wise, i.e.
male antennae vs. female antennae; male terminal abdomen vs. female terminal abdomen; male antennae vs. male terminal abdomen; and female antennae vs. female terminal abdomen.
the male and female antennal transcriptomes had very similar global expression profiles, whereas male and female terminal abdomens were more different, both from each other as well as from the two antennal transcriptomes .
when divergence based on fpkm values was analyzed using jensen-shannon  distance, male and female antennae showed a low level of differentiation .
in contrast, the distances in all other pair-wise comparisons were larger, ranging from ca.
<dig>  to  <dig>  .
the global expression profile similarity between male and female antennae, and their distinctiveness from the other transcriptomes were confirmed in a principle component analysis , and in pair-wise regression analyses.
the latter analyses indicated that the female and male terminal abdominal segments had an overall higher expression than the female and male antennae, respectively .
the male antennae had a similar global expression profile as the female antennae, as shown by a close to 1: <dig> relationship between the expression levels .
in general, higher expression was found in the female terminal abdomen as compared to the male terminal abdomen .
expression of the chemosensory genes
overall expression profiles
based on our manual annotation of the hessian fly genome, a total of  <dig> chemosensory genes from five gene families  were predicted.
predicted amino acid sequences of these genes are available in additional file
<dig>  and phylogenies will be included in the forthcoming publication of the hessian fly genome .
in contrast to the global comparison where the male and female antennal transcriptomes demonstrated the most similar expression profiles, the female and male terminal abdomens had the most similar profiles when only the  <dig> chemosensory-related genes were included .
regression analysis demonstrated a close to 1: <dig> relationship between the expression levels of the chemosensory genes in the two terminal abdominal tissues , although the overall expression level was low in these tissues.
most of the chemosensory genes had the highest expression in the female or male antennae as compared to the female or male terminal abdomens , respectively, although a few genes had higher expression in the terminal abdominal tissues .
in contrast to the global analysis, the expression profiles of the male and female antennal tissues were clearly different when only the  <dig> genes related to chemosensation were included in the analysis .
fpkm values for the chemosensory genes in the four tissues are presented in additional file
<dig>
odorant receptors
based on our annotation, the current assembly of the hessian fly genome contains  <dig> or genes, including a gene for the conserved co-receptor orco and six pseudogenes.
our transcriptome analysis showed that  <dig> of the or genes were expressed in at least one of the analyzed tissues.
thus, almost 20%  of the or genes were not expressed .
a substantially larger number of or genes were represented by reads in the female antennae  as compared to in the male antennae  .
thus, while a large proportion  of the or genes in the genome was expressed in the antennae of the female, only ca.
half  of them were expressed in the antennae of the male.
female and male terminal abdomens expressed  <dig> and  <dig> or genes, respectively, including orco .
thus, relatively large numbers of or genes were expressed outside of the main olfactory organ in both sexes, although the expression level was low in general .
the orco gene had the highest level of expression in both sexes, followed by or <dig> for which the fpkm value in the male antennae was 90% of the fpkm value for orco.
gene family
a
b
anumbers within brackets refer to the number of annotated genes of each chemosensory gene family in the hessian fly genome.
bexpression defined as genes having at least  <dig> mapped read-pairs.
in the pair-wise comparison between the female antennae and the female terminal abdomen,  <dig> or genes had higher expression in the antenna , whereas two or genes, or <dig> and or <dig>  had higher expression in the terminal abdomen .
in the corresponding comparison for the male ,  <dig> or genes were expressed at a higher level in the antennae, whereas four or genes  had higher expression in the terminal abdomen.
when comparing expression in the antennae of males vs. females,  <dig> or genes showed higher expression in females, whereas  <dig> were more highly expressed in males.
in addition, several of the or subfamilies in our dendrogram contained receptors of which the majority showed either female or male antennal up-regulation .
for instance, the genes for most of the receptors within the clades formed by or86- <dig>  or98- <dig>  and or102- <dig> were more highly expressed in the female antennae as compared to the male antennae, whereas most of the receptor genes within the clade formed by or111- <dig> were more highly expressed in the male antennae.
furthermore,  <dig> of the or genes showed no expression in the antennae of either males or females, while  <dig> or genes were similarly expressed in the antennae of the two sexes.
finally, five or genes showed higher level of expression in the female terminal abdomen as compared to the male terminal abdomen, whereas seven or genes had higher expression in the male terminal abdomen in relation to the female terminal abdomen .
gustatory receptors
we annotated  <dig> members of the gr gene family in the genome, including genes for the conserved carbon dioxide receptors
<cit> , three putative sugar receptor genes , and three pseudogenes.
however, only  <dig> of the gr genes were expressed in the tissues analyzed here .
these included the genes for the three putative carbon dioxide receptors and the three sugar receptors.
only slightly larger numbers of gr genes were expressed in the antennae as compared to the terminal abdominal tissues .
three of the gr genes  had similar levels of expression in all four tissues .
in the females, six gr genes had up-regulated expression in the antennae as compared to the terminal abdomen.
these included the genes for the carbon dioxide receptors  that had the highest expression of all the gr genes, and the putative sugar receptor gene gr <dig>  in contrast, another putative sugar receptor gene, gr <dig>  had up-regulated expression in the female terminal abdomen as compared to the female antennae .
similar to the females, gr1- <dig> and gr <dig> had higher expression in the male antennae as compared to the male terminal abdomen, whereas gr <dig> had higher expression in the terminal abdomen .
in contrast to the female, the expression of gr <dig> was also up-regulated in the male antennae as compared to terminal abdomen, whereas gr <dig>  was expressed at a similar level in the male antennae and terminal abdomen.
gr <dig> was also the only gr gene that was differentially expressed between the two antennal transcriptomes.
finally, four gr genes, including the putative sugar receptor genes gr4- <dig> were more highly expressed in the male terminal abdomen than in the female terminal abdomen.
none of the gr genes were more highly expressed in the female terminal abdomen than in the male terminal abdomen.
ionotropic receptors
in total,  <dig> ir genes, including four pseudogenes were annotated in the hessian fly genome.
of these,  <dig> genes encoded "antennal" irs and  <dig> encoded "divergent" irs.
we found hessian fly orthologues for most of the conserved antennal ir genes identified in other diptera
<cit> , including ir8a, ir21a, ir25a, ir41a, ir60a, ir64a, ir68a, ir76b, ir93a, and seven members of the ir <dig> group.
however, ir31a, ir40a, and ir92a, which are found in other diptera, were not present in the hessian fly genome and could also not be identified in our transcriptome data.
the remaining  <dig> mdir genes were regarded as members of the divergent class of insect irs.
twelve of the identified ir genes  were not expressed, whereas all the antennal ir genes showed expression in at least one of the analyzed tissues .
larger numbers of ir genes were expressed in the antennal tissues as compared to the terminal abdominal tissues .
the gene for the broadly occurring co-receptor ir25a had the highest expression of all the mdir genes.
expression of this gene was high in the antennae of both sexes and in the male terminal abdomen.
expression of  <dig> ir genes was higher in the female antennae as compared to the female terminal abdomen.
of these genes,  <dig> encoded antennal irs and four encoded divergent irs .
although ir60a is classified as an "antennal" ir, its gene showed a non-significant tendency for higher expression in the female terminal abdomen than in the female antennae.
apart from that, only one of the divergent ir genes  was significantly up-regulated in the female terminal abdomen.
in the males,  <dig> ir genes  had higher expression in the antennae as compared to the terminal abdomen, while three ir genes had higher expression in the terminal abdomen.
similar to the females, ir116a was one of them, but also ir87b and ir87c.
only one of the ir genes, ir75c, had up-regulated expression in the female antennae compared to those of the male, representing the only example of differential ir gene expression between the antennae of the two sexes.
comparing the two terminal abdominal tissues, five ir genes, including both antennal and divergent receptor genes, had higher expression in the male terminal abdomen, whereas expression of two ir genes was higher in the female terminal abdomen .
sensory neuron membrane proteins
the hessian fly genome assembly contained  <dig> orthologues of snmp <dig> and one snmp <dig>  all of them were expressed in the antennae, and six of them  were expressed in both terminal abdominal tissues .
in general, the snmp genes had high expression, and the most highly expressed gene  had a higher fpkm value than orco.
in the females, snmp1a-c, e and f had up-regulated expression in the antennae as compared to the terminal abdomen, while one  was up-regulated in the terminal abdomen, both in relation to the female antennae and to the male terminal abdomen .
furthermore, snmp <dig> had slightly higher expression in the female terminal abdomen than in the female antennae, but the difference was not significant.
the highest expression of snmp <dig> was found in the male terminal abdomen.
in the males, snmp1b-f had up-regulated expression in the antennae as compared to the terminal abdomen.
comparing expression of snmp genes between the antennae of the two sexes, snmp1a was greatly up-regulated  in the female.
in contrast, snmp1b and 1f had male-biased antennal expression, while snmp1b, 1c, and 1e were similarly expressed in the antennae of the two sexes .
odorant binding proteins
we annotated  <dig> obp genes in the hessian fly genome assembly, including two pseudogenes.
generally, obps are classified into different phylogenetic groups depending on the number of conserved cysteine residues.
classic obps have six cysteine residues, whereas members of the plus-c class have additional cysteines and one characteristic proline
<cit> .
five of the hessian fly obps  belonged to the plus-c class.
all but one  of the  <dig> mdobp genes showed expression in at least one of the analyzed tissues .
essentially all obp genes were expressed in the antennae of both sexes, while slightly lower numbers of genes were expressed in the terminal abdominal tissues .
overall, the most highly expressed obp genes had much higher expression levels  than the most highly expressed or or snmp genes .
the plus-c class of obp genes  had comparatively low expression levels in all four tissues .
in the female, most  of the obp genes were more highly expressed in the antennae as compared to the terminal abdomen.
only one obp gene  showed a higher level of expression in the female terminal abdomen than in the female antennae , but the overall expression level of this gene was comparatively low.
the males were similar to the females by having  <dig> obp genes with higher expression in the antennae as compared to the terminal abdomen.
these included most of the genes that also were antennally up-regulated in the female.
three obp genes  had higher expression in the male terminal abdomen as compared to the male antennae.
however, despite the fact that most obp genes had their highest expression in the two antennal tissues, their expression level in the terminal abdomen tissues was in general much higher than the expression of the other chemosensory gene families in the abdominal tissues .
a few obp genes showed differential expression between the male and female antennae .
only two obp genes  showed higher levels of expression in the female antennae, whereas six genes  had higher expression in the male antennae.
nine of the obp genes showed differential expression levels between the two terminal abdominal tissues; five of them  were more highly expressed in the female, while four  were more highly expressed in the male.
discussion
chemosensory genes in the hessian fly in relation to drosophila and mosquitoes
the chemosensory gene families have been extensively studied in various drosophila and mosquito species
<cit> , but had not been analyzed in plant-feeding dipterans, such as gall midges.
the hessian fly in the cecidomyiidae family is both phylogenetically
<cit>  and ecologically
<cit>  well separated from flies and mosquitoes.
the specific habits and highly evolved relationship to its host also makes studies of the hessian fly chemosensory genes interesting for a more complete understanding of olfaction in diptera.
the current hessian fly genome assembly contains genes for  <dig> ors,  <dig> grs,  <dig> irs,  <dig> snmps, and  <dig> obps.
in comparison to drosophila and mosquitoes, the number of hessian fly or genes is relatively large, whereas the number of genes encoding taste receptors  is small
<cit> , suggesting that olfaction is the dominant chemical sense in this species.
this might be related to the short adult lifespan, with olfaction being a more time-efficient modality for host and mate assessment, occurring during flight at a distance rather than after landing via contact chemoreception.
while clear orthologues of ors are rarely found in distantly related species, the antennal irs appear conserved across insect species
<cit> .
however, the hessian fly genome lacked genes for three antennal irs  that are present in other diptera, suggesting that also some olfactory functions have been lost in this herbivorous species.
insects generally have two snmp genes , but some species have multiple paralogues of snmp <dig> <cit> .
snmp <dig> has been shown to be important for pheromone detection in drosophila <cit> .
the hessian fly genome contains six snmp <dig> genes, which is the largest number identified in any insect so far.
it is possible that this large number reflects the use of a multi-component sex pheromone in the hessian fly
<cit> , and that the detection of the individual pheromone components requires specific or-snmp <dig> combinations.
sex-specific antennal expression of chemosensory genes
transcriptome profiling using rnaseq data has been shown to accurately predict gene expression levels
<cit> , and has recently been used to analyze the expression of chemosensory genes in several other insect species, including d. melanogaster and the a. gambiae mosquito
<cit> .
here, we used a transcriptome profiling approach to analyze sex- and tissue specific expression in the hessian fly.
our global transcriptome analysis that included all the ca.
<dig>  <dig> predicted genes indicated that the gene expression profiles of the male and female antennae were highly similar to each other, while expression profiles of the antennae were different from the terminal abdominal tissues.
in addition, the go annotation also indicated very similar frequencies of transcript-associated go terms in the two antennal transcriptomes.
in fact, the results from the latter analysis were also similar to previous results from go annotation of antennal transcriptomes of two bark beetle species and the moth manduca sexta <cit> , suggesting that "transcript frequencies" in terms of presence/absence are conserved across insect orders.
when only the  <dig> chemosensory genes were included in the analysis, the expression profile differed widely between the four tissues.
the majority of the chemosensory genes had up-regulated expression levels in the antennae as compared to the terminal abdominal tissues.
furthermore and in contrast to the global analysis, a clear difference in gene expression between the antennae of males and females was found.
although all analyzed chemosensory gene families contained some genes with expression levels that differed between the sexes, the most striking difference was found among the or genes, for which ca.
62% of the genes with any expression showed at least a three-fold significant difference between male and female antennae, and ca.
50% of them showed at least a ten-fold difference, with the majority being female-biased.
thus, the proportion of or genes with sex-biased expression was substantially larger than previously found in d. melanogaster <cit> , a species in which both sexes feed  on the same food, but similar to a. gambiae <cit> , where males and females feed on entirely different food.
similar to the mosquito, the sex-biased expression in the hessian fly is likely to reflect the different ecological interests of males and females
<cit> .
in both sexes, the gene for the co-receptor orco
<cit>  had the highest expression levels of all or genes.
this is consistent with its broad expression and role in forming heteromeric complexes with the conventional ors, required for the odor response of the neuron and localization of the or in the cell membrane
<cit> .
during their ca.
18 h adult existence, male hessian flies do just one thing, namely search for and mate with virgin females, carrying enough sperm to fertilize up to  <dig> females
<cit> .
finding females is accomplished by long-range attraction to a female-produced sex pheromone, comprised of five to seven compounds
<cit> .
similar to moths e.g.
<cit> , it is likely that the subfamily formed by the receptors with male-biased expression in the hessian fly  includes the receptors for the pheromone compounds.
we currently test this hypothesis using heterologous assays.
interestingly, two of the mdor genes  in the clade with putative pheromone receptors had high expression levels also in the females, and other receptors in this clade had moderate female expression.
thus, if these ors are the pheromone receptors, it is likely that females also are able to detect some of their pheromone compounds.
in contrast to the males, for which there is no evidence in cecidomyiids of attraction to host plants
<cit> , adult females depend on plant chemicals to find and select host plants for their offspring
<cit> .
females, during their ca.
9–24 h existence, crawl from an eclosion site in the soil to a perch where they extend the abdominal segments to release a sex pheromone and, after mating, fly to find host plants for their eggs.
attraction to host plants is mediated by an array of volatile chemicals
<cit>  and discrimination between host and non-host plants would presumably require a larger and different set of ors as compared to detection of the pheromone.
thus, it is likely that the  <dig> or genes with female-biased expression are involved in host plant recognition.
similar to the hessian fly, most of the ors with sex-biased expression in a. gambiae were up-regulated in the female, i.e.
the sex that blood-feeds and finds sites for oviposition
<cit> .
interestingly, the or gene expression profiles in the hessian fly antennae were often correlated to the phylogenetic position of the ors in the tree, i.e.
the majority of the receptors within a phylogenetic subfamily showed male- or, in most cases, female-biased antennal expression.
it is possible that ors with similar sequences detect structurally similar chemicals.
thus, this result might also reflect the relative importance of host odors for females versus the pheromone compounds for males.
a correlation between sex-specific expression bias and position in the phylogeny was also found for the ors in a recent study on two ant species
<cit> .
this is likely to represent adaptations that enable individuals to carry out the tasks that are specific to their sex or caste.
the majority of the conserved antennal mdir genes had relatively high expression in the antennae, but only one of them  showed a clear difference in expression level between the sexes.
this suggests that in contrast to the ors, receptors in this more ancient family
<cit>  detect compounds that might be of general importance for insects, regardless of sex.
of the antennal ir genes, mdir60a was not expressed in the antenna.
similarly, benton et al.
<cit>  did not detect expression of ir60a in antennal samples of d. melanogaster.
to our knowledge, a ligand for this receptor has not been identified yet in any species.
the carbon dioxide receptor genes mdgr1- <dig> had high expression in the antennae , suggesting that this organ is responsible for carbon dioxide detection, which is similar to drosophila
<cit> , but different from mosquitoes that detect this gas with the maxillary palps
<cit> .
only one gr gene, mdgr <dig>  had male antennal-biased expression.
in a phylogenetic analysis , this receptor  grouped together with dmgr43a, which has been shown to monitor fructose levels in the fly brain
<cit> .
the reason why a tentative fructose receptor would have high expression in the antenna of male hessian fly is not clear, but this receptor could possibly be important if males do feed on nectar
<cit> .
since snmp <dig> is important for pheromone detection in d. melanogaster <cit> , we expected to find the highest expression levels of the mdsnmp <dig> genes in the male antennae of the hessian fly.
however, only two  of the six snmp <dig> orthologues present in the hessian fly genome had higher expression in the male antennae as compared to the female antennae.
these two snmps were also relatively highly expressed in the female antennae.
interestingly, snmp1a had >200-fold higher expression in the female antennae as compared to the male antennae.
expression of snmp genes in females might, again, suggest that females also are able to detect components of the sex pheromone, or alternatively, that the role of snmp <dig> might not be restricted to detection of sex pheromones.
expression of snmp <dig> in the antennae of both sexes has also been found in other dipteran and non-dipteran species
<cit> .
in general, the mdobp genes had the highest expression levels of all analyzed chemosensory gene families, and the majority of them were most highly expressed in the antennae.
only a few obp genes were differentially expressed between the antennae of males and females, which is similar to a. gambiae and the yellow fever mosquito, aedes aegypti <cit> .
similar to the other gene families, the mdobp genes that did show sex-biased expression might be important for the detection of compounds that elicit sex-specific behaviors.
expression in terminal abdominal tissues
all chemosensory gene families in the hessian fly contained several genes that were represented by reads also in the terminal abdominal samples.
the obp gene family was the most broadly expressed, followed by the snmps, while the or gene family had the most antennal-specific expression.
nevertheless, all gene families contained genes that were significantly more highly expressed in the male or female terminal abdominal segments as compared to the antennae.
in the male, these included four ors, two grs, three irs, one snmp, and three obps.
in the female, two or genes and one gene from each of the other families were more highly expressed in the terminal abdomen than in the antennae.
these observations suggest that chemosensation is not restricted to the antennae of gall midges.
female hessian flies examine the leaf surface with the tip of the abdomen prior to oviposition
<cit> , and oviposition is induced by chemical cues from host plants
<cit> , consistent with the expression of chemoreceptors in this tissue.
similarly, the gustatory sense of d. melanogaster in not restricted to head appendages and taste sensilla are present on the legs, wings, and female genitalia
<cit> .
in addition, or59b in d. melanogaster is predominantly expressed in the male accessory gland , and a recent study showed that several or genes are expressed in the testes of a. gambiae and that orco is involved in chemical-induced sperm activation
<cit> .
these results are consistent with the expression of several mdors including mdorco in the male terminal abdominal tissue analyzed here.
similar to the present study, other studies have reported that snmp <dig> and snmp <dig> as well as the obps can be expressed in variety of tissues in addition to the chemosensory organs, suggesting that these proteins can have physiological roles independent of olfaction
<cit> .
the snmps and obps that were expressed in the female hessian fly terminal abdominal segments where the pheromone gland is found
<cit> , might be involved in the binding and transport of pheromone molecules before these are volatilized and released.
the gr and divergent ir families contained the largest proportions of genes for which we did not detect expression in the present study.
similarly, the five genes  of the plus-c class of obps showed only low levels of expression.
it is possible that these genes are expressed mainly in legs or wings, or in yet other tissues
<cit> , which we did not analyze.
indeed, several ir genes and plus-c obp genes showed high expression in the labella and tarsi of a. aegypti <cit> .
alternatively, these and the other "non-expressed" chemosensory genes might be expressed only in the larvae
<cit> , which are likely to use contact chemical cues in order to locate a suitable feeding site inside the wheat plant
<cit> .
CONCLUSIONS
by conducting a comprehensive manual gene annotation and transcriptome analysis we found clear sex- and tissue-specific expression patterns of the chemosensory genes in the hessian fly.
while many studies of these genes have been performed on diptera, especially drosophila and mosquitoes, this is the first of a fly that feeds on plants.
the expression profile difference between male and female hessian fly was more pronounced than that observed in d. melanogaster, but similar to a. gambiae, a species that also displays sex-divergent chemically-mediated behaviors.
we found the most extreme sex-specific expression difference among the or genes, which is likely to reflect the sex-divergent olfactory behaviors of the short-lived adults.
in addition, the relatively large number of or genes in the genome in combination with the few gr genes and divergent ir genes, suggest that long-range attraction to hosts and mates is more important than contact assessment in this species, in which the adults have reduced mouthparts and feeding behavior has not been documented.
we also found that several genes of each gene family had expression that was not restricted to the antennae, indicating that they serve multiple and tissue-specific physiological roles.
future studies should focus on characterizing the function of the proteins encoded by these differentially expressed chemosensory genes in order to consolidate the link between the genetics and ecology of hessian fly olfaction.

 insects comprise the largest group of animals on earth and most species are heavily dependent on chemoreception for survival and reproduction. due to ecological adaptations, the gene families involved in chemoreception have evolved, both in terms of gene numbers and sequences, to become highly divergent across different insect taxa. furthermore, in a particular species, expression of the genes is regulated differently depending on, for instance, sex, tissue, life stage, or physiological state. the sequencing of insect genomes and more recently also of transcriptomes from chemosensory appendages, have paved the way for identification and characterization of the gene families that are important for chemoreception e.g.. several multi-gene families encode proteins with crucial roles in olfaction and taste, including both receptors and non-receptor proteins. receptors encoded by three large gene families are involved in chemoreception and expressed in sensory neurons, located mainly on the antennae and to a lesser extent on other sensory appendages. these families encode the odorant receptors, gustatory receptors, and ionotropic receptors. the ors and grs belong to the same receptor superfamily. while the ors detect volatile chemicals, the grs are responsible for contact chemoreception and also for detection of carbon dioxide. the irs constitute a more ancient family of chemoreceptors, including both conserved "antennal" irs involved in olfaction, as well as "divergent" irs that seem to play a role in taste. the sensory neuron membrane proteins are expressed in pheromone-responding neurons in drosophila and lepidoptera, and have in some cases been shown to be important for proper pheromone responses. finally, the odorant binding proteins are small soluble proteins, highly abundant in the sensillum lymph. obps bind odor molecules, most of which are hydrophobic, and transport them through the hydrophilic environment in the sensillum to the membrane-bound receptor. obps have also been shown to improve the specificity of pheromone-detecting neurons, but see. chemosensory gene families of insects have been most extensively studied in two model organisms that both belong to the order diptera, namely drosophila melanogaster, which feeds on yeast and sugar both as larva and adult, and anopheles gambiae, which feeds on organic material as an aquatic larvae and on blood as an adult female or nectar as a male e.g.. the herbivores are a third major group within the diptera in which the chemosensory genes so far had not been studied. here we performed the first analysis of the chemosensory gene families in a herbivorous dipteran species, namely the hessian fly, mayetiola destructor say. the large gall midge family cecidomyiidae is an example of a group of flies that uses plants as food and contains a large number of important pests of agriculture and forestry. cecidomyiids belong to the suborder nematocera and are phylogenetically interesting because they are positioned between mosquitoes, also in the nematocera, and drosophila, in the other suborder of the diptera, the brachycera. plant-feeding cecidomyiid larvae have highly specialized antagonistic relationships with host plants and have been called ‘reprogrammers of plant genomes’ because of their ability to manipulate the plant to create a specialized nutritive tissue that provides a superior diet for the larva, a trait that is thought to have contributed to adaptive radiation. most gall midges have a very narrow host range. other characteristics of cecidomyiids are a very short adult lifespan of 1– days or less, during which all reproductive activities must be completed, and mate location by the male is mediated by species-specific female-produced sex pheromones. cecidomyiids are also special in that adults either do not feed or only rarely feed on water or perhaps nectar. the hessian fly is a serious pest on wheat in many of the world’s major wheat-growing regions. it is one of the most thoroughly studied plant-feeding insect species, with particular attention on behavioral and sensory ecology, genetics contributing to plant interactions, pheromone communication, pest management, and plant resistance. the hessian fly larva has an interaction with plants that is common for plant pathogens but rare for insects, being the target of a ‘gene-for-gene’ plant defense, mediated by resistance genes to which it can adapt via mutations in a matching avirulence gene. like other gall midges, adults of m. destructor have a maximum lifespan of 1– days, which limits their behavioral repertoire. essentially, males emerge from pupation sites in or near the soil and take flight to find virgin females, whereas females emerge with a full complement of mature eggs, release a pheromone to attract a mate, and then continuously search for and lay eggs on a large number of plants until death occurs. chemical cues are crucial for mate finding, host finding, and oviposition behavior and are, thus, clearly essential for the reproductive success of the hessian fly. by having an interesting phylogenetic position between drosophila and mosquitoes, as well as very different behavioral and ecological traits, the hessian fly is an attractive model for studies of the chemosensory gene families that underlie pheromone detection and host plant finding and assessment. the hessian fly genome sequence has recently become available. within the framework of the hessian fly genomics consortium, we annotated a total of genes involved in chemoreception, corresponding to the first identification of the chemosensory gene families in the cecidomyiidae. in the present study, we used illumina rnaseq to sequence transcriptomes of the male and female antennae and terminal abdominal segments, the latter being involved in mating and in the female also in sex pheromone production and host selection behavior. the generated reads were then mapped to the genome sequence to examine the expression profiles of the chemosensory genes in these four tissues. we also analyzed global expression profiles of predicted transcripts and performed gene ontology annotation. in particular, our results reveal that most of the chemosensory genes have up-regulated expression in the antennae, while a few have up-regulated expression in the terminal abdominal tissues. in addition, a strikingly large number of genes, especially or genes, have sex-specific antennal expression, which is likely to reflect the sex-divergent behaviors of adult hessian flies.  mapping and transcript prediction we performed illumina sequencing and subsequently mapped the sequenced reads to the current hessian fly genome assembly to predict transcripts and analyze gene expression profiles. the sequencing yielded ca. 314 m read-pairs, of which ca. 82 m were derived from the male antennal library, 75 m from female antennae, 77 m from female terminal abdomens, and 79 m from male terminal abdomens. the illumina reads have been submitted to the sequence read archive at ncbi. the proportion of reads that successfully mapped to the hessian fly genome assembly was high in all four transcriptomes, ranging from 75% to 80%, despite strict mapping parameters that allowed for only 2 bp mismatches to account for the high sequence similarity among recently duplicated chemosensory genes. the mapping results were used by cufflinks to predict transcripts expressed in female antennae, transcripts in male antennae, in female terminal abdomen, and in male terminal abdomen. based on these transcripts, genes and isoforms were predicted. the gc content was 34-35% in all transcriptomes. global transcriptome profiling gene ontology annotation go annotation, using blast2go, was performed to compare the four transcriptomes with regards to the frequency of transcript-associated go terms. the proportion of transcript sequences with significant blast hits was similar in the four transcriptomes, i.e. 57% for female antennae, 59% for male antennae, 65% for female terminal abdomen, and 60% for male terminal abdomen. the top-hit species distribution showed that most hits were represented by sequences from the yellow-fever mosquito, followed by the southern house mosquito and the malaria mosquito. the vinegar fly came fourth despite the fact that ca. five times as many protein sequences from d. melanogaster have been deposited in the nr database as compared to sequences from a. aegypti. the top blast hits for all predicted transcripts as well as their fpkm values and genome locations are given in additional file for the female antennae, 41% of the total transcripts had assigned go terms in addition to blast hits. corresponding numbers in the other transcriptomes were 43% for male antennae, 50% for female terminal abdomen, and 44% for male terminal abdomen. the four transcriptomes turned out to be very similar, both with respect to "molecular function" and "biological process" go annotation. the only notable differences were found in the molecular function annotation where the "odorant binding" term only was present in the two antennal tissues, and the "signal receptor activity" term was almost twice as abundant in the two antennal tissues as compared to the two terminal abdominal tissues. this result is likely to reflect the -sensory role of the main sensory organ, the antenna. however, it also indicates that genes coding for proteins with functions common to all tissue types dominate in this type of analysis that provides no information on the expression levels of the individual genes themselves. pair-wise transcriptome comparisons the expression levels of the genes predicted by cufflinks were used for global transcriptome profiling. transcriptomes were compared pair-wise, i.e. male antennae vs. female antennae; male terminal abdomen vs. female terminal abdomen; male antennae vs. male terminal abdomen; and female antennae vs. female terminal abdomen. the male and female antennal transcriptomes had very similar global expression profiles, whereas male and female terminal abdomens were more different, both from each other as well as from the two antennal transcriptomes. when divergence based on fpkm values was analyzed using jensen-shannon distance, male and female antennae showed a low level of differentiation. in contrast, the distances in all other pair-wise comparisons were larger, ranging from ca. to. the global expression profile similarity between male and female antennae, and their distinctiveness from the other transcriptomes were confirmed in a principle component analysis, and in pair-wise regression analyses. the latter analyses indicated that the female and male terminal abdominal segments had an overall higher expression than the female and male antennae, respectively. the male antennae had a similar global expression profile as the female antennae, as shown by a close to 1: relationship between the expression levels. in general, higher expression was found in the female terminal abdomen as compared to the male terminal abdomen. expression of the chemosensory genes overall expression profiles based on our manual annotation of the hessian fly genome, a total of chemosensory genes from five gene families were predicted. predicted amino acid sequences of these genes are available in additional file and phylogenies will be included in the forthcoming publication of the hessian fly genome. in contrast to the global comparison where the male and female antennal transcriptomes demonstrated the most similar expression profiles, the female and male terminal abdomens had the most similar profiles when only the chemosensory-related genes were included. regression analysis demonstrated a close to 1: relationship between the expression levels of the chemosensory genes in the two terminal abdominal tissues, although the overall expression level was low in these tissues. most of the chemosensory genes had the highest expression in the female or male antennae as compared to the female or male terminal abdomens, respectively, although a few genes had higher expression in the terminal abdominal tissues. in contrast to the global analysis, the expression profiles of the male and female antennal tissues were clearly different when only the genes related to chemosensation were included in the analysis. fpkm values for the chemosensory genes in the four tissues are presented in additional file odorant receptors based on our annotation, the current assembly of the hessian fly genome contains or genes, including a gene for the conserved co-receptor orco and six pseudogenes. our transcriptome analysis showed that of the or genes were expressed in at least one of the analyzed tissues. thus, almost 20% of the or genes were not expressed. a substantially larger number of or genes were represented by reads in the female antennae as compared to in the male antennae. thus, while a large proportion of the or genes in the genome was expressed in the antennae of the female, only ca. half of them were expressed in the antennae of the male. female and male terminal abdomens expressed and or genes, respectively, including orco. thus, relatively large numbers of or genes were expressed outside of the main olfactory organ in both sexes, although the expression level was low in general. the orco gene had the highest level of expression in both sexes, followed by or for which the fpkm value in the male antennae was 90% of the fpkm value for orco. gene family a b anumbers within brackets refer to the number of annotated genes of each chemosensory gene family in the hessian fly genome. bexpression defined as genes having at least mapped read-pairs. in the pair-wise comparison between the female antennae and the female terminal abdomen, or genes had higher expression in the antenna, whereas two or genes, or and or had higher expression in the terminal abdomen. in the corresponding comparison for the male, or genes were expressed at a higher level in the antennae, whereas four or genes had higher expression in the terminal abdomen. when comparing expression in the antennae of males vs. females, or genes showed higher expression in females, whereas were more highly expressed in males. in addition, several of the or subfamilies in our dendrogram contained receptors of which the majority showed either female or male antennal up-regulation. for instance, the genes for most of the receptors within the clades formed by or86- or98- and or102- were more highly expressed in the female antennae as compared to the male antennae, whereas most of the receptor genes within the clade formed by or111- were more highly expressed in the male antennae. furthermore, of the or genes showed no expression in the antennae of either males or females, while or genes were similarly expressed in the antennae of the two sexes. finally, five or genes showed higher level of expression in the female terminal abdomen as compared to the male terminal abdomen, whereas seven or genes had higher expression in the male terminal abdomen in relation to the female terminal abdomen. gustatory receptors we annotated members of the gr gene family in the genome, including genes for the conserved carbon dioxide receptors, three putative sugar receptor genes, and three pseudogenes. however, only of the gr genes were expressed in the tissues analyzed here. these included the genes for the three putative carbon dioxide receptors and the three sugar receptors. only slightly larger numbers of gr genes were expressed in the antennae as compared to the terminal abdominal tissues. three of the gr genes had similar levels of expression in all four tissues. in the females, six gr genes had up-regulated expression in the antennae as compared to the terminal abdomen. these included the genes for the carbon dioxide receptors that had the highest expression of all the gr genes, and the putative sugar receptor gene gr in contrast, another putative sugar receptor gene, gr had up-regulated expression in the female terminal abdomen as compared to the female antennae. similar to the females, gr1- and gr had higher expression in the male antennae as compared to the male terminal abdomen, whereas gr had higher expression in the terminal abdomen. in contrast to the female, the expression of gr was also up-regulated in the male antennae as compared to terminal abdomen, whereas gr was expressed at a similar level in the male antennae and terminal abdomen. gr was also the only gr gene that was differentially expressed between the two antennal transcriptomes. finally, four gr genes, including the putative sugar receptor genes gr4- were more highly expressed in the male terminal abdomen than in the female terminal abdomen. none of the gr genes were more highly expressed in the female terminal abdomen than in the male terminal abdomen. ionotropic receptors in total, ir genes, including four pseudogenes were annotated in the hessian fly genome. of these, genes encoded "antennal" irs and encoded "divergent" irs. we found hessian fly orthologues for most of the conserved antennal ir genes identified in other diptera, including ir8a, ir21a, ir25a, ir41a, ir60a, ir64a, ir68a, ir76b, ir93a, and seven members of the ir group. however, ir31a, ir40a, and ir92a, which are found in other diptera, were not present in the hessian fly genome and could also not be identified in our transcriptome data. the remaining mdir genes were regarded as members of the divergent class of insect irs. twelve of the identified ir genes were not expressed, whereas all the antennal ir genes showed expression in at least one of the analyzed tissues. larger numbers of ir genes were expressed in the antennal tissues as compared to the terminal abdominal tissues. the gene for the broadly occurring co-receptor ir25a had the highest expression of all the mdir genes. expression of this gene was high in the antennae of both sexes and in the male terminal abdomen. expression of ir genes was higher in the female antennae as compared to the female terminal abdomen. of these genes, encoded antennal irs and four encoded divergent irs. although ir60a is classified as an "antennal" ir, its gene showed a non-significant tendency for higher expression in the female terminal abdomen than in the female antennae. apart from that, only one of the divergent ir genes was significantly up-regulated in the female terminal abdomen. in the males, ir genes had higher expression in the antennae as compared to the terminal abdomen, while three ir genes had higher expression in the terminal abdomen. similar to the females, ir116a was one of them, but also ir87b and ir87c. only one of the ir genes, ir75c, had up-regulated expression in the female antennae compared to those of the male, representing the only example of differential ir gene expression between the antennae of the two sexes. comparing the two terminal abdominal tissues, five ir genes, including both antennal and divergent receptor genes, had higher expression in the male terminal abdomen, whereas expression of two ir genes was higher in the female terminal abdomen. sensory neuron membrane proteins the hessian fly genome assembly contained orthologues of snmp and one snmp all of them were expressed in the antennae, and six of them were expressed in both terminal abdominal tissues. in general, the snmp genes had high expression, and the most highly expressed gene had a higher fpkm value than orco. in the females, snmp1a-c, e and f had up-regulated expression in the antennae as compared to the terminal abdomen, while one was up-regulated in the terminal abdomen, both in relation to the female antennae and to the male terminal abdomen. furthermore, snmp had slightly higher expression in the female terminal abdomen than in the female antennae, but the difference was not significant. the highest expression of snmp was found in the male terminal abdomen. in the males, snmp1b-f had up-regulated expression in the antennae as compared to the terminal abdomen. comparing expression of snmp genes between the antennae of the two sexes, snmp1a was greatly up-regulated in the female. in contrast, snmp1b and 1f had male-biased antennal expression, while snmp1b, 1c, and 1e were similarly expressed in the antennae of the two sexes. odorant binding proteins we annotated obp genes in the hessian fly genome assembly, including two pseudogenes. generally, obps are classified into different phylogenetic groups depending on the number of conserved cysteine residues. classic obps have six cysteine residues, whereas members of the plus-c class have additional cysteines and one characteristic proline. five of the hessian fly obps belonged to the plus-c class. all but one of the mdobp genes showed expression in at least one of the analyzed tissues. essentially all obp genes were expressed in the antennae of both sexes, while slightly lower numbers of genes were expressed in the terminal abdominal tissues. overall, the most highly expressed obp genes had much higher expression levels than the most highly expressed or or snmp genes. the plus-c class of obp genes had comparatively low expression levels in all four tissues. in the female, most of the obp genes were more highly expressed in the antennae as compared to the terminal abdomen. only one obp gene showed a higher level of expression in the female terminal abdomen than in the female antennae, but the overall expression level of this gene was comparatively low. the males were similar to the females by having obp genes with higher expression in the antennae as compared to the terminal abdomen. these included most of the genes that also were antennally up-regulated in the female. three obp genes had higher expression in the male terminal abdomen as compared to the male antennae. however, despite the fact that most obp genes had their highest expression in the two antennal tissues, their expression level in the terminal abdomen tissues was in general much higher than the expression of the other chemosensory gene families in the abdominal tissues. a few obp genes showed differential expression between the male and female antennae. only two obp genes showed higher levels of expression in the female antennae, whereas six genes had higher expression in the male antennae. nine of the obp genes showed differential expression levels between the two terminal abdominal tissues; five of them were more highly expressed in the female, while four were more highly expressed in the male. discussion chemosensory genes in the hessian fly in relation to drosophila and mosquitoes the chemosensory gene families have been extensively studied in various drosophila and mosquito species, but had not been analyzed in plant-feeding dipterans, such as gall midges. the hessian fly in the cecidomyiidae family is both phylogenetically and ecologically well separated from flies and mosquitoes. the specific habits and highly evolved relationship to its host also makes studies of the hessian fly chemosensory genes interesting for a more complete understanding of olfaction in diptera. the current hessian fly genome assembly contains genes for ors, grs, irs, snmps, and obps. in comparison to drosophila and mosquitoes, the number of hessian fly or genes is relatively large, whereas the number of genes encoding taste receptors is small, suggesting that olfaction is the dominant chemical sense in this species. this might be related to the short adult lifespan, with olfaction being a more time-efficient modality for host and mate assessment, occurring during flight at a distance rather than after landing via contact chemoreception. while clear orthologues of ors are rarely found in distantly related species, the antennal irs appear conserved across insect species. however, the hessian fly genome lacked genes for three antennal irs that are present in other diptera, suggesting that also some olfactory functions have been lost in this herbivorous species. insects generally have two snmp genes, but some species have multiple paralogues of snmp. snmp has been shown to be important for pheromone detection in drosophila. the hessian fly genome contains six snmp genes, which is the largest number identified in any insect so far. it is possible that this large number reflects the use of a multi-component sex pheromone in the hessian fly, and that the detection of the individual pheromone components requires specific or-snmp combinations. sex-specific antennal expression of chemosensory genes transcriptome profiling using rnaseq data has been shown to accurately predict gene expression levels, and has recently been used to analyze the expression of chemosensory genes in several other insect species, including d. melanogaster and the a. gambiae mosquito. here, we used a transcriptome profiling approach to analyze sex- and tissue specific expression in the hessian fly. our global transcriptome analysis that included all the ca. predicted genes indicated that the gene expression profiles of the male and female antennae were highly similar to each other, while expression profiles of the antennae were different from the terminal abdominal tissues. in addition, the go annotation also indicated very similar frequencies of transcript-associated go terms in the two antennal transcriptomes. in fact, the results from the latter analysis were also similar to previous results from go annotation of antennal transcriptomes of two bark beetle species and the moth manduca sexta, suggesting that "transcript frequencies" in terms of presence/absence are conserved across insect orders. when only the chemosensory genes were included in the analysis, the expression profile differed widely between the four tissues. the majority of the chemosensory genes had up-regulated expression levels in the antennae as compared to the terminal abdominal tissues. furthermore and in contrast to the global analysis, a clear difference in gene expression between the antennae of males and females was found. although all analyzed chemosensory gene families contained some genes with expression levels that differed between the sexes, the most striking difference was found among the or genes, for which ca. 62% of the genes with any expression showed at least a three-fold significant difference between male and female antennae, and ca. 50% of them showed at least a ten-fold difference, with the majority being female-biased. thus, the proportion of or genes with sex-biased expression was substantially larger than previously found in d. melanogaster, a species in which both sexes feed on the same food, but similar to a. gambiae, where males and females feed on entirely different food. similar to the mosquito, the sex-biased expression in the hessian fly is likely to reflect the different ecological interests of males and females. in both sexes, the gene for the co-receptor orco had the highest expression levels of all or genes. this is consistent with its broad expression and role in forming heteromeric complexes with the conventional ors, required for the odor response of the neuron and localization of the or in the cell membrane. during their ca. 18 h adult existence, male hessian flies do just one thing, namely search for and mate with virgin females, carrying enough sperm to fertilize up to females. finding females is accomplished by long-range attraction to a female-produced sex pheromone, comprised of five to seven compounds. similar to moths e.g., it is likely that the subfamily formed by the receptors with male-biased expression in the hessian fly includes the receptors for the pheromone compounds. we currently test this hypothesis using heterologous assays. interestingly, two of the mdor genes in the clade with putative pheromone receptors had high expression levels also in the females, and other receptors in this clade had moderate female expression. thus, if these ors are the pheromone receptors, it is likely that females also are able to detect some of their pheromone compounds. in contrast to the males, for which there is no evidence in cecidomyiids of attraction to host plants, adult females depend on plant chemicals to find and select host plants for their offspring. females, during their ca. 9–24 h existence, crawl from an eclosion site in the soil to a perch where they extend the abdominal segments to release a sex pheromone and, after mating, fly to find host plants for their eggs. attraction to host plants is mediated by an array of volatile chemicals and discrimination between host and non-host plants would presumably require a larger and different set of ors as compared to detection of the pheromone. thus, it is likely that the or genes with female-biased expression are involved in host plant recognition. similar to the hessian fly, most of the ors with sex-biased expression in a. gambiae were up-regulated in the female, i.e. the sex that blood-feeds and finds sites for oviposition. interestingly, the or gene expression profiles in the hessian fly antennae were often correlated to the phylogenetic position of the ors in the tree, i.e. the majority of the receptors within a phylogenetic subfamily showed male- or, in most cases, female-biased antennal expression. it is possible that ors with similar sequences detect structurally similar chemicals. thus, this result might also reflect the relative importance of host odors for females versus the pheromone compounds for males. a correlation between sex-specific expression bias and position in the phylogeny was also found for the ors in a recent study on two ant species. this is likely to represent adaptations that enable individuals to carry out the tasks that are specific to their sex or caste. the majority of the conserved antennal mdir genes had relatively high expression in the antennae, but only one of them showed a clear difference in expression level between the sexes. this suggests that in contrast to the ors, receptors in this more ancient family detect compounds that might be of general importance for insects, regardless of sex. of the antennal ir genes, mdir60a was not expressed in the antenna. similarly, benton et al. did not detect expression of ir60a in antennal samples of d. melanogaster. to our knowledge, a ligand for this receptor has not been identified yet in any species. the carbon dioxide receptor genes mdgr1- had high expression in the antennae, suggesting that this organ is responsible for carbon dioxide detection, which is similar to drosophila, but different from mosquitoes that detect this gas with the maxillary palps. only one gr gene, mdgr had male antennal-biased expression. in a phylogenetic analysis, this receptor grouped together with dmgr43a, which has been shown to monitor fructose levels in the fly brain. the reason why a tentative fructose receptor would have high expression in the antenna of male hessian fly is not clear, but this receptor could possibly be important if males do feed on nectar. since snmp is important for pheromone detection in d. melanogaster, we expected to find the highest expression levels of the mdsnmp genes in the male antennae of the hessian fly. however, only two of the six snmp orthologues present in the hessian fly genome had higher expression in the male antennae as compared to the female antennae. these two snmps were also relatively highly expressed in the female antennae. interestingly, snmp1a had >200-fold higher expression in the female antennae as compared to the male antennae. expression of snmp genes in females might, again, suggest that females also are able to detect components of the sex pheromone, or alternatively, that the role of snmp might not be restricted to detection of sex pheromones. expression of snmp in the antennae of both sexes has also been found in other dipteran and non-dipteran species. in general, the mdobp genes had the highest expression levels of all analyzed chemosensory gene families, and the majority of them were most highly expressed in the antennae. only a few obp genes were differentially expressed between the antennae of males and females, which is similar to a. gambiae and the yellow fever mosquito, aedes aegypti. similar to the other gene families, the mdobp genes that did show sex-biased expression might be important for the detection of compounds that elicit sex-specific behaviors. expression in terminal abdominal tissues all chemosensory gene families in the hessian fly contained several genes that were represented by reads also in the terminal abdominal samples. the obp gene family was the most broadly expressed, followed by the snmps, while the or gene family had the most antennal-specific expression. nevertheless, all gene families contained genes that were significantly more highly expressed in the male or female terminal abdominal segments as compared to the antennae. in the male, these included four ors, two grs, three irs, one snmp, and three obps. in the female, two or genes and one gene from each of the other families were more highly expressed in the terminal abdomen than in the antennae. these observations suggest that chemosensation is not restricted to the antennae of gall midges. female hessian flies examine the leaf surface with the tip of the abdomen prior to oviposition, and oviposition is induced by chemical cues from host plants, consistent with the expression of chemoreceptors in this tissue. similarly, the gustatory sense of d. melanogaster in not restricted to head appendages and taste sensilla are present on the legs, wings, and female genitalia. in addition, or59b in d. melanogaster is predominantly expressed in the male accessory gland, and a recent study showed that several or genes are expressed in the testes of a. gambiae and that orco is involved in chemical-induced sperm activation. these results are consistent with the expression of several mdors including mdorco in the male terminal abdominal tissue analyzed here. similar to the present study, other studies have reported that snmp and snmp as well as the obps can be expressed in variety of tissues in addition to the chemosensory organs, suggesting that these proteins can have physiological roles independent of olfaction. the snmps and obps that were expressed in the female hessian fly terminal abdominal segments where the pheromone gland is found, might be involved in the binding and transport of pheromone molecules before these are volatilized and released. the gr and divergent ir families contained the largest proportions of genes for which we did not detect expression in the present study. similarly, the five genes of the plus-c class of obps showed only low levels of expression. it is possible that these genes are expressed mainly in legs or wings, or in yet other tissues, which we did not analyze. indeed, several ir genes and plus-c obp genes showed high expression in the labella and tarsi of a. aegypti. alternatively, these and the other "non-expressed" chemosensory genes might be expressed only in the larvae, which are likely to use contact chemical cues in order to locate a suitable feeding site inside the wheat plant.  by conducting a comprehensive manual gene annotation and transcriptome analysis we found clear sex- and tissue-specific expression patterns of the chemosensory genes in the hessian fly. while many studies of these genes have been performed on diptera, especially drosophila and mosquitoes, this is the first of a fly that feeds on plants. the expression profile difference between male and female hessian fly was more pronounced than that observed in d. melanogaster, but similar to a. gambiae, a species that also displays sex-divergent chemically-mediated behaviors. we found the most extreme sex-specific expression difference among the or genes, which is likely to reflect the sex-divergent olfactory behaviors of the short-lived adults. in addition, the relatively large number of or genes in the genome in combination with the few gr genes and divergent ir genes, suggest that long-range attraction to hosts and mates is more important than contact assessment in this species, in which the adults have reduced mouthparts and feeding behavior has not been documented. we also found that several genes of each gene family had expression that was not restricted to the antennae, indicating that they serve multiple and tissue-specific physiological roles. future studies should focus on characterizing the function of the proteins encoded by these differentially expressed chemosensory genes in order to consolidate the link between the genetics and ecology of hessian fly olfaction. 
in addition, the large number of ors in the genome but the reduced set of grs and divergent irs suggest that the short-lived adults rely more on long-range olfaction than on short-range gustation. we then mapped illumina-sequenced transcriptome reads to the genome to explore gene expression in male and female antennae and terminal abdominal segments. we identified genes encoding odorant receptors, gustatory receptors, ionotropic receptors, odorant binding proteins, and sensory neuron membrane proteins in the hessian fly genome. we show that a large number of the chemosensory genes in the hessian fly genome have sex- and tissue-specific expression profiles. our findings provide the first insights into the molecular basis of chemoreception in plant-feeding flies, representing an important advance toward a more complete understanding of olfaction in diptera and its links to ecological specialization. the order diptera contains two model organisms, the vinegar fly drosophila melanogaster and the mosquito anopheles gambiae, whose chemosensory genes have been extensively studied. our results reveal that a large number of chemosensory genes have up-regulated expression in the antennae, and the expression is in many cases sex-specific. we also report up-regulated expression of some genes from all chemosensory gene families in the terminal segments of the abdomen, which play important roles in reproduction. representing a third dipteran lineage with an interesting phylogenetic position, and being ecologically distinct by feeding on plants, the hessian fly genome sequence has recently become available. among plant-feeding insects, the hessian fly is unusual in ‘reprogramming’ the plant to create a superior food and in being the target of plant resistance genes, a feature shared by plant pathogens. hessian flymayetiola destructorgene expressiongenometranscriptomeodorant receptorionotropic receptorgustatory receptorodorant binding proteinsensory neuron membrane protein chemoreception is essential for reproductive success, including detection of sex pheromone and plant-produced chemicals by males and females, respectively.
 the chemical senses of insects mediate behaviors that are closely linked to survival and reproduction. the order diptera contains two model organisms, the vinegar fly drosophila melanogaster and the mosquito anopheles gambiae, whose chemosensory genes have been extensively studied. representing a third dipteran lineage with an interesting phylogenetic position, and being ecologically distinct by feeding on plants, the hessian fly genome sequence has recently become available. among plant-feeding insects, the hessian fly is unusual in ‘reprogramming’ the plant to create a superior food and in being the target of plant resistance genes, a feature shared by plant pathogens. chemoreception is essential for reproductive success, including detection of sex pheromone and plant-produced chemicals by males and females, respectively.  we identified genes encoding odorant receptors, gustatory receptors, ionotropic receptors, odorant binding proteins, and sensory neuron membrane proteins in the hessian fly genome. we then mapped illumina-sequenced transcriptome reads to the genome to explore gene expression in male and female antennae and terminal abdominal segments. our results reveal that a large number of chemosensory genes have up-regulated expression in the antennae, and the expression is in many cases sex-specific. sex-specific expression is particularly evident among the or genes, consistent with the sex-divergent olfactory-mediated behaviors of the adults. in addition, the large number of ors in the genome but the reduced set of grs and divergent irs suggest that the short-lived adults rely more on long-range olfaction than on short-range gustation. we also report up-regulated expression of some genes from all chemosensory gene families in the terminal segments of the abdomen, which play important roles in reproduction.  we show that a large number of the chemosensory genes in the hessian fly genome have sex- and tissue-specific expression profiles. our findings provide the first insights into the molecular basis of chemoreception in plant-feeding flies, representing an important advance toward a more complete understanding of olfaction in diptera and its links to ecological specialization. hessian flymayetiola destructorgene expressiongenometranscriptomeodorant receptorionotropic receptorgustatory receptorodorant binding proteinsensory neuron membrane protein 