- Research article
- Open Access
The Odorant Binding Protein Gene Family from the Genome of Silkworm, Bombyx mori
BMC Genomics volume 10, Article number: 332 (2009)
Chemosensory systems play key roles in the survival and reproductive success of insects. Insect chemoreception is mediated by two large and diverse gene superfamilies, chemoreceptors and odorant binding proteins (OBPs). OBPs are believed to transport hydrophobic odorants from the environment to the olfactory receptors.
We identified a family of OBP-like genes in the silkworm genome and characterized their expression using oligonucleotide microarrays. A total of forty-four OBP genes were annotated, a number comparable to the 57 OBPs known from Anopheles gambiae and 51 from Drosophila melanogaster. As seen in other fully sequenced insect genomes, most silkworm OBP genes are present in large clusters. We defined six subfamilies of OBPs, each of which shows lineage-specific expansion and diversification. EST data and OBP expression profiles from multiple larvae tissues of day three fifth instars demonstrated that many OBPs are expressed in chemosensory-specific tissues although some OBPs are expressed ubiquitously and others exclusively in non-chemosensory tissues. Some atypical OBPs are expressed throughout development. These results reveal that, although many OBPs are chemosensory-specific, others may have more general physiological roles.
Silkworms possess a number of OBPs genes similar to other insects. Their expression profiles suggest that many OBPs may be involved in olfaction and gustation as well as general carriers of hydrophobic molecules. The expansion of OBP gene subfamilies and sequence divergence indicate that the silkworm OBP family acquired functional diversity concurrently with functional constraints. Further investigation of the OBPs of the silkworm could give insights in the roles of OBPs in chemoreception.
Olfactory and gustatory systems play crucial roles for insects in finding food, mates, and oviposition sites. Exquisitely sensitive chemosensory systems of insects can recognize and discriminate diverse chemicals. Olfaction is mediated by specific olfactory sensory neurons, which project their dendrites into a lymphatic cavity where odorant binding proteins (OBP) are present at high concentrations. It has been suggested that OBPs have key functions in recognizing and delivering hydrophobic odorants to olfactory receptors (OR) on dendritic membranes [1–4].
Insect OBPs are a class of small water-soluble extracellular proteins with molecular masses of approximately 14 kDa . They are believed to serve multiple functions. For example, OBPs may act as solubilizers and carriers of the lipophilic odorants in sensillum lymph and also as semi-selective filters in odorant discrimination. Furthermore, OBPs may present odorants to activate neuronal receptors or to sequester and deactivate odorants after stimulation . Experimental evidence has demonstrated that OBPs could selectively bind odorants or pheromones [6–8]. Silkworm pheromone-binding protein 1 (BmPBP1) is capable of enhancing sensitivity and selectively mediating the response to bombykol rather than bombykal [9, 10]. Recently, several studies have shown that OBPs are required for correct recognition of some odors. In Drosophila, LUSH mutants are defective for avoidance of concentrated alcohols or benzaldehyde and have complete loss of sensitivity to the pheromone 11-cis vaccenyl acetate [11–13]. LUSH undergoes a pheromone-specific conformational change that triggers the firing of pheromone-sensitive neurons . Two other odorant binding proteins for Drosophila, OBP57d and OBP57e, are not only involved in taste perception but can also change the behavioral response to toxins contained in fruit . Four polymorphisms in three OBP genes in the Obp99 cluster are associated with variation in olfactory response to benzaldehyde in Drosophila . In the fire ant, Solenopsis invicta, allelic variantion in pheromone-binding protein GP-9 regulates social organization of colonies. Single-queen colonies are always homozygous for the B allele, whereas multiple-queen colonies possess at least one copy of the b variant .
Recent progress in whole-genome sequencing of insects has provided insights into the molecular mechanisms of olfaction. In the dipteran species D. melanogaster and A. gambiae, more than 50 OBPs and approximately 70 Ors have been identified [18–22]. Subsequently, 21 OBPs and 170 Ors have been annotated in the bee [23, 24]. Insects have far fewer ORs but more OBPs than do nematodes and mammals . In mammals, OBPs show fairly broad binding spectra and seem capable of acting only as general carriers . It has been proposed that odorant discrimination in insects might be due to combined usage of ORs and OBPs . Nevertheless, insect OBPs have also been found in non-chemosensory tissues, implicating them in broader non-olfactory functions [27–29].
The silkworm, Bombyx mori is an oligophagous, economically important member of the Lepidoptera, a taxonomic group that includes numerous pests of agriculture and forestry. Silkworm is a well-established model for studying insect olfaction . In addition, the olfactory system of the giant sphinx moth Manduca sexta has been investigated in depth by the Hildebrand lab . Studies of the role of OBPs in molecular olfaction in silkworm may help us to understand and manipulate olfactory-driven food preference in this and other Lepidoptera. Prior to this study, only four silkworm OBP genes had been reported .
Our annotation of forty-four candidate BmorOBP genes in the silkworm genome revealed lineage-specific subfamily expansions in the silkworm OBP family. Based on a genome-wide oligonucleotide microarray, the expression profiles of 32 candidate BmorOBP genes were detected in different tissues of day three 5th instar and through development from fifth instars to adult moths. Some BmorOBPs were specific to olfactory tissues and others were expressed broadly in non-olfactory tissues. We found several OBP genes whose expression patterns were sexually dimorphic.
The OBP gene family in the silkworm
We have identified a total of 44 candidate OBP genes in the newly assembled silkworm genome (Table 1). These include the four previously identified OBPs: PBP1 (pheromone binding protein), GOBP1 and GOBP2(general odorant binding proteins), and ABPX (antennal binding protein) . Additionally, the PBP2 and PBP3 identified herein have been deposited in GenBank (accession numbers AM403100 and AM403101). Thirty-five OBP-like genes were predicted by the GLEAN algorithm, 17 of which had their intron-exon boundaries corrected manually. Genomic regions containing another eight OBP genes were predicted by FGENESH+ and exon/intron splice site prediction software. We found many ESTs to support the prediction of 24 of the OBP genes. Because of the poor similarity of the signal peptide sequences, the first exons were difficult to predict for seven of the OBPs. OBP24 was incomplete, lacking the N-terminal region. This gene may be intact in the genome or a pseudogene. OBP9 was identified from EST data, but the entire gene is missing from the current assembly of the genome. It is possible that the current assembly does not cover the entire genome. The corresponding annotations on the Silkworm Genome Database (SilkDB) have been updated http://silkworm.swu.edu.cn/silkdb/.
OBP genes are clustered in the genome
42 OBP genes were distributed across 10 chromosomes, BmOBP28 is located on scaffold896, but cannot be mapped to a chromosome based on the current genome assembly (Figure 1, 2). More than three-quarters of the silkworm OBP genes are located within clusters as also seen in D. melanogaster, A. gambiae and A. mellifera [19, 22, 23]. Twenty OBP genes are organized into five clusters on five chromosomes. Of these gene clusters, only Cluster 5 includes intervening genes. The largest cluster (Cluster4) contains 12 OBP genes, which occur in both orientations within a 90 kb region on chromosome 18. Cluster5 (OBP1–6), including PBP1, GOBP1 and GOBP2 reported in previous studies, is located on chromosome 19 in the same orientation. Three non-OBP genes are located between OBP1 and OBP2. Cluster1 (OBP22–27) is located on chromosome 5. OBP26 is present in the reverse orientation to the other members of the cluster. Three smaller clusters (Cluster2, Cluster3 and Cluster6) containing two or four genes are present on chromosome 7, 14 and 23. This genome organization, and especially the presence of several large clusters, indicates a relatively recent expansion of the silkworm OBP family.
Characteristics of the silkworm OBP family
Insect OBPs are generally quite divergent and the overall pairwise sequence identity is modest . The alignment of the predicted silkworm OBP-like proteins (Figure 3) shows low average pairwise sequence identity between OBP family members. The predicted proteins have low molecular masses (14–22 kDa) and signal sequences are predicted at the hydrophobic N-terminus.
A six-cysteine signature is the most typical feature of classical insect OBPs . The spacing pattern of conserved cysteines in the silkworm OBP family is similar to that in Drosophila. Following the naming system proposed by Hekmat-Scafe et al. , we refer to OBPs missing C2 and C5 as Minus-C OBPs, and those carrying more than six conserved cysteine residues as Plus-C OBPs. All six cysteine residues are completely conserved among twenty-nine typical silkworm OBPs. The spacing pattern of conserved cysteines in these typical OBPs is C1-X25–68-C2-X3-C3-X31–46-C4-X8–29-C5-X8-C6 (where X is any amino acid). There are ten Minus-C OBPs (OBP22–31) that are missing the second and the fifth cysteines. Five Plus-C OBP members (OBP40–44) carry additional conserved cysteines located between C1 and C2 and after C6.
The majority of silkworm OBP genes carry 0–4 introns that are located in conserved positions (Figure 3). Most introns are inserted in phase 0 and 1. Generally, the first intron is always present in phase 0, near the cleavage site of the predicted signal peptide. Six classes of conserved splice sites have been identified in the honey bee, D. melanogaster, A. gambiae, and T castaneum (Figure 3) . The splice sites in most silkworm OBPs belong to one of the six classes. However, several genes appear to have introns inserted in nonconserved positions or phases, such as Cluster5 (Figure 3).
Phylogenetic analysis of the silkworm OBP family
The phylogenetic tree of the silkworm OBPs, constructed using the neighbor-joining method (Figure 4), indicates six possible protein subfamilies. Following the description by Hekmat-Scafe, we named these subfamilies as PBP/GOBP, CRLBP, ABPI and ABPII as well as the two atypical families Plus-C and Minus-C. High bootstrap values support many terminal relationships and three subfamilies: PBP/GOBP, Minus-C and ABPI. However there was weak support for the other three subfamilies: ABPII, Plus-C and CRLBP and the overall tree architecture. The groupings are supported by a number of additional features.
First, the spacing pattern of conserved cysteines is similar within each subfamily. The spacings of C1–C2 and C4–C5 in the PBP/GOBP and ABPI subfamilies are larger than in other subfamilies. By contrast, the spacing between C3 and C4 of Minus-C and ABPII is smaller than in the others. In all of the members of the Plus-C subfamily, C2 and C3 are separated by four residues, while C5 and C6 are separated by seven residues.
Second, the pairwise identity within each subfamily is higher than that between members of different subfamilies. The PBP/GOBP subfamily has the highest average pairwise sequence identity (36%), with a range from 22% to 55%. The average identities for the ABPII and Minus-C subfamilies are 35% and 29%, respectively. The other three subfamilies have lower internal sequence identities. Genes within the ABPI subfamily have lower average identity values than those in the ABPII subfamily.
Third, subfamilies are supported by the chromosomal clustering of OBP genes. The PBP/GOBP subfamily comprises the six members in Cluster 5. The Minus-C subfamily comprises nine members, of which seven occur in Cluster 1 and Cluster 2. OBP28, which is located on a small scaffold, shares high identity (78%) with OBP22. The Plus-C subfamily comprises five members, of which four members are located on chromosome 23. The CRLBP subfamily comprises eight members. Six of these are on chromosome 14 and OBP32–35 form Cluster 3. Gene Cluster 4 is divided into two subfamilies: the ABPI subfamily comprising seven members and the ABPII subfamily containing five genes. Two additional ABP family members are located on another scaffold. OBP20 has a single exon. In addition, we found the transcription terminating signal (AAACAAAA) in the 3' UTR. Two direct repeat sequences (TAATGAAATAAAATTA) are present in the 5'UTR and the 3'UTR. OBP20 may have moved to new genomic positions by retroposition.
Finally, all members within a subfamily share certain common intron insertion sites, which differ among subfamilies. The PBP/GOBP subfamily contains two intron insertion sites which were not found in D. melanogaster, A. gambiae, A. mellifera and T. castaneum. The ABPI and ABPII subfamilies have lost the conserved splice sites at S6 and S2, respectively. In the Minus-C subfamily, Cluster 1 and OBP28 have only one intron at the N-terminus, whereas Cluster 2 and OBP29 have additional introns at non-conserved sites. In the Plus-C subfamily, three common intron insertion sites are located at S1, S3 and S6 sites. In the CRLBP subfamily, OBP32, OBP35 and OBP36 have only one single exon each. OBP37 and OBP38 have conserved splice sites.
It is notable that the relatedness of the ABPII and Minus-C subfamilies are supported by a better bootstrap value than that for the ABPII and ABPI subfamilies. This interesting feature has also been found in bee and Drosophila. This suggests that the Minus-C subfamily may be derived from an ancestor with six conserved cysteines. To better understand the high degree of OBP sequence divergence, we analyzed the evolutionary constraints that were acting on this gene family. The average pairwise ratio of nonsynonymous to synonymous substitutions (d N/d S) for sequences in each subfamily was < 1. This indicates that there is strong negative selection for silkworm OBP genes. However, we observed that pairwise d N/d S values for several members of the Minus-C family are > 1 (see Additional file 1). This suggests that the members of the Minus-C subfamily are undergoing positive selection.
We also built a phylogenetic tree based on the alignment of OBP sequences in five species, B. mori, D. melanogaster, A. gambiae, T. castaneum and A. mellifera, representing four orders (Figure 5). The six subfamilies for the silkworm defined above also form clades. Despite little bootstrap support for the clade, the Minus-C and ABP subfamilies are grouped together with OBPs in A. mellifera and T. castaneum. Only a few orthologies could be found among these five species (e.g., BmorOBP38 having two fly orthologues). The other subfamilies show no obvious relationship across species. This suggests that significant lineage-specific expansion and divergence have occurred in these insects.
Expression patterns of the silkworm OBP genes
A numbers of EST libraries have been constructed for silkworm and more than 238,000 ESTs are available in GenBank. ESTs corresponding to 24 BmorOBP genes were identified using tBLASTn with BmorOBP protein sequences. The coding regions of 16 genes were covered completely by ESTs (Table 1). ESTs for 17 of the OBP genes were recovered broadly from chemosensory libraries, including larval maxillary galea, epidermis, brain and adult antennae. Interestingly, ESTs for seven of the OBP genes were only found in maxillary galea. OBP23 is present in silkgland and wing disk besides maxillary galea. It is noteworthy that OBP23 and OBP11 are highly expressed in maxillary galea, with 37 and 100 ESTs, respectively. Three OBP genes (OBP1–3), together with OBP14, are only observed in the antennal library. OBP27 is present in multiple tissues (silkgland, brain, malpighian tubule, fat body, midgut, wing disk and testis) as well as an antennal library. Three genes (OBP27, OBP30 and OBP31) were found in brain, with OBP31 represented by two ESTs in the compound eye. In larval epidermis, only one EST was found for OBP39. Meanwhile, ten BmorOBPs were recovered from non-chemosensory tissue libraries, such as silkgland, malpighian tubule, fat body, midgut, testis, ovary, compound eye, hemocyte, and wing disk. Most ESTs were also found in the fat body, ovary, testis, silkgland, and wing disk.
We designed and constructed a genome-wide microarray using 70-mer oligonucleotides based on the draft silkworm genome sequence database. There are four probes, each of which contains a 20-bases stretch mismatching with their target genes. The probe sw11831 is for OBP7 and OBP8, and sw20121 is for OBP25–27. The gene expression patterns of OBPs were surveyed in multiple silkworm tissues on day 3 of the fifth instar and also in whole insects at 15 different time points from day 3 of the fifth instar larva through to the adult moth. A list of the thirty-one silkworm OBP probes used in this microarray is provided in Table 1. The expression data are visualized in Figure 6.
OBP gene expression in multiple tissues on day 3 of the fifth instar silkworm is consistent with EST representation in the database. We found fifteen genes with significant levels of expression (Figure 6). The expression profiles of OBP genes differ markedly even among members of the same gene cluster. The majority of OBPs are expressed in testis, ovary, brain, epidermis and fat body. Three OBPs (OBP23, 25 and 31) gave stronger signals in brain and epidermis than in other tissues. Five genes (OBP1, 2, 13, 19, and 42) are restricted to brain and have low expression levels. OBP40 and OBP41 share a similar expression pattern in six tissues. However, overall, OBP41 is expressed at higher levels than is OBP40. OBP43 is expressed at low levels in testis, ovary, epidermis and fat body, which is a different pattern than that for OBP42. Sex-biased expression was examined based on two-fold differences in expression level between the sexes. OBP29, which is expressed at the highest level in testis and at low levels in ovary and fat body of males, was the most interesting case. However, the expression of the majority of OBPs does not appear to have an obvious sexual bias on day 3 of the fifth instar.
Our whole-organism array data failed to detect significant expression levels of 11 of the OBP genes at any time point in either sex. It is possible that some of these genes are expressed during these life stages but at levels below the detection limit. For example OBPs 13, 19 and 27, although not detected at the whole organism level in fifth instars, are detected in specific tissues at that life stage. Three classes of developmental expression profile for 18 members of the silkworm OBP-like family are shown in Figure 7. Expression of the first class of OBPs was restricted to adults and pupae near eclosion. A very faint signal was observed in male moths for OBP14. Three members of the same gene cluster (OBP1, 2, 3) and OBP13 were detected at moderate levels in adult. OBP3 was only observed in male. OBP42 was more abundant in males than females nine days after spinning. OBP30, OBP32 and OBP36 were detected with weak signals prior to the emergence and reached the highest level in male moths.
The second class of OBPs with three members (OBP31, 40, 41) is strongly expressed throughout all stages. Expression of OBP31 gradually rises and reaches its highest level in late pupae. OBP40 and OBP41 show three obvious expression peaks at larva, 60 h after spinning and adult. The expression of OBP29 in males, which is expressed at high levels in all lifes stages and reaches the highest level in adults, also follows this pattern.
Members of the third class of OBPs were expressed in several distinct phases. OBP23 is a good example with expression peaks in the larva, four days after spinning and again with a weak peak in the adult moth. The expression peaks of OBP25 are at 0 h, 60 h after spinning and in the adult. Furthermore, the transcripts were more abundant in males than in females at the late pupae stages. The highest expression level of OBP43 was at 12 h and gradually weakened until 6D after spinning. Only a weak signal was detected in the adult female. OBP8 shows a weak signal at several time points. OBP29 expression in females also followed this general pattern with high values 6d and 9d after spinning. This contrasts with its expression pattern in males, as described above.
In addition, we determined expression levels of OBP3–6 in moth antennae by QRT-PCR (Figure 8). Consistent with previous studies , BmorOBP3 was predominantly expressed in antennae of male moths. BmorOBP4 was equally expressed in antennae of both sexes. In contrast, BmorOBP5 showed substantially higher expression in female antennae. BmorOBP6 was marginally more highly expressed female antennae.
We have identified an OBP gene family comprising 44 members in the silkworm. This number is comparable to that in D. melanogaster and more than twice that in A. mellifera. Previously, eighteen chemosensory proteins (CSPs) have been identified in the silkworm genome . Although there is no conclusive experimental support for olfactory functions, it was suggested that the CSPs may be a second class of OBPs. Moderate numbers of olfactory and gustatory receptors have been reported in the silkworm genome [35–37]. Recent evidence has demonstrated that OBPs are required for recognition of odorants or pheromones in a number of species [9–17]. The silkworm OBPs could play a chemosensory role in chemoreception by a combinatorial interaction with chemoreceptors.
Our phylogenetic analysis revealed that the silkworm OBP family comprises six subfamilies. Evidence of significant expansion of, and divergence among, OBP subfamilies indicates that the OBP family has undergone rapid evolution following a complex set of gene duplication events. This may have been required to enhance the ability to detect diverse sets of odorants. For example, silkworms at the larval stage can accurately distinguish volatile compounds and tastants released from mulberry leaves and other plants. The Plus-C OBP genes have only been found in silkworm, flies, and mosquitoes; therefore, we conclude that these genes arose after the separation of the Mecopteria. The honey bee possesses the smallest OBP repertoire. The silkworm ABP and Minus-C subfamilies share a high degree of sequence similarity to their homologues in honey bee. The divergence of the OBP families between silkworm and bee might be due to differences in their social organization, foraging behaviour and life cycles. Both species find food resources using plant volatiles, however, the range of food sources exploited by the honey bee and the sophistication of its chemical communication is substantially greater than for the silkworm. OBP subfamilies that are common to both species may be required for detection of similar odorants, such as plant volatiles. The silkworm-specific OBP subfamily may be important for oviposition, mate finding and so on. For example, the PBP\GOBP subfamily forms a monophyletic group specific for Lepidoptera.
The diversity of the OBP gene family suggests a role for positive selection in the rapid evolution and functional diversification of these genes. We found evidence for positive selection in the Minus-C OBP subfamily. Nevertheless, the relatively low dN/dS value in other OBP subfamilies suggests a purifying selection due to functional constraints. This conclusion is consistent with that obtained by investigating nucleotide variation in two OBPs (OS-E and OS-F) and comparative analysis of the OBP family in 12 Drosophila genomes [38–40]. Subfamily members arising by duplication may acquire subtle functional differences. Nevertheless, small changes in sequence may have profound functional consequences . Moreover, the diversity of expression patterns and their expression as heterodimers or homodimers increase the potential function of the OBP family. Several studies have demonstrated that some OBPs might form homodimers or heterodimers [7, 42–48].
The expression profiles of silkworm OBP genes have been characterized for both sexes and several tissues and a number of developmental stages. One group of OBP genes is expressed only in olfactory tissues, whereas others are expressed more broadly, including in tissues with no known olfactory function. The majority of members in the same subfamily show a similar expression pattern. In this study, we found that only members of the PBP/GOBP subfamily are antenna-specific. OBP1 and OBP2 were found to be expressed at approximately equal levels in the antennae of female and male moths by Northern blot . The well characterized BmorPBP1 is predominantly expressed in the antennae of male moths. In vitro studies demonstrated that PBP1 could selectively bind the pheromone component bombykol [8, 9]. The presence of four candidate BmorPBPs in the silkworm genome offers the possibility of selective transport of the three female pheromone components to receptor neurons [49, 50]. This hypothesis is supported by binding studies with two PBPs from olfactory sensilla of the silkmoth species, Antheraea polyphemus and Antheraea pernyi . Female moths release bombykal to repel potential mates. Based on the sequence conservation of the PBP gene family and the uniqueness of the PBP lineage in Lepidoptera, we speculate that BmorOBP4–6 may be candidate binding proteins for bombykal. In contrast to BmorOBP3, BmorOBP5 and BmorOBP6 seem to be expressed at higher levels in female moths. The two OBPs may bind other odors, including oviposition cues or as yet uncharacterized odors that may be released by male moths .
Another typical class of OBP is the ABP subfamily, most members of which are expressed specifically in chemosensory organs. Intriguingly, we found that OBP8–12 are expressed in larval maxillary galea which contains taste sensilla. Some Drosophila OBP genes are expressed in both olfactory and gustatory tissues and some are exclusively expressed in gustatory organs [18, 53]. For example, Drosophila OBP57d and OBP57e are involved in taste perception and the response to toxins demonstrating that OBPs are important for gustation . The ABPI subfamily in silkworm might be involved in the perception of the taste of mulberry leaves.
Several members of the Minus-C and Plus-C subfamilies are expressed in multiple tissues, including several non-chemosensory organs, such as fat body, testis, and ovary. Most of these are detectable throughout development and show three obvious peaks of expression. The expression patterns characteristic of this OBP subfamily are similar to those of the CSP family in silkworm and in other species. Although pupae are dormant, the transformation from larva to moth involves profound metabolic changes. Larval organs and appendages are digested internally and replaced by adult structures. The fat body supplies considerable energy for this metamorphosis. We hypothesize that some members of the OBP and CSP families may play general physiological roles as carriers or may mediate responses to ligands that are important for metamorphosis and development.
The silkworm OBP will assist the identification of OBPs from other insect species by similarity screening, especially in the Lepidoptera. The comparison of OBP families across many different insect species may shed light on evolutionary divergence among OBPs and on insect chemosensory mechanisms of host and environmental adaptation.
44 OBPs have been identified in the genome of the silkworm. These may represent the entire repertoire of silkworm OBPs. Modest numbers of OBPs might interact with chemoreceptors to enhance the capabilities of chemoreception. The remarkable sequence divergence and subfamily expansion suggests that silkworm OBP family members bind to diverse sets of odorants. The family shows evidence of purifying selection, likely due to functional constraints. The expression profile of the OBP family suggests that these proteins might be involved in olfaction and gustation, as well as having general transport roles in non-chemosensory tissues.
Identification of the silkworm OBP family members
Known insect OBP sequences were downloaded from GenBank [19, 22, 23, 54]. These insect OBP sequences were used to search for similar genes in the silkworm genome sequence with TBLASTN [55, 56]. Silkworm genomic regions containing OBP genes were predicted using FGENESH+ . Gene prediction was revised by comparing with the EST database. Candidate OBP genes were checked for three universal features of the insect OBP family: a conserved cysteine pattern, a predicted size (~14 kDa), and a signal sequence predicted by SignalP .
Nomenclature of the silkworm OBPs
We adopted nomenclature for silkworm OBPs that is analogous to those proposed for the Anopheles OBP and the honey bee OBPs [22, 23]. We use the prefix BmorOBP to reflect that the gene is a putative member belonging to the silkworm Odorant Binding Protein-like family. The previously published silkworm general binding protein (BmorGOBP) has been renamed BmorOBP1. OBP genes organized into a cluster were given consecutive numbers. The classical OBP members are listed prior to the atypical members (Table 1).
The forty-four conceptually translated protein sequences from silkworm OBP genes identified in this study, along with the OBPs from four other insect species (D. melanogaster, A. gambiae, A. mellifera and T. castaneum), were used to construct a phylogenetic tree. Sequences were aligned using ClustalX . The α-helices identified in BmorPBP1 were indicated under the alignment . The highly divergent signal peptide sequences at the N-terminus were truncated. Neighbor-joining trees were produced using the Phylip package . Bootstrap analysis was performed using 1000 neighbor-joining replicates. The tree was displayed using MEGA4 . All pairwise dn/ds values were calculated using the program KaKs_Calculator. The program adopts model selection and model averaging to calculate nonsynonymous (Ka) and synonymous (Ks) substitution rates . In addition, several existing methods (NG, LWL, LPB, MLWL, MLPB, YN, MYN, and GY) for calculating Ka and Ks are also incorporated into KaKs_Calculator. Fisher's exact test for small samples was applied to justify the validity of Ka and Ks calculated.
The silkworm strain Dazao was reared on mulberry leaves and pupae were maintained at room temperature until eclosion. For the silkworm genome-wide oligonucleotide microarray, the anterior/median silk gland (A/MSG), posterior silk gland (PSG), testis, ovary, fat body, midgut, integument, hemocyte, malpighian tubule, and head on day 3 of the fifth instar were hand-dissected on ice. We also collected individuals at 15 different time points from day 3 of the fifth instar to moth. Moth antennae, heads, thoraces, abdomens, legs, wings were dissected and all tissues were immediately frozen with liquid nitrogen and stored at -70°C until use.
A genome-wide microarray with 22,987 70-mer oligonucleotides was designed and constructed by Southwest University (Chongqing, China) and CapitalBio Corporation (Beijing, China) . Briefly, the 70-mer probes representing 21,375 predicted genes from the silkworm WGS and 1,612 ESTs of interest that were not contained in the predicted genes were synthesized by MWG-Biotech Inc (Ebersberg, Germany). Arrays were fabricated using a SmartArrayer™ (CapitalBio Corporation, Beijing, China). cDNA labeled with fluorescent dye (Cy5 and Cy3-dCTP) was produced by the linear RNA amplification method . The resulting labeled cDNAs were denatured in hybridization solution (3×SSC, 0.2% SDS, 5×Denhardt's solution and 25% formamide) at 95°C for 3 min before hybridization. Hybridizations were performed in a hybridization chamber which was placed in a three-phase tilting agitator (BioMixer™; CapitalBio). After hybridization, slides were washed twice with washing solution. We used a dual-dye experiment to analyze the expression patterns and each experiment was performed as a dye-swap. In pilot experiments, we performed several self-to-self hybridizations to evaluate system noise. All arrays were scanned with a confocal LuxScan™ scanner and the images obtained were analyzed using LuxScan™ 3.0 software (CapitalBio Corporation, Beijing, China). For the individual channel data extracts, we removed any faint spots from both channels (Cy3 and Cy5), whose signal intensities were below 400 units after subtracting the background. The linear normalization method was used to normalize individual channel data, based on the expression levels of four confirmed housekeeping genes. The normalized signal intensity values were further analyzed using one-way analysis of variance (ANOVA), with the significance level set at a P value of less than 0.001 (P < 0.001) across all investigated tissues. We visualized the cluster data using the Treeview program .
Total RNA was isolated from adult tissues using Trizol Reagent (Invitrogen). First-strand cDNA was synthesized using the superscript first-strand synthesis system for RT-PCR (Invitrogen). Gene-specific primers for RT-PCR were synthesized commercially and silkworm actin3 was used as an internal control. For OBP3, OBP4, OBP5 and OBP6, the QRT-PCR reactions were performed using an ABI PRISM™ 7000 Sequence Detection System (Applied Biosystems) and SYBR Green I (SYBR® Premix Ex Taq™, TaKaRa). The PCR cycles were as follows: a 10 min denaturation at 95°C, followed by 40 cycles with 5s at 95°C and 30s at 60°C. The expression levels of the PBP candidate genes were calculated relative to the control gene BmActin3 according to Livak & Schmittgen .
Pelosi P: Perireceptor events in olfaction. J Neurobiol. 1996, 30: 3-19. 10.1002/(SICI)1097-4695(199605)30:1<3::AID-NEU2>3.0.CO;2-A.
Pelosi P: Odorant-binding proteins: structural aspects. Ann NY Acad Sci. 1998, 855: 281-293. 10.1111/j.1749-6632.1998.tb10584.x.
Steinbrecht RA: Odorant-binding proteins: expression and function. Ann NY Acad Sci. 1998, 855: 323-332. 10.1111/j.1749-6632.1998.tb10591.x.
Krieger J, Breer H: Olfactory reception in invertebrates. Science. 1999, 286: 720-723. 10.1126/science.286.5440.720.
Pelosi P: Odorant-binding proteins. Crit Rev Biochem Mol Biol. 1994, 29 (3): 199-228. 10.3109/10409239409086801.
Vogt RG, Riddiford LM: Pheromone binding and inactivation by moth antennae. Nature. 1981, 293: 161-163. 10.1038/293161a0.
Plettner E, Lazar J, Prestwich EG, Prestwich GD: Discrimination of pheromone enantiomers by two pheromone binding proteins from the gypsy moth Lymantria dispar. Biochemistry. 2000, 39: 8953-8962. 10.1021/bi000461x.
Pophof B: Pheromone-binding proteins contribute to the activation of olfactory receptor neurons in the silkmoths Antheraea polyphemus and Bombyx mori. Chem Senses. 2004, 29: 117-125. 10.1093/chemse/bjh012.
Grosse-Wilde E, Svatos A, Krieger JA: pheromone binding protein mediates the bombykol-induced activation of a pheromone receptor in vitro. Chem Senses. 2006, 31: 547-555. 10.1093/chemse/bjj059.
Syed Z, Ishida Y, Taylor K, Kimbrell DA, Leal WS: Pheromone reception in fruit flies expressing a moth's odorant receptor. Proc Natl Acad Sci USA. 2006, 103: 16538-16543. 10.1073/pnas.0607874103.
Kim MS, Repp A, Smith DP: LUSH odorant-binding protein mediates chemosensory responses to alcohols in Drosophila melanogaster. Genetics. 1998, 150: 711-721.
Wang Y, Wright NJD, Guo HF, Xie Z, Svoboda K, Malinow R, Smith DP, Zhong Y: Genetic manipulation of the odor-evoked distributed neural activity in the Drosophila mushroom body. Neuron. 2001, 29: 267-276. 10.1016/S0896-6273(01)00196-9.
Xu PX, Atkinson R, Jones DN, Smith DP: Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron. 2005, 45: 193-200. 10.1016/j.neuron.2004.12.031.
Laughlin JD, Ha TS, Jones DN, Smith DP: Activation of pheromone-sensitive neurons is mediated by conformational activation of Pheromone-Binding Protein. Cell. 2008, 133: 1255-1265. 10.1016/j.cell.2008.04.046.
Matsuo T, Sugaya S, Yasukawa J, Aigaki T, Fuyama Y: Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol. 2007, 5: e118-10.1371/journal.pbio.0050118.
Wang P, Lyman RF, Shabalina SA, Mackay TF, Anholt RR: Association of Polymorphisms in Odorant-Binding Protein Genes With Variation in Olfactory Response to Benzaldehyde in Drosophila. Genetics. 2007, 177: 1655-1665. 10.1534/genetics.107.079731.
Krieger MJ, Ross KG: Identification of a major gene regulating complex social behavior. Science. 2002, 295: 328-332. 10.1126/science.1065247.
Galindo K, Smith DP: A large family of divergent odorant-binding proteins expressed in gustatory and olfactory sensilla. Genetics. 2001, 159: 1059-1072.
Hekmat-Scafe DS, Scafe CR, McKinney AJ, Tanouye MA: Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002, 12: 1357-1369. 10.1101/gr.239402.
Hill CA, Fox AN, Pitts RJ, Kent LB, Tan PL, Chrystal MA, Cravchik AF, Collins H, Robertson HM, Zwiebel LJ: G-protein-coupled receptors in Anopheles gambiae. Science. 2002, 298: 176-178. 10.1126/science.1076196.
Robertson H, Warr C, Carlson J: Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc Natl Acad Sci USA. 2003, 100: 14537-14542. 10.1073/pnas.2335847100.
Xu P, Zwiebel L, Smith D: Identification of a distinct family of genes encoding atypical odorant-binding proteins in the malaria vector mosquito, Anopheles gambiae. Insect Mol Biol. 2003, 12: 549-560. 10.1046/j.1365-2583.2003.00440.x.
Forêt S, Maleszka R: Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006, 16: 1404-1413. 10.1101/gr.5075706.
Robertson HM, Wanner KW: The chemoreceptor superfamily in the honey bee Apis mellifera: Expansion of the odorant, but not gustatory, receptor family. Genome Res. 2006, 16: 1395-1403. 10.1101/gr.5057506.
Prasad B, Reed R: Chemosensation: Molecular mechanisms in worms and mammals. Trends Genet. 1999, 15: 150-153. 10.1016/S0168-9525(99)01695-9.
Löbel D, Jacob M, Volkner M, Breer H: Odorants of different chemical classes interact with distinct odorant binding protein subtypes. Chem Senses. 2002, 27: 39-44. 10.1093/chemse/27.1.39.
Paesen G, Happ G: The B proteins secreted by the tubular accessory sex glands of the male mealworm beetle, Tenebrio molitor, have sequence similarity to moth pheromone-binding proteins. Insect Biochem Mol Biol. 1995, 25: 401-408. 10.1016/0965-1748(94)00085-V.
Thymianou S, Mavroidis M, Kokolakis G, Komitopoulou K, Zacharopoulou A, Mintzas A: Cloning and characterization of a cDNA encoding a male-specific serum protein of the Mediterranean fruit fly, Ceratitis capitata, with sequence similarity to odourant-binding proteins. Insect Mol Biol. 1998, 7: 345-353. 10.1046/j.1365-2583.1998.740345.x.
Paiva-Silva G, Sorgine M, Benedetti C, Meneghini R, Almeida I, Machado E, Dansa-Petretski M, Yepiz-Plascencia G, Law J, Oliveira P, Masuda H: On the biosynthesis of Rhodnius prolixus heme-binding protein. Insect Biochem Mol Biol. 2002, 32: 1533-1541. 10.1016/S0965-1748(02)00074-7.
Butenandt A, Beckmann R, Stamm D, Hevker E: On the sex pheromone of the silkworm moth Bombyx mori. Isolation and structure. Z Naturforsch B. 1959, 14: 283-284.
Hildebrand JG, Shepherd GM: Mechanisms of olfactory discrimination: converging evidence for common principles across phyla. Annu Rev Neurosci. 1997, 20: 595-631. 10.1146/annurev.neuro.20.1.595.
Krieger J, von Nickisch-Rosenegk E, Mameli M, Pelosi P, Breer H: Binding proteins from the antennae of Bombyx mori. Insect Biochem Mol Bio. 1996, 26: 297-307. 10.1016/0965-1748(95)00096-8.
Pelosi P, Maida R: Odorant-binding proteins in insects. Comp Biochem Physiol B. 1995, 111: 503-514. 10.1016/0305-0491(95)00019-5.
Gong DP, Zhang HJ, Zhao P, Lin Y, Xia QY, Xiang ZH: Identification and expression pattern of the chemosensory protein gene family in the silkworm, Bombyx mori. Insect Biochemistry and Molecular Biology. 2007, 37: 266-277. 10.1016/j.ibmb.2006.11.012.
Wanner KW, Anderson AR, Trowell SC, Theilmann DA, Robertson HM, Newcomb RD: Female-biased expression of odourant receptor genes in the adult antennae of the silkworm, Bombyx mori. Insect Mol Biol. 2007, 16: 107-119. 10.1111/j.1365-2583.2007.00708.x.
Wanner KW, Robertson HM: The gustatory receptor family in the silkworm moth Bombyx mori is characterized by a large expansion of a single lineage of putative bitter receptors. Insect Mol Biol. 2008, 17: 621-629. 10.1111/j.1365-2583.2008.00836.x.
Nakagawa T, Sakurai T, Nishioka T, Touhara K: Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science. 2005, 307: 1638-1642. 10.1126/science.1106267.
Sánchez-Gracia A, Aguadé M, Rozas J: Patterns of nucleotide polymorphism and divergence in the odorant-binding protein genes OS-E and OS-F: analysis in the melanogaster species subgroup of Drosophila. Genetics. 2003, 165: 1279-1288.
Sánchez-Gracia A, Rozas J: Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans. Genetics. 2007, 175: 1923-1935. 10.1534/genetics.106.068015.
Vieira FG, Sánchez-Gracia A, Rozas J: Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: Purifying selection and birth-and-death evolution. Genome Biol. 2007, 8: R235-10.1186/gb-2007-8-11-r235.
Leal WS, Chen AM, Ishida Y, Chiang VP, Erickson ML, Morgan TI, Tsuruda JM: Kinetics and molecular properties of pheromone binding and release. Proc Natl Acad Sci USA. 2005, 102: 5386-5391. 10.1073/pnas.0501447102.
Campanacci V, Longhi S, Nagnan-Le MP, Cambillau C, Tegoni M: Recombinant pheromone binding protein 1 from Mamestra brassicae (MbraPBP1). Eur J Biochem. 1999, 264: 707-716. 10.1046/j.1432-1327.1999.00666.x.
Danty E, Briand L, Michard-Vanhee C, Perez V, Arnold G, Gaudemer O, Huet D, Huet JC, Ouali C, Masson C, Pernollet JC: Cloning and expression of a queen pheromone-binding protein in the honeybee: an olfactory-specific, developmentally regulated protein. J Neurosci. 1999, 19: 7468-7475.
Sandler BH, Nikonova L, Keal WS, Clardy J: Sexual attraction in the silkworm moth: structure of the pheromone-binding protein-bombykol complex. Chem Biol. 2000, 7: 143-151. 10.1016/S1074-5521(00)00078-8.
Ban L, Scaloni A, D'Ambrosio C, Zhang L, Yan Y, Pelosi P: Biochemical characterization and bacterial expression of an odorant-binding protein from Locusta migratoria. Cell Mol Life Sci. 2003, 60: 390-400. 10.1007/s000180300032.
Honson N, Johnson MA, Oliver JE, Prestwich GD, Plettner E: Structure-activity studies with pheromone-binding proteins of the gypsy moth, Lymantria dispar. Chem Senses. 2003, 28: 479-489. 10.1093/chemse/28.6.479.
Kruse SW, Zhao R, Smith DP, Jones DNM: Structure of a specific alcohol-binding site defined by the odorant binding protein LUSH from Drosophila melanogaster. Nat Struct Biol. 2003, 10: 694-700. 10.1038/nsb960.
Andronopoulou E, Labropoulou V, Douris V, Woods DF, Biessmann H, Iatrou K: Specific interactions among odorant-binding proteins of the African malaria vector Anopheles gambiae. Insect Mol Biol. 2006, 15: 797-811. 10.1111/j.1365-2583.2006.00685.x.
Kasang G, Kaissling KE, Vostrowsky O, Bestmann J: Bombykal, the second component of the Bombyx mori L. silkworm pheromone. Angew Chem. 1978, 90: 74-75. 10.1002/ange.19780900132.
Kasang G, Schneider D, Schäfer W: The silkworm moth Bombyx mori. Presence of the (E, E)-Stereoisomer of Bombykol in the female pheromone gland. Naturwissenschaften. 1978, 65: 337-338. 10.1007/BF00368376.
Maida J, Krieger T, Gebauer U, Lange U, Ziegelberger G: Three pheromone-binding proteins in olfactory sensilla of the two silkmoth species Antheraea polyphemus and Antheraea pernyi. Eur J Biochem. 2000, 267: 2899-2908. 10.1046/j.1432-1327.2000.01303.x.
Alisha RA, Wanner KW, Trowell SC, Warr CG, Jaquin-Joly E, Zagatti P, Robertson H, Newcomb RD: Molecular basis of female-specific odorant responses in Bombyx mori. Insect Biochemistry and Molecular Biology. 2009, 39: 189-197. 10.1016/j.ibmb.2008.11.002.
Shanbhag S, Park S, Pikielny C, Steinbrecht R: Gustatory organs of Drosophila melanogaster: Fine structure and expression of the putative odorant-binding protein PBPRP2. Cell Tissue Res. 2001, 304: 423-437. 10.1007/s004410100388.
Zhou J, Huang W, Zhang G, Pickett J, Field L: "Plus-C" odorant-binding protein genes in two Drosophila species and the malaria mosquito Anopheles gambiae. Gene. 2004, 327: 117-129. 10.1016/j.gene.2003.11.007.
Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, Zhao P, Zha X, Cheng T, Chai C, Pan G, Xu J, Liu C, Lin Y, Qian J, Hou Y, Wu Z, Li G, Pan M, Li C, Shen Y, Lan X, Yuan L, Li T, Xu H, Yang G, Wan Y, Zhu Y, Yu M, Shen W, Wu D, Xiang Z, Yu J, Wang J, Li R, Shi J, Li H, Li G, Su J, Wang X, Li G, Zhang Z, Wu Q, Li J, Zhang Q, Wei N, Xu J, Sun H, Dong L, Liu D, Zhao S, Zhao X, Meng Q, Lan F, Huang X, Li Y, Fang L, Li C, Li D, Sun Y, Zhang Z, Yang Z, Huang Y, Xi Y, Qi Q, He D, Huang H, Zhang X, Wang Z, Li W, Cao Y, Yu Y, Yu H, Li J, Ye J, Chen H, Zhou Y, Liu B, Wang J, Ye J, Ji H, Li S, Ni P, Zhang J, Zhang Y, Zheng H, Mao B, Wang W, Ye C, Li S, Wang J, Wong GK, Yang H: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306: 1937-1940. 10.1126/science.1102210.
Salamov A, Solovyev V: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10: 516-522. 10.1101/gr.10.4.516.
Bendtsen J, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.
Lee D, Damberger F, Peng G, Horst R, Guntert P, Nikonova L, Leal W, Wuthrich K: NMR structure of the unliganded Bombyx mori pheromone-binding protein at physiological pH. FEBS Lett. 2002, 531: 314-318. 10.1016/S0014-5793(02)03548-2.
Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.65. 2005, Department of Genome Sciences, University of Washington, Seattle, WA
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Zhang Z, Li J, Wang J, Wong GK, Yu J: KaKs_Calculator: Calculating Ka and Ks through Model Selection and Model Averaging. Geno Prot Bioinfo. 2006, 4: 259-263. 10.1016/S1672-0229(07)60007-2.
Xia QY, Cheng DJ, Duan J, Wang GH, Cheng TC, Zha XF, Liu C, Zhao P, Dai FY, Zhang Z, He NJ, Zhang L, Xiang ZH: Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori. Genome Biology. 2007, 8: R162-10.1186/gb-2007-8-8-r162.
Guo Y, Guo HY, Zhang L, Xie HY, Zhao X, Wang FX, Li Z, Wang YH, Ma SL, Tao JP, Wang WX, Zhou YX, Yang WP, Cheng J: Genomic analysis of anti-hepatitis B virus (HBV) activity by small interfering RNA and lamivudine in stable HBV-producing cells. J Virol. 2005, 79: 14392-14403. 10.1128/JVI.79.22.14392-14403.2005.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time PCR and the the 2(-Delta Delta C(T)) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This work was supported by grants from the National Basic Research Program of China (No. 2005CB121000), Program for Changjiang Scholars and Innovative Research Team in University (No. IRT0750), and the National Natural Science Foundation (No. 30571407). We thank Dr. Stephen Trowell for critically editing the manuscript.
DPG carried out the analysis of genome and microarray data and drafted the manuscript. HJZ performed microarray analysis and revised the manuscript. PZ helped to perform the experiment of QRT-PCR experiment. QYX participated in conception and design of the study, and helped write the manuscript. ZHX contributed to the interpretation of results. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Ka/Ks estimation based on model selection using KaKs_Calculator. The data represent the average pairwise ratio of nonsynonymous to synonymous substitutions (d N/d S) for sequences in six subfamily. (DOC 92 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Gong, DP., Zhang, HJ., Zhao, P. et al. The Odorant Binding Protein Gene Family from the Genome of Silkworm, Bombyx mori. BMC Genomics 10, 332 (2009). https://doi.org/10.1186/1471-2164-10-332
- Malpighian Tubule
- Wing Disk
- Male Moth
- Mulberry Leave