- Research article
- Open Access
Evidence of peripheral olfactory impairment in the domestic silkworms: insight from the comparative transcriptome and population genetics
BMC Genomics volume 19, Article number: 788 (2018)
The insect olfactory system is a highly specific and sensitive chemical detector, which plays important roles in feeding, mating and finding an appropriate oviposition site. The ecological niche of Bombyx mori has changed greatly since domestication from B. mandarina, and its olfactory response to environmental odorants clearly decreased. However, the mechanisms that result in the olfactory impairment are largely unknown.
The antennal transcriptomes were compared between the domestic and wild silkworms. Comparison of the same sex between the domestic and wild silkworms revealed 1410 and 1173 differentially expressed genes (DEGs) in males and females, respectively. To understand the olfactory impairment, we mainly focused on the olfactory-related genes. In total, 30 olfactory genes and 19 odorant-degrading enzymes (ODEs) showed differential expression in the two comparisons, in which 19 and 14 were down-regulated in the domestic silkworm, respectively. Based on population genomic data, the down-regulated odorant receptors (ORs) showed a higher ratio of unique non-synonymous polymorphisms to synonymous polymorphisms (N/S ratio) in the domestic populations than that in the wild silkworms. Furthermore, one deleterious mutation was found in OR30 of the domestic population, which was located in transmembrane helix 6 (TM6).
Our results suggested that down-regulation of the olfactory-related genes and relaxed selection might be the major reasons for olfactory impairment of the domestic silkworm reared completely indoor environment. Reversely, wild silkworm may increase expression and remove deleterious polymorphisms of olfactory-related genes to retain sensitive olfaction.
The olfactory system is mainly responsible for the sense of smell. Insects rely on a wide range of olfactory senses to locate mates, find an appropriate oviposition site as well as to avoid predators and other dangers. Insect antenna is served as an important periphery olfactory system, which contained several sensillum types, such as sensilla trichodea (medium-length and long), sensilla basiconica and sensilla coeloconica . It plays an improtant role in reception and processing of semiochemicals, and rapid inactivation of the odorants once they have conveyed information .
Odor perception is a selective and sensitive process, which is heavily dependent on the various receptors expressed on olfactory sensory neurons (OSNs) in antennae [2, 3]. Olfactory receptors include odorant receptors (ORs), ionotropic receptors (IRs), and sensory neuron membrane proteins (SNMPs). In insects, odor molecules are first captured and transported to the receptors by water-soluble extracellular proteins that are located in the fluid surrounding the sensory dendrite of antennal sensilla, including odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) . Once the message is conveyed, the chemical signal would be rapidly inactivated to prevent the accumulation of residual stimulant . Those antenna-biased odorant-degrading enzymes (ODEs) play important roles in inactivating sex pheromones and plant volatile elicitors [5, 6]. Thus, olfaction is orchestrated by a series of olfactory proteins and odorant-degrading enzymes.
Bombyx mori is one of the model organisms for the study of insect olfaction, especially in pheromone perception [7, 8] and functional diversities of odorant receptors [9, 10]. About 5000 years ago, B. mori was domesticated from the wild silkworm B. mandarina . It was found that the domestic silkworm’s response to environmental odorants was clearly decreased when compared with the wild silkmoth . Previous studies suggested that high olfactory sensitivity has metabolic costs and would be a waste of energy if the organism’s physiology were not ready for the final behavioral output . In the wild silkmoth, maintaining olfactory sensitivity should be an important evolutionary strategy for survival and reproduction. The comparison between B. mori and its wild counterpart might be interesting for studying environment-dependent olfactory adaptation.
To cope with the variable conditions, insect OSNs can adapt to their environment through olfactory plasticity . Some studies have demonstrated that reduced expression of olfactory genes can lower the sensitivity of a sensory neuron [13,14,15]. In Culex quinquefasciatus, silencing the CquiOBP1 gene showed significantly lower electrophysiological responses to known mosquito oviposition pheromones than the antennae of water-injected, control mosquitoes . In this regard, expression plasticity of olfactory genes might be an important reason for the variation of olfaction phenotype. In addition, genetic variation of olfactory genes is another mechanism for altering the sensitivity of OSNs. In Drosophila melanogaster, single-nucleotide polymorphisms (SNPs) in OBP genes contribute to individual variation in chemosensory behavior [16, 17]. Furthermore, olfactory genes may have undergone different selective pressures in domestic and wild populations . To retain sensitive olfaction, evolutionary processes observed in the wolf have had the effect of removing deleterious polymorphisms and accumulating tolerant polymorphisms in olfactory genes, when compared with dogs .
In this study, we analyzed the antennal transcriptomes of the domestic and wild silkmoths, and mainly focused on gene families that have been implicated in olfaction. Furthermore, the genomes of eight domestic silkworms and seven wild silkworm samples have been resequenced in our lab. These genomic data enable us to detect evolutionary rates and deleterious mutations of differentially expressed chemosensory genes. The objective of this study was to detect expression difference of olfactory-related genes in the antennae of the domestic and wild silkworms and to deepen understanding of how the domestic silkworms to impair olfactory sensitivity. Comparison of the indoor species and its corresponding wild species would help us understand the potential mechanisms of olfactory adaptation in wild condition.
Materials and methods
Silkworm collection and sample preparation
The domestic silkworm strain Dazao was reared on fresh mulberry leaves at 25 ± 1 °C and 75% ± 3% relative humidity (14 h lights : 10 h dark) in an indoor chamber. The wild silkworms were collected from Chongqing City, China. The larvae of the wild silkworms were reared on mulberry trees in open-air field chambers. Two days after pupation, the pupae of the domestic and wild silkworms were transferred to a same indoor chamber until dissection of adult antennae. The intact antennae of 30 single-sex individuals were dissected from virgin moths at 24 h after adult eclosion and used for one sample. Two replicates were taken for the antennae of female (W_F) and male (W_M) of the wild silkworm, female (D_F) and male (D_M) of the domestic silkworm, respectively. All the samples were preserved in RNAlater (Ambion, Austin, USA) and stored at − 80 °C for RNA isolation.
RNA sequencing and assembly
RNA isolation, library preparation, and sequencing. Total RNA of the antennae was extracted with Trizol reagent (Invitrogen, Burlington, ON, Canada) according to the manufacturer’s instructions. RNA purity was checked with the NanoPhotometer spectrophotometer (IMPLEN, Westlake Village, CA, USA). The quality of the RNA samples was checked using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Illumina mRNA sequencing libraries were run for paired-end reads sequencing on an Illumina HiSeq™ Genome Analyzer platform (Novogene, Beijing, China).
Quality control and assembly. The raw reads were filtered by removing adaptor sequences and low-quality sequences containing > 10% poly-N or > 50% of bases whose Phred quality scores ≤5 with NGS QC Toolkit v2.3.3 . The reference genome of B. mandarnia is unavailable. The reads of the whole genome and transcriptome sequencing were often mapped to the domestic silkworm reference genome [20, 21]. The index of silkworm reference genome (http://silkworm.genomics.org.cn/) was built using Bowtie version 2.2.8 . The clean reads of B. mori and B. mandarnia were aligned to the reference genome using TopHat version 2.1.1 . The uniquely mapped reads were retained and used for further analysis. Finally, transcripts were assembled by Cufflinks version 2.2.1 .
Identification of differentially expressed genes (DEGs)
Gene expression levels were estimated using FPKM (fragments per kilobase of transcript per million mapped reads) values with Cufflinks . HTSeq v0.5.4  was used to count the number of reads mapped to each gene. Identification of DEGs was performed using the DESeq package . The P-value adjusted by Benjamini-Hochberg method of 0.05 and fold-change of 3 was set as the threshold values for significant differential expression.
GO annotation and KEGG enrichment analyses of DEGs
To predict the function, all the unigenes were used to BLASTX search against the nr (non-redundant) protein database in NCBI with an E-value <1e-5. The BLASTX results were converted into functional annotations by gene ontology (GO) terms using Blast2GO software . The statistical significance of the functional GO enrichment was evaluated using a false discovery rate (FDR < 0.05). Pathways from Kyoto Encyclopedia of Genes and Genomes (KEGG) were assigned using the online KEGG Automatic Annotation Server (KAAS).
Identification of novel olfactory-related genes
All the antenna-expressed genes (FPKM ≥1 in at least one of the four antennal samples) were translated using TransDecoder (http://transdecoder.github.io) assisted by Pfam domain information. Hidden Markov Model (HMM) files were downloaded from Pfam database (http://pfam.xfam.org/), including iGluR ligand-gated ion channel (PF00060), odorant receptor (PF02949/PF13853). HMM files were used to screen the translated proteins of the antenna-expressed genes using HMMER 3.0 (E-value <1e-5). Using known olfactory-related protein sequences as queries, BLASTP was also used to search against the predicted protein database of the antenna-expressed genes (E-value <1e-5). All the significant hits were subsequently checked by BLASTP against nr database in NCBI.
The putative amino acid sequences of ORs and IRs were aligned using MUSCLE . Positions that had a high percentage of gaps (> 70%) were trimmed. The VT + G + F was selected as the most suitable model of evolution by ProtTest 3.2  based on the Akaike information criterion (AIC). Maximum-likelihood (ML) trees were reconstructed using RAxML version 8.2.12  with “PROTGAMMAVTF” implementation, four discrete rate categories, and 100 bootstrap replicates. The neighbor-joining (NJ) and maximum parsimony (MP) trees were reconstructed by MEGA X . NJ trees were reconstructed with 500 bootstrap replicates, Jones–Taylor–Thornton (JTT) model, and a gamma distribution (shape parameter = 2.21), which were identified as relatively good models by ProtTest. MP trees were carried out with default settings, namely: the subtree pruning and regrafting (SPR) algorithm, random addition of sequences with ten replicates, and bootstrap test with 500 replicates. Based on the phylogenetic trees, novel IR genes were named according to a unified nomenclature system and Croset’s method .
Population genetics and molecular evolution of the olfactory genes
In this study, eight domestic silkworm individuals and seven wild silkworms were used for whole-genome sequencing with Illumina’s HiSeq 4000 system. The domestic samples were from eight silkworm strains, 7532, S03, S02, Xianghui, HB05, Yanjinhuang, Xiaoshiwan, and Jianpuzhai. One individual per geographical location in China was captured, including Beibei District in Chongqing City, Hongya County in Sichuan Province, Anyue County in Sichuan Province, Nanchong City in Sichuan Province, Ziyang City in Sichuan Province, Wuhan City in Hubei Province, and Suzhou City in Jiangsu Province. The clean reads of each sample were mapped to the silkworm reference genome using BWA (https://sourceforge.net/projects/bio-bwa/). Picard (http://broadinstitute.github.io/picard/) was used for sorting the BAM file by reading position and removing the highly repetitive reads. The single nucleotide polymorphism (SNP) and INDEL in each sample were identified by GATK v2.7 , Samtools  and FreeBayes . The overlapped SNPs and INDELs detected by the three tools were used for further analysis. The consensus sequences of the candidate genes were extracted in the genome sequence of each sample, including the coding sequences (CDS) and gene sequences comprised of CDS, introns, and 2-Kb 3′ and 5′ flanking regions.
In the domestic and wild populations, genetic diversities (π) of the differentially expressed olfactory genes were estimated by DnaSP 5.1 . Tajima’s D, Fu, and Li’s D*, and Fu and Li’s F* tests were used to detect whether the genes were evolving neutrally using DnaSP. The ratios of non-synonymous (dN) and synonymous (dS) substitution rates were calculated by the YN00 program implemented in the PAML 4.5 package . In addition, we counted the numbers of non-synonymous (N) and synonymous (S) SNPs unique to the domestic and wild populations, respectively. And then, the population-unique N/S ratios were calculated for the differentially expressed olfactory genes.
Quantitative real-time PCR (qPCR)
The qPCR validation experiment was performed, which the method was introduced in our previous study . The primers and annealing temperature of each gene were listed in Additional file 1: Table S1. Gene expression levels were normalized against the corresponding ribosomal protein L3 (RpL3) expression levels. The relative expression level of each gene was calculated by the relative quantification (R = 2-ΔΔCt) method .
Overview of the antennal transcriptomes in the silkworms
To explore the mechanisms of olfactory adaptation, we collected the intact antennae from the domestic silkworm strain Dazao and wild silkworm. The RNA of each sample from 30 individuals was sequenced using an Illumina Genome Analyzer (II). After filtering low-quality reads and trimming adapters (Additional file 2: Table S2), the clean reads were mapped to the B. mori reference genome (Additional file 3: Table S3). All of the mapped reads were merged and assembled using Cufflinks . Totally, 22,767 unigenes were assembled from the antennal transcriptomes of the domestic and wild silkworms (Additional file 4: Table S4). To discard transcript models that had no read coverage or low coverage (FPKM < 1) in all samples, transcripts with FPKM ≥1 in at least one sample were considered for expression. The number of genes with FPKM ≥1 varied from 13,080 to 14,022 among the four samples (Additional file 5: Table S5). KEGG assignments were used to classify the functions of the antenna-expressed genes. Relatively, signal transduction has a higher gene number (n = 550) than the other pathways (Additional file 6: Figure S1). In addition, we found some pathways related to olfactory and detoxification functions, such as the sensory system (n = 56), xenobiotics biodegradation and metabolism (74), and environmental adaptation (n = 61). In the antennae, we found that 61–74 genes were highly expressed (FPKM > 1000), in which 46 genes were overlapped among the four samples (Additional file 5: Table S5). In addition, six genes (BmGOBP1, BmGOBP2, BmOBP27, BmPBP1, BmOBP20, and BmCSP1) were extremely highly expressed (FPKM > 10,000) in all the four samples (Additional file 4: Table S4). These results were in accord with the olfactory function of insect antennae.
Identification and expression of olfactory-related genes in the antennal transcriptomes
In the silkworm, olfactory-related gene families have been widely annotated in the whole genome. Using the known gene sequences, BLAST and HMMER searches were used to identify novel olfactory-related genes in the assembled unigenes of the antennal transcriptomes (Table 1). Seven novel OR genes were characterized and assigned new names numbering from 74 upwards to 80 (Fig. 1, Additional file 7: Table S6). Ionotropic receptor (IR) is a new family of olfactory receptors and a variant subfamily of iGluRs . In the previous study, 18 IR members were identified in the silkworm genome . In this study, six novel ionotropic receptors were characterized (Fig. 2, Additional file 7: Table S6). In insects, some important gene families related to odorant degrading function were revealed, such as carboxylesterases (COEs), cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), aldehyde oxidases (AOXs), and UDP-glucuronosyltransferases (UGTs) . These important families were detected in this study, and only few novel genes were identified for COEs (1) and P450s (1) (Table 1).
Our results indicated that more than half of the whole-genome OR, IR, and OBP genes were expressed (FPKM ≥1) in the antennae of the domestic and wild silkworms (Table 1). Two sensory neuron membrane protein (SNMP) genes were identified in the domestic silkworm genome , and both of them were expressed in the adult antennae (Table 1). In addition, almost all of the CSPs showed expression signals in the antennae. For the five potential ODE families, nearly three-quarters of the GSTs, almost half of the UGTs, COEs, and P450s and 4 out of 6 AOXs were shown to be expressed. Within the ten olfactory-related families of interest here, more than 60% of the antennal-expressed OBPs, CSPs and AOXs showed high expressions (FPKM > 100) (Additional file 4: Table S4). In contrast, about 60% of the antenna-expressed IRs, ORs, and COEs and about 40% of P450 and UGTs genes presented low expressions (FPKM < 10) in all the four samples.
Identification of DEGs and primary candidates related to olfactory impairment
For a more global view, the differentially expressed genes of all comparisons were identified by DESeq . From the total of 22,767 unigenes, we identified 2197 DEGs among the four samples (Additional file 8: Table S7). The numbers of DEGs among the six comparisons are shown in Fig. 3a and b. Based on FPKM values, hierarchical clustering of all the DEGs was conducted, which indicated that more genes were up-regulated in the domestic silkworm (Fig. 3c). To validate the RNA-seq data, qPCR was performed for 13 DEGs related to olfaction. The results of qPCR and Illumina sequencing data were consistent with each other (Fig. 4). In order to understand the functions of the DEGs, GO enrichment analysis was performed in BLAST2GO . We found that some GO terms were related to olfaction (Additional file 9: Figure S2; Additional file 10: Table S8), such as sensory perception (GO:0007600; n = 49) and odorant binding (GO:0005549; n = 58).
To understand the olfactory impairment, comparison of the same sex between the domestic and wild silkworms may be more meaningful. Totally, 1410 and 1173 DEGs were identified in D_M vs. W_M and D_F vs. W_F, respectively (Additional file 8: Table S7). GO enrichment analysis suggested that metabolic process and catalytic activity categories represent the highest numbers within the DEGs (Additional file 10: Table S8). Through the functional annotations based on GO and BLAST homology searches, except for the olfactory-related genes, no more interesting gene types were found. Thus, we only focused on the olfactory-related genes in further analysis. Compared with the domestic silkworms, 13 out of 19 olfactory genes and 13 out of 15 putative ODEs were up-regulated in females of the wild silkworm (Fig. 5a), 14 out of 24 olfactory genes and 12 out of 17 ODEs were up-regulated in males (Fig. 5b). In the comparisons of the same sex between the domestic and wild silkworms, overall 30 olfactory genes and 19 ODEs were differentially expressed, in which 19 olfactory genes and 14 ODEs were up-regulated in the wild silkworm. This indicated that the decreased expressions of the olfactory-related genes may lower the sensitivity of OSNs in the domestic silkworms [13, 14].
Molecular population genetics of the olfactory genes down-regulated in the domestic silkworm
Compared with the wild silkworm, the down-regulated olfactory genes (n = 19) might be the primary candidates for impairing olfactory sensitivity to plant volatiles in B. mori, in which 18 genes belonged to ORs and OBPs (Fig. 5). To further understand olfactory impairment, molecular evolution of the down-regulated ORs and OBPs was investigated. Based on the resequencing genome data of eight domestic silkworm strains and seven wild silkworm samples in our lab, gene sequences were obtained for four OBPs and 13 ORs except for pseudogene OR48 . For the coding sequences (CDS) and whole gene sequences, almost all the genes in the domestic silkworm strains showed lower nucleotide diversities (π) than the wild populations (Table 2). This is in accord with B. mori representing a bottlenecked subpopulation of B. mandarina. To identify the type of selection that may have occurred in the different groups (domestic, wild silkworms), Tajima’s D, Fu, and Li’s D*, and Fu and Li’s F* tests were performed. In the wild population, negative values for almost all of the genes indicated that there were greater numbers of rare alleles. In the domestic strains, most of the differentially expressed olfactory genes contained positive values from those tests, which may be due to a population bottleneck or balancing selection.
Evolutionary rate and intolerant mutations of the olfactory genes
Non-synonymous nucleotide changes do alter protein sequences and could be subject to adaptive selection or relaxation, whereas synonymous changes are more or less neutral. Totally, we identified 463 and 52 synonymous/non-synonymous polymorphism sites for the 13 ORs and 4 OBPs within 15 genome samples, respectively (Table 3). The global dN/dS ratios were estimated by the YN00 program in the PAML . The results indicated that the dN/dS ratios of all the 17 olfactory genes were < 1 in the domestic and wild populations (Table 2), suggesting that these genes might have undergone purifying selection. Interestingly, the dN/dS ratios of ORs and OBPs were different from its corresponding genetic diversities (Table 2), which were comparable in the domestic and wild silkworms (Fig. 6). These results suggested that relaxed purifying selection may be driving the increase of non-synonymous evolutionary rate in the domestic silkworm.
The ratio of the numbers of non-synonymous SNPs to the numbers of synonymous SNPs (N/S) can be used to assess the relaxation of purifying selection . We counted the population-unique non-synonymous (N) and synonymous (S) SNPs for the differentially expressed olfactory genes (Table 3). The N/S ratio for ORs in the domestic silkworm (0.82) was significantly higher than that in the wild silkworm (0.46) (P < 0.01, Fisher’s exact test one-tailed). The global dN/dS ratios and population-unique N/S ratios may not be enough to detect functional relaxation in the domestic silkworm. To assess the effect of a coding non-synonymous variant, we used ‘Sorting Tolerant From Intolerant’ (SIFT) algorithm  to predict whether population-unique non-synonymous SNP in the OR genes are tolerant or intolerant. It was indicated that one deleterious amino acid mutation was found in OR30 of the domestic population (Fig. 7a). The topology of the OR proteins was predicted using HMMTOP 2.1 (http://www.sacs.ucsf.edu/cgi-bin/hmmtop.py). The deleterious mutation of OR30 was located in the transmembrane helix 6 (TM6). Compared with B. mandarina, functional redundancy of olfaction might have the effect of increasing population-unique N/S ratio and accumulating deleterious polymorphisms in the domestic silkworm.
Due to long-term exposure to different environments, the domestic silkworm and its wild relative B. mandarina showed different olfactory sensitivity to the plant odorants . Comparative study of antennal transcriptomes in indoor and wild species may help us understand the potential mechanisms for impairing olfactory sensitivity under domestication. In the adult antennal transcriptomes of the domestic and wild silkworms, 22,767 unigenes were assembled, in which over 57.45% of the unigenes were expressed (FPKM ≥1) (Additional file 4: Table S4). In insects, olfactory genes and odorant degrading enzymes mediate the olfactory response [2, 41]. In this study, some novel olfactory-related genes were characterized, especially, ORs and IRs (Table 1). We found that more than half of the olfactory-related genes were expressed in the antennae of the domestic and wild silkworms (Table 1, Additional file 4: Table S4). These expressed olfactory-related genes may be involved in processes that are directly or indirectly connected to sensory perception.
During sexual reproduction, female-produced sex pheromones guide flying males to their mates. Some function-specific odorant receptors and binding proteins provide vital links between sex pheromones and sexual behavior. In the silkworm, BmOR1, BmOR3 and BmPBP1 are the most important sex pheromone receptors and binding proteins [7, 8, 45], which showed male-biased or male-abundant expression patterns (Additional file 11: Table S9). Interestingly, the three genes have similar expression levels and no significant differences in males between the domestic and wild silkworms (Additional file 11: Table S9). Previous studies indicated that the sex pheromones bombykol and bombykal evoked odorant-specific electroantennography (EAG) amplitudes were comparable between the domestic and wild silkworm males . Our result is in accord with the neurophysiological consequences, suggesting that pheromone communication is equally important for indoor and wild species. This may be because sexual behavior might have been under strong stabilizing selection to maintain population reproduction both in domestic species, B. mori, and wild silkworm .
As a part of olfaction response, odorant-degrading enzymes are crucial to inactivation of stimulus molecules to avoid the continuous stimulation of the receptors [2, 41, 46]. For instance, the antennal GSTs can modify trans-2-hexenal, a plant-derived green leaf aldehyde known to stimulate the olfactory system of Manduca sexta . It was indicated that an antennal COE, esterase-6, showed a wide range of functions in degrading for many bioactive food esters in D. melanogaster [6, 48]. Previous studies have found that a number of putative ODEs were expressed in the antennae of the D. melanogaster , Spodoptera littoralis , Holotrichia parallela , and etc. In this study, 127 putative ODEs were expressed in the antennae of the domestic and wild silkworms (Additional file 10: Table S4). Furthermore, we identified 19 differentially expressed ODEs in D_M vs. W_M and D_F vs. W_F, of which 14 ODEs were down-regulated in the domestic silkworm (Fig. 5). Although the function of those differentially expressed ODEs has not been validated in the silkworm yet, they may have a similar function in degrading plant odorants with other species. Due to down-regulation of a great deal of the differentially expressed ODEs, it might affect the inactivation efficiency of plant stimulus molecules, and decrease the perception in domestic silkmoth . Compared with the domestic silkworm, the wild silkmoth would encounter more complex volatiles in the field, such as non-host plant odorants. The up-regulation of 14 out of 19 differentially expressed ODEs may play roles in inactivate extra volatiles and help wild silkmoth maintain its olfactory sensitivities for locating mates and an appropriate oviposition site. In addition, insects are often exposed to direct and residual contacts with toxic volatiles, especially insecticides, which may impair olfactory performances involved in scent recognition and neural treatment [41, 51]. Odorant-degrading enzymes may play important roles in detoxification of volatile xenobiotics [2, 41, 52]. Relatively, wild silkmoth would encounter much more toxic volatiles than indoor B. mori. It was suggested that the up-regulated ODEs might help wild silkworm detoxify toxic odorants to keep sensory processing sensitivity too.
In order to retain olfactory sensitivity, a proper expression level of olfactory genes, such as ORs and OBPs, is a key mechanism. In Dendroctonus armandi, RNA interference (RNAi) was used to reduce the expression of an odorant receptor . It was found that antennae of RNAi-treated D. armandi showed significantly lower electrophysiological responses to 11 major volatiles of its host . Similar studies for OBPs were conducted in D. melanogaster and Culex quinquefasciatus [13, 14]. Because the ecological niche of B. mori has changed greatly since domestication from B. mandarina. B. mori has become entirely dependent on humans for survival and does not need to seek oviposition sites and food in a complex environment. B. mandarina lives in an outdoor environment and need sensitive olfaction for mating and finding oviposition sites . In this study, the pupae of the domestic and wild silkworms were put in an identical indoor chamber until dissection of adult antennae. The differential expressions of the candidate genes might be caused by the divergence of the genomic background during long-term adaptation. Through the comparisons of the same sex between the domestic and wild silkworms, 30 differentially expressed olfactory genes were identified (Fig. 5). In addition, seven CSPs (CSP3, CSP8, CSP10, CSP11, CSP12, CSP18, and CSP19) were also identified as DEGs (Additional file 8: Table S7). Unexceptionally, all of them were up-regulated the domestic silkworm. Based on the microarray data of the tissues and developmental stages in the silkworms [54, 55], CSP3, CSP8, CSP11, CSP12 showed high expressions in various tissues and developmental stages (Additional file 12: Figure S3). We suspect that the up-regulated expressions of CSPs were not the reason for impairing olfactory sensitivity in the domestic silkworm, while its roles need to be explored in the future. For the 30 differentially expressed olfactory genes, 19 were down-regulated in the domestic silkworm, and almost all of them were ORs and OBPs (Fig. 5). In the silkworm, some of the ORs have been characterized for ligand responsiveness [9, 10]. The previous study indicated that BmOR56 showed a high sensitivity to cis-jasmone, which is a potent attractant in mulberry leaves for silkworm . In the domestic silkworm, it might be no need for highly sensitive olfaction to plant volatiles, which may result in down-regulation of BmOR56 and the other olfactory genes. However, the high background expression levels of the olfactory genes should be beneficial for wild silkworm to maintain a high perception of plant volatiles.
Except for expression level of the olfactory genes, sequence polymorphisms would also contribute to olfactory sensitivity [16, 18]. Nineteen olfactory genes were down-regulated in the domestic silkworm, of which 18 were ORs and OBPs (Fig. 5). We focused on the 17 functional OBPs and ORs excluded the pseudogene OR48 . Molecular population parameters and functional constraints were examined in the domestic and wild populations. In coding regions, the nucleotide diversities of the 17 olfactory genes showed lower nucleotide diversity in the domestic silkworms than wild silkworms (Table 2). This reduction in nucleotide diversity is likely due to inbreeding or the bottleneck experienced by domesticated strains . Neutral test on the olfactory genes displayed negative values appeared in the tested wild population (Table 2), which may be caused by positive selection or negative selection . Thus, natural selection may play important roles in maintaining olfactory sensitivity for host plant volatiles in the field.
Due to inbreeding and bottleneck effect, genetic diversities (π) of the 17 differentially expressed olfactory genes were lower in the domestic silkworm (Table 2). However, the global dN/dS ratios of the differentially expressed olfactory genes were comparable between the domestic and wild silkworms (Fig. 6). Thus, the relaxed purifying selection may drive an increase of non-synonymous evolutionary rate in the domestic silkworms. To refine this result, we found that the domestic silkworm has a higher population-unique non-synonymous/synonymous (N/S) ratio for the ORs than the wild silkworms (Table 3). Furthermore, we focused on the population-unique non-synonymous sites of the 13 OR genes. The potentially tolerated/intolerant mutations were predicted by SIFT tool . In the domestic population, one deleterious mutation was found in OR30, which was located in transmembrane helix 6 (TM6) (Fig. 7). In insects, the previous studies indicated that TM5–7 region of ORs was formed as a central part of the ion-conducting channel, in which many mutations have general deleterious effects on ion channel function and ion selectivity [56, 57]. For instance, mutation of Y464 (in TM7) would affect K+ selectivity in Bombyx mori ORCO . Except for the silkworms, accumulation of deleterious and intolerant mutations in olfactory genes was also found in domestic dogs . In future, the odorant ligands of OR30 and effects of the deleterious mutations on olfaction sensitivity need to be validated. It would help us understand whether the accumulation of deleterious mutations driven by relaxed selection is one of the mechanisms for impairing olfactory sensitivity in domestic silkmoth.
Comparative analysis of the antennal transcriptomes was conducted in the domestic and wild silkworms. The differentially expressed genes related to odor perception were identified, including 30 olfactory genes and 19 ODEs, most of which were down-regulated in the domestic silkworm. Furthermore, the down-regulated ORs showed a higher population-unique N/S ratio in the domestic populations than that in the wild silkworms. Especially, one deleterious mutation was found in OR30 in the domestic populations. It was indicated that down-regulation and functional relaxation of the olfactory-related genes might impair the olfactory sensitivity to environmental odorants in domestic silkworms. This study provides insights into the molecular mechanisms of olfactory impairment in the domestic silkworms. More importantly, it may help us understand the potential mechanisms for retaining olfactory sensitivity in wild conditions.
Differentially expressed genes
Fragments Per Kilobase of transcript per Million fragments mapped
Kyoto Encyclopedia of Genes and Genomes
Bisch-Knaden S, Daimon T, Shimada T, Hansson BS, Sachse S. Anatomical and functional analysis of domestication effects on the olfactory system of the silkmoth Bombyx mori. P Roy Soc B-Biol Sci. 2014;281(1774):20132582.
Leal WS. Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu Rev Entomol. 2013;58:373–91.
Fleischer J, Pregitzer P, Breer H, Krieger J. Access to the odor world: olfactory receptors and their role for signal transduction in insects. Cell Mol Life Sci. 2018;75(3):485–508.
Zhou JJ. Odorant-binding proteins in insects. Vitam Horm. 2010;83:241–72.
Ishida Y, Leal WS. Rapid inactivation of a moth pheromone. Proc Natl Acad Sci U S A. 2005;102(39):14075–9.
Chertemps T, Younus F, Steiner C, Durand N, Coppin CW, Pandey G, et al. An antennal carboxylesterase from Drosophila melanogaster, esterase 6, is a candidate odorant-degrading enzyme toward food odorants. Front Physiol. 2015;6:315.
Nakagawa T, Sakurai T, Nishioka T, Touhara K. Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science. 2005;307(5715):1638–42.
Sakurai T, Nakagawa T, Mitsuno H, Mori H, Endo Y, Tanoue S, et al. Identification and functional characterization of a sex pheromone receptor in the silkmoth Bombyx mori. Proc Natl Acad Sci U S A. 2004;101(47):16653–8.
Tanaka K, Uda Y, Ono Y, Nakagawa T, Suwa M, Yamaoka R, et al. Highly selective tuning of a silkworm olfactory receptor to a key mulberry leaf volatile. Curr Biol. 2009;19(11):881–90.
Anderson AR, Wanner KW, Trowell SC, Warr CG, Jaquin-Joly E, Zagatti P, et al. Molecular basis of female-specific odorant responses in Bombyx mori. Insect Biochem Mol Biol. 2009;39(3):189–97.
Sun W, Yu H, Shen Y, Banno Y, Xiang Z, Zhang Z. Phylogeny and evolutionary history of the silkworm. Sci China Life Sci. 2012;55(6):483–96.
Gadenne C, Barrozo RB, Anton S. Plasticity in insect olfaction: to smell or not to smell? Annu Rev Entomol. 2016;61:317–33.
Swarup S, Williams TI, Anholt RRH. Functional dissection of odorant binding protein genes in Drosophila melanogaster. Genes Brain Behav. 2011;10(6):648–57.
Pelletier J, Guidolin A, Syed Z, Cornel AJ, Leal WS. Knockdown of a mosquito odorant-binding protein involved in the sensitive detection of oviposition attractants. J Chem Ecol. 2010;36(3):245–8.
Liu Q, Liu W, Zeng B, Wang G, Hao D, Huang Y. Deletion of the Bombyx mori odorant receptor co-receptor (BmOrco) impairs olfactory sensitivity in silkworms. Insect Biochem Mol Biol. 2017;86:58–67.
Wang P, Lyman RF, Mackay TFC, Anholt RRH. Natural variation in odorant recognition among odorant-binding proteins in Drosophila melanogaster. Genetics. 2010;184(3):759–67.
Wang P, Lyman RF, Shabalina SA, Mackay TFC, Anholt RRH. Association of polymorphisms in odorant-binding protein genes with variation in olfactory response to benzaldehyde in Drosophila. Genetics. 2007;177(3):1655–65.
Chen R, Irwin DM, Zhang YP. Differences in selection drive olfactory receptor genes in different directions in dogs and wolf. Mol Biol Evol. 2012;29(11):3475–84.
Patel RK, Jain M. NGS QC toolkit: a toolkit for quality control of next generation sequencing data. PloS one. 2012;7(2):e30619.
Xia Q, Guo Y, Zhang Z, Li D, Xuan Z, Li Z, et al. Complete resequencing of 40 genomes reveals domestication events and genes in silkworm (Bombyx). Science. 2009;326(5951):433–6.
Fang SM, Hu BL, Zhou QZ, Yu QY, Zhang Z. Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms. BMC Genomics. 2015;16:60.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011;27(17):2325–9.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27(8):1164–5.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, et al. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010;6(8):e1001064.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nature Genet. 2011;43(5):491–8.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv preprint arXiv:12073907 [q-bioGN]; 2012.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
Yu QY, Fang SM, Zhang Z, Jiggins CD. The transcriptome response of Heliconius melpomene larvae to a novel host plant. Mol Ecol. 2016;25(19):4850–65.
Vogt RG, Grosse-Wilde E, Zhou JJ. The lepidoptera odorant binding protein gene family: gene gain and loss within the GOBP/PBP complex of moths and butterflies. Insect Biochem Mol Biol. 2015;62:142–53.
Koenig C, Hirsh A, Bucks S, Klinner C, Vogel H, Shukla A, et al. A reference gene set for chemosensory receptor genes of Manduca sexta. Insect Biochem Mol Biol. 2015;66:51–63.
Engsontia P, Sangket U, Chotigeat W, Satasook C. Molecular evolution of the odorant and gustatory receptor genes in lepidopteran insects: implications for their adaptation and speciation. J Mol Evol. 2014;79(1-2):21–39.
Vogt RG, Miller NE, Litvack R, Fandino RA, Sparks J, Staples J, et al. The insect SNMP gene family. Insect Biochem Mol Biol. 2009;39(7):448–56.
Kulmuni J, Havukainen H. Insights into the evolution of the CSP gene family through the integration of evolutionary analysis and comparative protein modeling. PloS one. 2013;8(5):e63688.
Yu Q, Lu C, Li B, Fang S, Zuo W, Dai F, et al. Identification, genomic organization and expression pattern of glutathione S-transferase in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2008;38(12):1158–64.
Ahn SJ, Vogel H, Heckel DG. Comparative analysis of the UDP-glycosyltransferase multigene family in insects. Insect Biochem Mol Biol. 2012;42(2):133–47.
Yu QY, Lu C, Li WL, Xiang ZH, Zhang Z. Annotation and expression of carboxylesterases in the silkworm, Bombyx mori. BMC Genomics. 2009;10:553.
Ai J, Zhu Y, Duan J, Yu Q, Zhang G, Wan F, Xiang ZH. Genome-wide analysis of cytochrome P450 monooxygenase genes in the silkworm, Bombyx mori. Gene. 2011;480(1-2):42–50.
Pelletier J, Bozzolan F, Solvar M, Francois MC, Jacquin-Joly E, Maibeche-Coisne M. Identification of candidate aldehyde oxidases from the silkworm Bombyx mori potentially involved in antennal pheromone degradation. Gene. 2007;404(1-2):31–40.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001;25(4):402–8.
Benton R, Vannice KS, Gomez-Diaz C, Vosshall LB. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell. 2009;136(1):149–62.
Vogt RG. Molecular basis of pheromone detection in insects. In: Gilbert LI, Iatrou K, Gill SS, editors. Comprehensive insect physiology, biochemistry, pharmacology and molecular biology. London: Elsevier; 2005. p. 753–804.
Wang MS, Zhang RW, Su LY, Li Y, Peng MS, Liu HQ, et al. Positive selection rather than relaxation of functional constraint drives the evolution of vision during chicken domestication. Cell Res. 2016;26(5):556–73.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nature Protoc. 2009;4(7):1073–81.
Krieger J, von Nickisch-Rosenegk E, Mameli M, Pelosi P, Breer H. Binding proteins from the antennae of Bombyx mori. Insect Biochem Mol Biol. 1996;26(3):297–307.
Younus F, Chertemps T, Pearce SL, Pandey G, Bozzolan F, Coppin CW, et al. Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster. Insect Biochem Mol Biol. 2014;53:30–43.
Rogers ME, Jani MK, Vogt RG. An olfactory-specific glutathione-S-transferase in the sphinx moth Manduca sexta. J Exp Biol. 1999;202:1625–37.
Younus F, Fraser NJ, Coppin CW, Liu JW, Correy GJ, Chertemps T, et al. Molecular basis for the behavioral effects of the odorant degrading enzyme Esterase 6 in Drosophila. Sci Rep. 2017;7:46188.
Belzunces LP, Tchamitchian S, Brunet J. Neural effects of insecticides in the honey bee. Apidologie. 2012;43:348–70.
Durand N, Carot-Sans G, Chertemps T, Montagne N, Jacquin-Joly E, Debernard S, et al. A diversity of putative carboxylesterases are expressed in the antennae of the noctuid moth Spodoptera littoralis. Insect Mol Biol. 2010;19(1):87–97.
Wang S, Liu Y, Zhou JJ, Yi JK, Pan Y, Wang J, et al. Identification and tissue expression profiling of candidate UDP-glycosyltransferase genes expressed in Holotrichia parallela motschulsky antennae. Bull Entomol Res. 2018;5:1–10.
Li X, Schuler MA, Berenbaum MR. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol. 2007;52:231–53.
Zhang R, Gao G, Chen H. Silencing of the olfactory co-receptor gene in Dendroctonus armandi leads to EAG response declining to major host volatiles. Sci Rep. 2016;6:23136.
Xia Q, Cheng D, Duan J, Wang G, Cheng T, Zha X, et al. Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori. Genome Biol. 2007;8(8):R162.
Wang GH, Jiang L, Zhu L, Cheng TC, Niu WH, Yan YF, et al. Characterization of Argonaute family members in the silkworm, Bombyx mori. Insect Sci. 2013;20(1):78–91.
Hopf TA, Morinaga S, Ihara S, Touhara K, Marks DS, Benton R. Amino acid coevolution reveals three-dimensional structure and functional domains of insect odorant receptors. Nat Commun. 2015;6:6077.
Nakagawa T, Pellegrino M, Sato K, Vosshall LB, Touhara K. Amino acid residues contributing to function of the heteromeric insect olfactory receptor complex. PloS one. 2012;7(3):e32372.
We thank Kathy Dibley for polishing the writing of this article, and Hong-Bo Zhang for preparing the genomic resequencing data.
This study was supported by the Fundamental Research Funds for the Central Universities (No. 2018CDXYSM0022), Initiation Fund of China West Normal University (No. 15E022), Promoting Cooperation of Scientific Research with the Region of America and Oceania and High-level Talents Training in Ministry of Education of China (No. 20152062). The funders had no role in study design, data collection, and analysis, decision to publish, or the preparation of the manuscript.
Availability of data and materials
Illumina reads for each of the libraries were deposited in Sequence Read Archive (SRA) database under the Accession no.: PRJNA407019.
Ethics approval and consent to participate
Experiments were conducted in accordance with the protocol approved by the Institutional Animal Care and Use Committee of the Chongqing University (permit number CBE-A201607020).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primer sequences used for the qPCR validation experiment. (DOCX 15 kb)
Table S2. Statistics of RNA-Seq data after quality control. R1 and R2 at the end of the sample name represent repeat 1 and 2, respectively. Gb: Giga base; Q20: percentage of bases with a Phred value of at least 20. Q30: percentage of bases with a Phred value of at least 30. (DOCX 16 kb)
Table S3. Summary of the clean reads mapped to the silkworm reference genome. (XLSX 11 kb)
Table S4. Expression levels of all genes across the four antennal samples. The accession number with “BGIBMGA” was from SilkDB (http://www.silkdb.org/silkdb/), and the others from GenBank in NCBI. A BLAST search was conducted using the all unique transcripts as query sequences against the protein non-redundant (nr) database in NCBI. The best BLAST hits and E-values are listed. (XLSX 2319 kb)
Table S5. Distribution of gene expressions in adult antennae of the domestic and wild silkworms. FPKM: Fragments Per Kilobase of transcript per Million fragments mapped. The female (W_F) and male (W_M) of the wild silkworm, female (D_F) and male (D_M) of the domestic silkworm were showed. (DOCX 15 kb)
Figure S1. Annotation of KEGG pathway for all the 22,767 unigenes. All the pathways were included in six categories: Metabolism, Genetic Information Processing, Environmental Information Processing, Cellular Processes, Organismal Systems, and Human Diseases. (PDF 232 kb)
Table S6. The sequences of the olfactory-related novel genes. (DOCX 28 kb)
Table S7. Differentially expressed genes between any two antennal samples. The P-values were adjusted for multiple testing using the Benjamini-Hochberg method. An adjusted P-value (padj) of 0.05 and fold-change of 3 was set as the threshold for significant differential expression. (XLSX 559 kb)
Figure S2. Scatterplot of enriched GO terms related to olfaction for all the differentially expressed genes. (PDF 156 kb)
Table S8. Gene Ontology enrichment analysis of the DEGs. D_Down and D_Up mean down- and up-regulated in the domestic silkworms, respectively. (XLSX 351 kb)
Table S9. The sex-biased genes in the domestic and wild silkworms. The putative functions were from annotations through BLAST against nr database in NCBI, and the detailed information of BLAST best hits was in Additional file 4: Table S4, Supporting information. (DOCX 16 kb)
Figure S3. Expression profiles of the differentially expressed CSPs in the tissues and developmental stages of the silkworm. The microarrays of tissues and development stages in the silkworm were retrieved from the previous studies [54, 55]. The expression signals were used to plot. When the expression signal was higher than 400, it was considered that this gene has expression evidence. (PDF 347 kb)
About this article
Cite this article
Qiu, C., Zhou, Q., Liu, T. et al. Evidence of peripheral olfactory impairment in the domestic silkworms: insight from the comparative transcriptome and population genetics. BMC Genomics 19, 788 (2018) doi:10.1186/s12864-018-5172-1
- Olfactory adaptation
- Functional constraint