The oriental armyworm Mythimna separata (Walk) is a serious migratory pest; however, studies on its olfactory response and its underlying molecular mechanism are limited. To gain insights to the olfactory mechanism of migration, olfactory genes were identified using antennal transcriptome analysis. The olfactory response and the expression of olfactory genes for 1-day and 5-day-old moths were respectively investigated by EAG and RT-qPCR analyses.
Putative 126 olfactory genes were identified in M. separata, which included 43 ORs, 13 GRs, 16 IRs, 37 OBPs, 14 CSPs, and 3 SNMPs. RPKM values of IR75d and 10 ORs were larger than co-receptors IR25a and ORco, and the RPKM value of PR2 was larger than that of other ORs. Expression of GR1 (sweet receptor) was higher than that of other GRs. Several sex pheromones activated evident EAG responses where the responses of 5-day-old male moths to the sex pheromones were significantly greater than those of female and 1-day old male moths. In accordance with the EAG response, 11 pheromone genes, including 6 PRs and 5 PBPs were identified in M. separate, and the expression levels of 7 pheromone genes in 5-day-old moths were significantly higher than those of females and 1-day-old moths. PR2 and PBP2 might be used in identifying Z11-16: Ald, which is the main sex pheromone component of M. separata. EAG responses to 16 plant volatiles and the expression levels of 43 olfactory genes in 1-day-old moths were significantly greater than that observed in the 5-day-old moths. Heptanal, Z6-nonenal, and benzaldehyde might be very important floral volatiles for host searching and recognized by several olfactory genes with high expression. Some plant volatiles might be important to male moths because the EAG response to 16 plant volatiles and the expression of 43 olfactory genes were significantly larger in males than in females.
The findings of the present study show the effect of adult age on olfactory responses and expression profile of olfactory genes in the migratory pest M. separate.
The oriental armyworm, Mythimna separata (Walk), is a serious pest of rice, maize, sorghum, and wheat in China, Japan, Southeast Asia, India, Eastern Australia, New Zealand, and some Pacific Islands [1, 2]. Long-distance movement of M. separata between overwintering sites and non-overwintering sites were responsible for six nationwide economic pests in China from 1970 to 1978 . The environmental, physiological, hormonal, and genetic control of individual M. separata migratory behavior has been systematically elucidated [3–7]. However, our understanding of olfaction systems in M. separate is limited.
Olfaction plays a key role in the interaction of moths with their environment such as foraging, aggregation, mating, and oviposition behaviors. M. separata moths utilize olfaction to find nectar, which is used as a supplementary nutrient for egg-ripening  and as an energy supplement for flight . Sex pheromones have been used to identify conspecific partners of M. separata for mating and its female sex pheromone has been identified [10, 11]. In addition, host plant volatiles but not light conditions significantly affect the level of nocturnal activity in M. separata caterpillars [12, 13]. Only a few studies on M. separata olfaction have been conducted, which include the antennal sensilla type , the EAG response to sex pheromone , and the remote-sensing sex pheromone trap for real-time monitoring of M. separata .
Recently, there has been significant progress in identifying olfactory genes in moths such as Bombyx mori [17–19], Manduca sexta , Helicoverpa armigera and H. assulta , and Spodoptera frugiperda  and S. litura . To date, only three ORs (two PRs and one ORco) and one PBP were identified in M. separata antennae . To gain more insights into the molecular mechanism underlying M. separata olfaction and into the olfactory mechanism of insect migration, the present study performed antennal de novo transcriptome analysis to identify olfactory genes and RT-qPCR to compare the expression profiles between sexes and different pre-mating statuses (1-day- and 5-day-old). The link between OR gene expression and chemosensory responses as measured by electroantennography is also discussed.
Eggs of M. separata were purchased from the Chinese Academy of Agricultural Sciences and were raised in odorless insect incubators with a temperature of 25 ± 1 °C, humidity of 70 ± 7%, and photoperiod (L: D) of 14 h:10 h (from 18:00). Larvae were fed with fresh maize leaves until pupation. Adults were checked at 17:00 daily. Male and female adults were divided into different insect rearing cages (30 cm × 30 cm × 30 cm) and were fed daily with fresh 10% glucose water.
Recordings of whole-antenna electrical activity in response to chemicals were performed according to the standard technique described elsewhere . Each chemical (Additional file 1: Table S1) (10 μL) diluted with paraffin oil to concentrations of 10−2 and 10−4 (V/V) was added onto a filter paper (30 mm × 4 mm) as stimulus, and paraffin oil was used as control. Six antennae of 1-day- and 5-day-old moths were tested after exposure to each chemical for 3.5 h in the dark. As the response of antennae declined during the course of experiment, response to the 1% concentration (V/V) of (Z)-3-hexenyl acetate was used as reference and responses to all tested chemicals were standardized.
Sample collection and total RNA extraction
Antennae of 1-day- and 5-day-old 25 adult M. separata (males and females separately) were collected after 3.5 h in the dark. Samples from each group were immediately homogenized in TRNzol-A+ (TIANGEN Biotech, Beijing, China) on ice, and total RNA was extracted according to the manual. The concentration and purity of total RNA were determined by using a spectrophotometer NanoDrop2000 (ThermoFisher, Waltham, MA, USA). RNA with an A260/A280 ratio between 1.75–2.05, an A260/A230 ratio > 1, and a concentration > 400 ng/μL were used in the following experiments. Extracts were treated with DNase I (Takara, Kusatsu, Shiga, Japan) to remove any DNA. RNA extractions were replicated three times.
Analysis of de novo transcriptome
The cDNA libraries for transcriptome analysis were prepared using TruSeq SBS Kit v3-HS (Illumina, San Diego, CA, USA) following the manufacturer’s recommendations. Briefly, oligo (dT) magnetic beads were used to isolate poly (A) mRNA from collected total RNA. A fragmentation buffer was added to break the mRNA into short fragments. A random hexamer primer was used to synthesize the first-strand cDNAs using short RNA fragments as templates. The second-strand cDNAs were synthesized using a buffer, dNTPs, RNase H, and DNA polymerase I. After purification, the short cDNAs were linked to sequencing adapters. The libraries were sequenced using a paired-end transcriptome platform Illumina NextSeq500 (Illumina, San Diego, CA, USA) with 2 × 150 bp. Adapters and low-quality bases (base quality < 20) containing reads were discarded from the raw reads to obtain clean reads for analysis. De novo transcriptome assembly was conducted with the short reads assembling program, Trinity . BLASTx alignment (E value < 0.00001) between unigenes and protein databases, including NCBI non-redundant protein sequences, Gene Ontology , Kyoto Encyclopedia of Genes and Genome , eggNOG , and SWISS-PROT was successively performed.
Analysis of olfactory genes
The putative olfactory genes were obtained from gene annotation. Amino acid sequence alignment was performed using CLUSTALX . For the phylogenetic analysis, the amino acid sequences of Drosophila melanogaster [31, 32], H. armigera , M. sexta , Ostrinia furnacalis [35, 36], Chilo suppressalis , B. mori [19, 38], S. litura , and Heliothis virescens  were used. Phylogenetic analyses were conducted using the maximum likelihood method of MEGA 6.0, which was based on the Jones-Taylor-Thornton substitution model, partial deletion gaps with 95% site coverage cutoff, and nearest neighbor interchanges heuristic search . Node support of phylogenetic tree was assessed using the bootstrap method with 500 bootstrap replicates.
Profiling analysis of gene expression based on antennal transcriptome
The gene expression level was calculated using the RPKM method based on the results of antennal transcriptome analysis , in which the number of mapped reads per million reads for a gene is divided by the length of that gene. When there were multiple transcript variants for a gene, the longest one was used in calculating its expression level. The identification of differentially expressed genes between male and female antennae was conducted by using DESeq (version 1.18.0) . P-values (p-value < 0.05) and fold-changes (|fold change| > 2) were calculated to determine the differentially expressed genes in the present study.
RT-qPCR of olfactory gene expression
Single-stranded cDNAs were synthesized from 1 μg of total RNA with ReverTra Ace qPCR RT Kit (Toyobo, Kita-ku, Osaka, Japan) following the manufacturer’s recommendations. qRT-PCR was performed with SsoFast™ EvaGreen® Supermix (Bio-Rad, Hercules, California, USA), following the manufacturer’s protocols, in a CFX-96™ PCR Detection System (Bio-Rad, Hercules, California, USA). The PCR primers used are listed in Additional file 1: Table S2. Actin and AK were used as reference genes. The difference in gene expression was measured by the 2-∆∆Ct algorithm . For every gene, expressions among all studied tissues were measured with 1-day-old female antennae as control (IR1 with 1-day-old male antennae as control). Then, the logarithm of all data to base 10 was used in generating a heat map. RNA extraction was repeated three times for each sample and two or more RT-qPCR replicates were prepared for each sample.
Data were analyzed using SPSS 17.0 (SPSS Inc., Chicago, IL, USA). Significance between two samples was determined by using an independent-samples t-test. The critical P value for each test was set at 0.05. A heat map was generated by using the software PermutMatrixEN .
Sequencing and unigene assembly
Using the Illumina NextSeq 500 sequencing system, a total of 75,705,108 (11,355,766,200 bp) and 88,836,476 (12,993,896,231 bp) raw reads were obtained from the female and male samples, respectively. After removing low-quality (<Q20, adaptor and contaminating sequence reads, female and male antennae yielded 69,495,894 (10,045,828,419 bp) and 75,444,012 (11,316,601,660 bp) clean reads, respectively. Pooled female and male data were assembled into 160,299 contigs (N50 = 676 bp), with an average length of 537 nt, and 41,056 unigenes (N50 = 1301) with a mean length of 763 bp (Additional file 2: Figure S1). For the evolutionary genealogy of genes, eggNOG annotation was used to assess the transcriptome (Fig. 1), and 5105 (10.75%) unigenes were assigned to signal transduction mechanisms. Among these unigenes, male antennae showed 28 genes with a higher expression, whereas 118 genes were upregulated in the female antennae (Fig. 2).
Analysis of olfactory and gustatory genes
Approximately 43 putative OR genes, 13 putative GR genes, 16 putative IR genes, 37 putative OBP genes, 14 CSP genes, and 3 SNMP genes were identified in M. separata (Tables 1, 2, 3, 4, 5, and 6).
The OR genes included one ORco, one OR18 ortholog, 6 PR genes, and 35 general OR genes (Table 1). Eleven ORs were likely full-length genes that encoded proteins of more than 390 amino acids, including PR2, PR3, ORco, OR1, OR3, OR4, OR6-9, and OR12. Four general OR genes (OR1, 5, 33, and 34) were not effectively clustered with other moths ORs (bootstrap values < 50) in the phylogenetic analysis (Fig. 3).
GR genes included 3 CO2 receptors (GR4-6), 5 sweet receptors (GR1, 3, 10, 11, and 13), and 5 bitter taste receptor genes (Table 2). The sequence sizes of GR1 were longer than 440 amino acids, indicating that it is a nearly full-length gene. Two bitter taste receptors (GR2 and GR8) were clustered with the GR genes of B. mori, with bootstrap values > 50 in the phylogenetic analysis (Fig. 4).
Fourteen conserved IR genes were identified in M. separata, including IR25a, IR40a, IR41a, IR64a, IR68a, IR76b, IR87a, IR8a, IR93a, IR75d, two IR75p genes (IR75p1 and IR75p2), and two IR75q genes (IR75q.1 and IR75q.2). Five IRs were nearly full-length genes that encoded proteins of more than 560 amino acids, including IR75p.2, IR75d, IR75p.1, IR8a, and IR25a. In addition, two diversified IR genes (IR1 and IR2) were identified in M. separata (Table 3), and these were not clustered with the IR genes used with bootstrap values > 50 (Fig. 5).
OBPs included 2 GOBPs (GOBP1-2), ABPX, 5 PBPs (PBP1-5), and 29 OBPs. Except for PBP3, OBP9, OBP4, PBP5, OBP13, and OBP29, other 31 OBPs were likely or nearly full-length genes that encoded proteins of more than 140 amino acids. Around 29 OBPs was subdivided into 16 classic OBP genes, 10 minus-C OBP genes (OBP1, 6–9, 11, 13, 17, 27, and 29), which lack two cysteine residues (C2 and C5), and 3 plus-C OBP genes (OBP12, 14, and 24), which have more than 6 cysteine residues (Table 4, Additional file 2: Figure S2). OBP2, 3, 7, 13, 17, 27, and 29 were not clustered with the moth OBP gene with bootstrap values > 50 in the phylogenetic analysis (Fig. 6).
For the 14 CSP genes, 11 CSPs were with full ORF and 4 conserved cysteines were detected in these genes (Table 5, Additional file 2: Figure S3). Except 4 CSPs (CSP3, 6, 12, and 14), all others were clustered with orthologs of other moths in the phylogenetic analysis (bootstrap values > 80) (Fig. 7). For the SNMP genes, 2 genes harbored full ORFs (Table 6). SNMP1 and SNMP2 were clustered with orthologous genes, with bootstrap values > 50, but SNMP3 was not clustered with the other SNMPs (Fig. 8).
RNA expression of olfactory genes in the antennae
In receptor genes, except for IR75d (509.1), the RPKM value of the other receptor genes were < 40 (Fig. 9). The RPKM values of 2 PRs (PR2 and PR3) and 8 ORs (OR3, OR4, OR6, OR10, OR11, OR12, OR15, and OR21) were larger than that in ORco (5.5), and the RPKM values of IR41a, IR76b, IR75p2, and IR2 were larger than that of IR8a (1.9, co-receptor gene of IRs) but smaller than that of IR25a (13.2, co-receptor gene of IRs). The RPKM value of PR2 (38.8) was larger than that of other ORs, and the RPKM value of PR3 (19.9) was the same as those of OR12 and OR10. The RPKM values of 2 PRs (PR5 and PR6), 15ORs (OR1, OR18, OR19, OR2, OR20, OR22, OR23, OR27, OR29, OR31, OR33, OR34, OR35, OR36, and OR5), and 8 IRs (IR40a, IR64a, IR68a, IR75p.1, IR75q.1, IR75q.2, IR87a, and IR93a) were < 1 (Fig. 9). Except for GR1 (2.3), the RPKM values of the gustatory receptor genes were < 0.6.
In contrast to the receptor genes, 16 OBPs and 8 CSPs had RPKM values > 40 (Fig. 9). The RPKM value of PBP2, 3 OBPs (GOBP2, ABPX, and OBP16), and CSP4 were > 10,000, and 2 PBPs (PBP1 and PBP3), 4 OBPs (GOBP1, OBP14, OBP18, and OBP21), and 3 CSPs (CSP1, CSP3, and CSP7) were > 1000. The RPKM values of 2 CSPs (CSP12 and CSP2) and 14 OBPs (OBP1, OBP13, OBP15, OBP2, OBP23, OBP24, OBP25, OBP29, OBP3, OBP4, OBP5, OBP6, OBP7, and OBP9) were < 4. The RPKM value of SNMP2 was > 1500 (Fig. 9), but the RPKM value of SNMP3 was very small (0.08).
Expression of olfactory genes between male and female antennae
The expression levels of 43 olfactory genes, including 10 CSPs, 6 GRs, 2 IRs, 16 OBPs and 9 ORs in male antennae were significantly higher than that in female antennae, and the difference among 3 OBPs (OBP13, OBP15, and OBP3) and IR1 was > 10-fold (Fig. 10). The expression level of 38 olfactory genes, including CSP9, 4 GRs, 8 IRs, 9 OBPs and 16 ORs in female antennae was significantly higher than that in male antennae and the difference in OBP9 was higher by 10-fold. The expression level of 4 IRs, 4 OBPs and 6 ORs was the same between male and female antennae. The difference in expression levels of 16 olfactory genes, including 3 GRs, IR2, 3 CSPs, 3 OBPs, and 6 ORs between male and female antennae was influenced by their pre-mating status (Fig. 10).
The expression levels of 4 PBPs (except for PBP3) and 5 PRs (except for PR6) were significantly higher in male antennae than that in female antennae, and the difference between PBP2 and the 3 PRs (PR1, PR2, and PR4) was 10-fold higher. The expression level of PBP3 in 5-day-old females was 10-fold higher than that in males. The expression level of PR6 in 1-day-old males was significantly higher than that in females, but that in 5-day-old males was significantly smaller than in the female counterpart (Fig. 10).
Expression of olfactory genes in different pre-mating status
The expression levels of 26 olfactory genes, including 3 CSPs, 2 GRs (GR1), 3 IRs, 9 OBPs (ABPx, OBP18), and 9 ORs (OR12, OR3, and OR21) in 5-day-old moths were significantly higher than that observed in 1-day-old moths, and the difference of OBP15 was 10-fold higher (Fig. 10). The expression levels of 43 olfactory genes, including 6 CSPs (CSP3, CSP4, and CSP7), 9 IRs (IR25a and IR2), 7 GRs, 8 OBPs (OBP14), and 13 ORs (OR6) in 5-day-old moths were significantly smaller than the 1-day-old moths, and difference among CSP9, OBP9, and OR22 was 10-fold higher. The expression levels of 14 olfactory genes, including GR12, 8 OBPs (GOBP1, OBP16, and OBP21), and 5 ORs (OR10, OR4, and ORco) in 1-day-old moths were similar to those of 5-day-old moths. The differential expression of the 28 olfactory genes, including 5 CSPs (CSP1), 3 GRs, 3 IRs (IR75d), 7 OBPs (GOBP2), and 10 ORs between different pre-mating statuses was influenced by gender, particularly for CSP11, CSP12, and OBP3 (Fig. 10).
Three PBPs (PBP2, PBP3, and PBP4) and four PRs (except for PR4 and PR5) were upregulated in 5-day-old moths relative to that in 1-day-old moths, and a 10-fold difference in expression between PBP2 and PR3 was observed (Fig. 10). The expression level of PR4 in 1-day-old females was significantly higher than that in 5-day-old females. The expression level of PBP5 was the same between 1-day-old and 5-day-old moths. The difference in the expression level of PBP1 and PR5 between different pre-mating statuses was influenced by gender.
EAG responses of antennae to sex pheromones of sympatric species
Antennal EAG response to Z11-16:Ald (1.84 mV) was the strongest in the 18 sex pheromone chemicals, whereas the EAG responses to 11 sex pheromone chemicals of sympatric species (except for Z11-16:Ald, Z9-16:Ald, Z13-18:OAc, Z11-16:OH, E\Z12-14:OAc, Z9E11-14:OAc and (Z)-11-tetradecenyl acetate) were relatively very small (<0.4 mV) (Fig. 11).
In 6 of 7 sex pheromones (except Z9E11-14:OAc) with evident EAG responses (>0.4 mV), the male responses were significantly greater than that in females, and the responses to the 10−2 concentration were significantly larger than that observed with the 10−4 concentration (Fig. 11). The EAG responses of 5-day-old males to Z11-16:Ald, Z9-16:Ald, and Z13-18:OAc were significantly greater than 1-day-old males, but the EAG responses of 5-day-old males to Z11-16:OH and E\Z12-14:OAc were significantly smaller than 1-day-old males. Female EAG responses to Z9E11-14:OAc were significantly greater than males, and the EAG responses of 5-day-old females were significantly larger than 1-day-old females (Fig. 11).
EAG responses of antennae to plant odorants
The EAG responses to heptanal, Z6-nonenal, and benzaldehyde were larger than 1.2 mV, but the EAG responses to 21 plant odorants (except for (Z)-3-Hexen-1-ol, (E)-2-Hexenal, E2-5:Ald, heptanal, nonanal, Z6-nonenal, nonanol, Z6-nonenol, benzaldehyde, (1)-linalool, isoamyl acetate, ethyl hexanoate, ethyl heptanoate, ethyl isovalerate, methyl salicylate, and phenylacetic acid ethyl ester) were relatively very small (<0.4 mV) (Fig. 12).
The EAG responses at 10−2 concentration were significantly greater than that using 10−4. Furthermore, the EAG responses of males were significantly greater than that in females, and 1-day-old males were significantly greater than that in 5-day-old males. At the same time, the EAG responses of 5-day-old females to E2-5:Ald and heptanal were significantly larger than that in 1-day-old females at 10−2 concentration (Fig. 12).
Olfactory genes of M. separata
We identified a total of 126 putative olfactory genes in the antennae of M. separata, including 43 OR genes, 13 GR genes, 16 IR genes, 37 OBPs, 14 CSPs, and 3 SNMPs. The number of ORs, IRs, OBPs, and CSPs in M. separate is smaller than that in B. mori [19, 45], which was conducted using the whole genome, but was the same as that of other moths that were studied using the same protocol (antennal transcriptome), which included Agrotis ipsilon , C. suppressalis , H. armigera , O. furnacalis [35, 36]. However, the number of GRs and SNMPs was larger than those of most moths studied.
The maxillary palps, which harbor considerably fewer sensilla than antennae, are believed to specialize in taste reception in various moth species . However, 14 gustatory receptor genes were identified in M. separate antennae, indicating that the antennae might be an important taste organ for M. separata. As shown in a former study, sensilla chaetica were distributed around each antennal segment of M. separata , which was with a terminal pore and with a shape that was higher and wider than that of other sensilla and have been suggested to have contact/chemoreceptor functions . Sugar as the supplementary nutrient is essential for egg-ripening  and as the energy supplement influences the duration of migration . Five sweet receptors were identified in M. separata that specifically recognize glucose, galactose, fructose, mannose, sucrose, maltose, trehalose, raffinose, glycerol, and mannitol . GR1, a sweet receptor with higher RPKM than other GRs might be used to recognize the most important sugar for M. separata. In addition, the expressions of GR1 in 1-day-old males were larger than that in females and 5-day-old males. Two conserved members of SNMP (SNMP1 and SNMP2) were detected in various moth species . A new SNMP member SNMP3 that has been identified in S. exigua recently and other four moths  was also reported in M. separata. SNMP3 shares only less than 30% identity with SNMP1 and 2 .
Different insects have evolved various feeding behaviors in host plants and the evolution of olfactory acuity and discriminatory power of insect must be consistent with the difference in plant odor . Approximately 5 ORs, 2 GRs, 2 IRs, 7 OBPs, and 4 CSPs of M. separata were not effectively clustered with those of Lepidoptera (bootstrap values <50). These olfactory genes might be correlated with the special habitat of M. separata, including attacking rice, maize, sorghum, wheat, and other gramineous plants. M. separata larvae have evolved a specific detoxification system for benzoxazinoid, which is part of the chemical defense system of graminaceous plants .
Co-receptor for odorant detection in M. separata
IRs act in combinations of up to three subunits, comprising individual odor-specific receptors and one or two broadly expressed co-receptors, and heteromeric IR complex formation is necessary and sufficient for mediating odor-evoked electrophysiological responses in vivo and in vitro . IR25a and IR8a function as co-receptors because of their higher sequence identity compared to the other IRs , expression in all or most OSN cell bodies, and co-expression with divergent IRs in former studies [56–59]. However, in the present study, the expression levels of IR8a (RPKM = 1.9) and IR25a (RPKM = 13.2) were evidently smaller than that of IR75d (RPKM = 509.1). IR75d has been identified in various insects such as Grapholita molesta , Athetis dissimilis , Rhynchophorus ferrugineus , B. mori, Acyrthosiphon pisum, Pediculus humanus humanus, and dipterans . In Heliconius melpomene, the IR75d gene is highly expressed across various tissues and sexes . In D. melanogaster, IR75d is expressed in three different OSNs that are housed in three coeloconic sensilla (ac1, ac2, and ac4)  and is the only ORN class that occurs in more than one sensilla subtype . In addition, IR75d is co-expressed with IR84a and IR76a in ac4 sensilla . These findings suggest that IR75d might be a co-expressed receptor in M. separate.
Orco is the common olfactory receptor co-receptor in insects and is present in apparently all OR-expressing ORNs . The expression level of Orco is always higher than that of other ORs . It has been proposed that Orco is a membrane localization protein which stabilizes ORs in the dendritic membranes, as well as a chaperon molecule that facilitates in the correct protein folding of ORs and in forming heteromers with the particular OR as a cation channel [68, 69]. Downregulation of the Orco transcript not only results in weaker EAG responses [58, 70–72], but also diminishes the significant behavior preference [72–74] for host volatiles and pheromones. However, the RPKM value of ORco (5.5) in M. separata was less than 9 general ORs and 2 PRs. The low expression level of ORco in M. separata provides allows odorant receptors to mediate odorant responses also in the absence of Orco. When ORs from different species were expressed without Orco in heterologous expression assays such as Drosophila Schneider 2 (S2) and S. frugiperda 9 (SF9) cells, these also evoked ligand-specific responses . In addition, the functional expression of BmOR-1 and BmOR-3 in modified HEK 293 cells was possible without co-receptor BmOR-2 and BmOR-3 being heterologously expressed. Furthermore, modified HEK 293 cells are highly sensitive and selective responsiveness to bombykal . The Orco agonist VUAA1 did not collectively alter mating efficiency, the number of eggs produced, and the number of hatched nymphs in the bed bug, Cimex lectularius L . In fact, Orco is not essential for odorant detection when ORs are successfully inserted into the plasma membrane .
Pheromone recognition of M. separata
The sex pheromone of M. separata isolated from Japan and Taiwan was identified as a blend of (Z)-11-hexadecenyl acetate (Z11-16:Ac) and (Z)-11-hexadecenol (Z11-16:OH) [78, 79]. However, Z11-16Ald is the main female sex pheromone in mainland China [10, 11], and male moths had the highest peak response to Z11-16:Ald . In the present study, the EAG response of male oriental armyworm to Z11-16:Ald was significantly greater than that to other pheromones. In addition, the male EAG response to Z11-16:Ald was significantly greater than that in females, and the EAG response to the 10−2 concentration was significantly greater than that to the 10−4 concentration. However, the EAG responses to Z11-16:Ac was very small (<0.2 mV) and were not significantly influenced by concentration, sex, and pre-mating status. Z11-16:Ald is the main female sex pheromone in M. separata. The calling behavior of females is the most frequent and the sex pheromone titer in the sex gland is the highest in 4–5 day old M. separata . In the present study, the EAG responses to Z11-16:Ald of 5-day-old males were significantly greater than that of 1-day-old males, which was consistent with the female calling period of M. separata. In all candidate genes for sex pheromone recognition, PR2 and PBP2 were the most abundant PR and PBP, respectively, and were specifically expressed in male antennae. In addition, the expression level of PR2 and PBP2 in 5-day-old males was 10-fold higher than those found in 1-day-old males. The consistency between the EAG response to Z11-16:Ald and the expression profile of PR2 and PBP2 suggests that these two genes might be used in recognizing Z11-16:Ald, the main female sex pheromone.
In addition to Z11-16:Ald, there were also EAG responses to other sex pheromones such as Z9-16:Ald, Z11-16:OH, and E\Z12-14:OAc. Male EAG responses to these pheromones were significantly greater than female EAG responses, and the responses to the 10−2 concentration were significantly greater than that to the 10−4 concentration. These sex pheromones might be identified by male M. separata. Besides PR2 and PBP2, there were also other 5 PRs and 4 PBPs in the antennae and the expression levels of 3 PBPs and 4 PRs in male antennae were significantly higher than the females. These genes might be used in recognizing former sex pheromones. For example, MsepPR1 is used as a specific pheromone receptor for Z11-16:Ac, which is the main sex pheromone component of the Japanese M. separata . EAG responses of 5-day-old males to Z9-16:Ald was significantly greater than that of 1-day-old males in the present study, and Z9-16:Ald might be a minor sex pheromone component in M. separata female glands . In accordance with the EAG response, the expression levels of 2 PBPs (PBP1 and PBP3) and PR3 were higher in 5-day-old male moths than in the 1-day-old moths. E\Z12-14:OAc is the main female sex pheromone of O. furnacalis, which is a major pest of maize . M. separata males might be capable of distinguishing conspecific females from other species by recognizing their sex pheromone components or might search host plants with female sex pheromone of moths that infest the same host as that of M. separata. For example, M. separata might find maize plants by recognizing E\Z12-14:OAc.
Host recognition in M. separata
At least four distinct large-scale and long-distance migration events of M. separata between overwintering sites in southern China and northern temperate zones occur annually, including two northward displacements in the spring and early summer and two southward displacements in summer and fall . The migration takes place in 1–2 days after eclosion for the orient armyworm . Adults need a sufficient amount of sugar as the energy supplement for flight, and sugar is a crucial factor that influences the duration of migration . In addition, adults also need a sufficient amount of sugar for survival; otherwise, these soon perish . Under natural conditions, the adult armyworm feeds on flower nectar or honey-dew of aphids . In the present study, EAG responses of 1-day-old moths to 16 plant volatiles were significantly larger than the responses of 5-day-old moths, and expression level of the 43 olfactory genes (6 CSPs, 9 IRs, 7 GRs, 8 OBPs, and 13 ORs) in 1-day-old moths were significantly higher than those in the 5-day-old moths. These plant volatiles might be used by an adult armyworm to search for food, and might be mediated by these olfactory genes, which were differentially expressed in 1- and 5-day-old moth antennae. In particular, the EAG responses to heptanal, Z6-nonenal, and benzaldehyde were greater than 1.2 mV, which is the EAG response to sex pheromones. Heptanal, Z6-nonenal, and benzaldehyde are flower volatile compounds [82, 83]. These three flower volatiles might be very important to adult armyworms as they search for food. The RPKM values of the 7 OBPs and 4 CSPs were > 1000, whereas that of the 6 ORs and 2 IRs were > 10. Therefore, these olfactory genes might play a critical role in the recognition of host volatiles.
Spermatids are observed at the time of male M. separate emergence , whereas oocytes ripen only after the female moth has consumed a sufficient amount of food . In accordance with the larger sugar requirement of the female, the expression levels of 38 olfactory genes (CSP9, 4 GRs, 8 IRs, 9 OBPs, and 16 ORs) in females were significantly higher than those in the males (Fig. 10). However, no plant volatiles with higher female EAG responses than that in the males were observed in the present study. The inconsistency between gene expression and EAG response might be attributable to the loss of important chemicals. In fact, identifying all stimulating chemicals may be a relatively difficult task. A panel of 480 odorants, including esters, acids, aldehydes, ketones, alcohols, pyrazines, aromatics, terpenes, and sulfur compounds, against 21 larval odorant receptors (more than 10,000 odorant–receptor combinations) were conducted in fly, but no strong odorants were identified for Or2a and Or49a . Other than sex pheromones, host plant volatiles might also be used by sexually mature males to find females for mating. In Lepidoptera, several cases of synergism between plant semiochemicals and pheromones have been observed both in the laboratory and field . In our study, the male responses to 16 plant volatiles with high EAG responses (>0.4 mV) were significantly greater than the female responses, and the expression levels of 43 olfactory genes (10 CSPs, 6 GRs, 2 IRs, 16 OBPs, and 9 ORs) in males were significantly greater than in the females. These results prove that some plant volatiles are very important to male moths, and using of these volatiles might be a strategy to optimize mating opportunities. However, our understanding of how insect responses to plant stimuli aid in mating are limited .
In the present study, the EAG response to volatile organic compounds of larvae host plants was very small, which include β-myrcene, (Z)-3-hexenyl acetate, nonanal, linalool, β-caryophyllene, (E)-2-hexenal, and (Z)-3-hexen-1-ol from wheat [88, 89]. These low EAG responses may be attributable to the pre-mating status of the moths used in the present study. Only 1- and 5-day-old unmated male and female moths were used in our study. The most important behavior of 1- and 5-day-old moths involved finding their nutrition for sex maturation as well as searching for a mating partner. Therefore, volatiles of larvae host plant might not be required by 1- and 5-day-old unmated moths.
Antennal binding protein X
Volutionary genealogy of genes: Non-supervised Orthologous Groups
General odorant binding proteins
Odorant binding proteins
Olfactory receptor co-receptor
Pheromone binding proteins
Reads per kb per million reads
Quantitative reverse transcription PCR
Sensory neuron membrane proteins
Sharma HC, Davies JC. The oriental armyworm, mythimna separata (Wlk.) distribution, biology and control: a literature review. Miscellaneous report No 59. Wrights Lane: Overseas Development Administration; 1983.
Jiang X, Luo L, Zhang L, Sappington TW, Hu Y. Regulation of migration in mythimna separata (walker) in china: a review integrating environmental, physiological, hormonal, genetic, and molecular factors. Environ Entomol. 2011;40(3):516–33.
Zhao XC, Feng HQ, Wu B, Wu XF, Liu ZF, Wu KM, McNeil JN. Does the onset of sexual maturation terminate the expression of migratory behaviour in moths? a study of the oriental armyworm, mythimna separata. J Insect Physiol. 2009;55(11):1039–43.
Quo F, Wu TN, Tsai HL, Liu CL. Studies on the reporduction of the armyworm leucania separata walker (Lepidoptera: noctuidae) I. The biological characteristics of adults. Acta Ecol Sin. 1963;12(5–6):565–77.
Sakurai T, Mitsuno H, Haupt SS, Uchino K, Yokohari F, Nishioka T, Kobayashi I, Sezutsu H, Tamura T, Kanzaki R. A single sex pheromone receptor determines chemical response specificity of sexual behavior in the silkmoth bombyx Mori. PLoS Genet. 2011;7(6):e1002115.
Van Hiel M, Van Wielendaele P, Temmerman L, Van Soest S, Vuerinckx K, Huybrechts R, Broeck J, Simonet G. Identification and validation of housekeeping genes in brains of the desert locust schistocerca gregaria under different developmental conditions. BMC Mol Biol. 2009;10(1):56.
Chang HT, Guo MB, Wang B, Liu Y, Dong SL, Wang GR. Sensillar expression and responses of olfactory receptors reveal different peripheral coding in two helicoverpa species using the same pheromone components. Sci Rep. 2016;6:18742.
Feng B, Lin X, Zheng K, Qian K, Chang Y, Du Y. Transcriptome and expression profiling analysis link patterns of gene expression to antennal responses in spodoptera litura. BMC Genomics. 2015;16(1):269.
Mitsuno H, Sakurai T, Murai M, Yasuda T, Kugimiya S, Ozawa R, Toyohara H, Takabayashi J, Miyoshi H, Nishioka T. Identification of receptors of main sex-pheromone components of three lepidopteran species. Eur J Neurosci. 2008;28(5):893–902.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotech. 2011;29(7):644–52.
Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, Gibson TJ, Benton R. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010;6(8):e1001064.
Liu NY, Xu W, Papanicolaou A, Dong SL, Anderson A. Identification and characterization of three chemosensory receptor families in the cotton bollworm helicoverpa armigera. BMC Genomics. 2014;15(1):597.
Zhang T, Coates BS, Ge X, Bai S, He K, Wang Z. Male- and female-biased gene expression of olfactory-related genes in the antennae of Asian corn borer, ostrinia furnacalis (guenee) (Lepidoptera: crambidae). PLoS ONE. 2015;10(6):e0128550.
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(6):621–9.
Gu SH, Sun L, Yang RN, Wu KM, Guo YY, Li XC, Zhou JJ, Zhang YJ. Molecular characterization and differential expression of olfactory genes in the antennae of the black cutworm moth Agrotis ipsilon. PLoS ONE. 2014;9(8):e103420.
Zhang J, Wang B, Dong S, Cao D, Dong J, Walker WB, Liu Y, Wang G. Antennal transcriptome analysis and comparison of chemosensory gene families in two closely related noctuidae moths, helicoverpa armigera and H. Assulta. PLoS ONE. 2015;10(2):e0117054.
Liu NY, Zhang T, Ye ZF, Li F, Dong SL. Identification and characterization of candidate chemosensory gene families from spodoptera exigua developmental transcriptomes. Int J Biol Sci. 2015;11(9):1036–48.
Antony B, Soffan A, Jakše J, Abdelazim MM, Aldosari SA, Aldawood AS, Pain A. Identification of the genes involved in odorant reception and detection in the palm weevil rhynchophorus ferrugineus, an important quarantine pest, by antennal transcriptome analysis. BMC Genomics. 2016;17:69.
Li Q, Barish S, Okuwa S, Maciejewski A, Brandt AT, Reinhold D, Jones CD, Volkan PC. A functionally conserved gene regulatory network module governing olfactory neuron diversity. PLoS Genet. 2016;12(1):e1005780.
Fan J, Zhang Y, Francis F, Cheng D, Sun J, Chen J. Orco mediates olfactory behaviors and winged morph differentiation induced by alarm pheromone in the grain aphid, sitobion avenae. Insect Biochem Mol Biol. 2015;64:16–24.
DeGennaro M, McBride CS, Seeholzer L, Nakagawa T, Dennis EJ, Goldman C, Jasinskiene N, James AA, Vosshall LB. Orco mutant mosquitoes lose strong preference for humans and are not repelled by volatile DEET. Nature. 2013;498(7455):487–91.
Hansen IA, Rodriguez SD, Drake LL, Price DP, Blakely BN, Hammond JI, Tsujimoto H, Monroy EY, Maio WA, Romero A. The odorant receptor co-receptor from the bed bug, cimex lectularius L. PLoS ONE. 2014;9(11):e113692.
Kou R, Chow Y, Ho HY. Chemical composition of sex pheromone gland extract in female oriental armyworm pseudaletia separata walker (Lepidoptera: noctuidae) in Taiwan. Bull Inst Zool Acad Sin. 1992;31:246–50.
Martin N, Moore K, Musto CJ, Linn CE. Flight tunnel response of male european corn borer moths to cross-specific mixtures of european and Asian corn borer sex pheromones: evidence supporting a critical stage in evolution of a new communication system. J Chem Ecol. 2016;42(1):51–4.
Han EN, Gatehouse AG. The effect of adult feeding on female pre-calling period and ovarian development in a migratory moth, the oriental armyworm mythimna separata (Lepidoptera: noctuidae). Bull Entomol Res. 1991;81(04):395–400.
Mathew D, Martelli C, Kelley-Swift E, Brusalis C, Gershow M, Samuel ADT, Emonet T, Carlson JR. Functional diversity among sensory receptors in a drosophila olfactory circuit. Proc Natl Acad Sci U S A. 2013;110(23):E2134–43.
We are grateful to Yingjie Du (University of California Los Angeles) for his help in the English corrections. This work was supported financially by the Special Fund for Agro-scientific Research in the Public Interest in China (Grant No. 201203036) to YD.
This work was supported by the Special Fund for Agro-scientific Research in the Public Interest in China (Grant No. 201203036) to YDU.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its supplementary information files].
BF performed the statistical analysis, and participated in the design of the study, the transcriptome analysis and wrote part of manuscript. YH participated in the design of the study and revised the early manuscript. QG performed the RT-qPCR and the electroantennogram (EAG) recording. YD proposed the project, designed the whole study, wrote and finalized the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Authors and Affiliations
Institute of Health and Environmental Ecology, Wenzhou Medical University, University Town, Wenzhou, 325035, China
Bo Feng, Qian-shuang Guo & Yongjun Du
Ningbo City College of Vocational Technology, Xuefu Road, Yinzhou High Educational Park, NingBo, 315100, ZheJiang, China
Size distribution of all 41,056 unigenes assembled from the pooled M. separata RNA extract. Figure S2. Aligned putative full ORF of OBP gene sequences of M. separata. Six conserved cysteines were highlighted by blue color. Figure S3. Aligned putative ORF of CSP gene sequences of M. separata, CSP13 and CSP14 excluded. Four conserved cysteines were highlighted by blue color. (DOC 683 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
He, Yq., Feng, B., Guo, Qs. et al. Age influences the olfactory profiles of the migratory oriental armyworm mythimna separate at the molecular level.
BMC Genomics18, 32 (2017). https://doi.org/10.1186/s12864-016-3427-2