- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
Three closely related species, Helicoverpa armigera, H. assulta, and Heliothis virescens from Lepidoptera Noctuidae, are used as a model system for exploring sexual communication and species isolation. Pheromone receptors (PRs) previously discovered in model moth species include seven in H. armigera, six in H. assulta, and six in H. virescens. PRs named OR6, OR13, and OR16 among these species were found to be functional, characterized by an in vitro Xenopus oocytes system. Using an in vivo transgenic fly system, functional assays of OR6, OR13, and OR16 clades from three closely related Noctuidae species showed that OR13 function was highly conserved, whereas OR6 and OR16 exhibited functional divergence. Similar results were produced from assays in the Xenopus oocytes system. Combined with earlier behavioral results and electrophysiological recordings, we found corresponding relationships among pheromones, PRs, and neurons at the periphery sensory system of each species. Our results provide vital information at the neuronal and molecular level, shedding insight into the sexual communication of closely related species in Lepidoptera.
Introduction
Moth sex pheromones are released by female moths to attract conspecific males, allowing for long-distance mating attraction. Reception of moth sex pheromones among closely related species is complicated by diverse releasing and receiving pheromone signals, as well as varying components, quantities, and ratios of moth sex pheromones (Berg et al., 1995; Hansson et al., 1995; Baker et al., 2004; Wang et al., 2010; Vasquez et al., 2011; Zhang and Löfstedt, 2015). Sexual communication of closely related moth species in Lepidoptera Noctuidae has been studied over a few decades as a model system for exploring sex pheromone recognition and species isolation mechanisms (Kehat and Dunkelblum, 1990; Almaas and Mustaparta, 1991; Berg et al., 1998; Baker et al., 2004; Krieger et al., 2004; Groot et al., 2006; Gould et al., 2010; Wang et al., 2010; Wu et al., 2015; Chang et al., 2016). However, there is not sufficient evidence to explain how intra- and interspecific sexual communication signals of closely related species are discriminated (Vasquez et al., 2011; Zhang et al., 2016).
Three Lepidoptera species across two genera, Helicoverpa armigera and H. assulta in Helicoverpa and Heliothis virescens in Heliothis, are phylogenetically closely related and have been thoroughly studied. H. armigera and H. assulta are sympatrically occurring species found throughout different regions of China, and Heliothis virescens is found in America and other countries (Wang et al., 2005; Cho et al., 2008). Sex pheromone blends found in females of these three species overlap in several sex pheromone components. The major component is (Z)-11-hexadecenal (Z11-16:Ald) in H. armigera and H. virescens, and (Z)-9-hexadecenal (Z9-16:Ald) in H. assulta (Vetter and Baker, 1983; Cork et al., 1992; Baker et al., 2004), occurring in different ratios with other minor components (Nesbitt et al., 1979; Cork et al., 1992; Chang et al., 2016). Five additional compounds were identified in gland extracts of H. armigera females: hexadecanal (16: Ald), hexadecanol (16: OH), (Z)-11-hexadecenol (Z11-16:OH), (Z)-7-hexadecenal (Z7-16:Ald), and (Z)-9-tetradecenal (Z9-14:Ald) (Nesbitt et al., 1979; Dunkelblum et al., 1980; Kehat and Dunkelblum, 1990). Similarly, seven compounds were identified from gland extracts of H. assulta females: 16:Ald, (Z)-9-hexadecenyl acetate (Z9-16:OAc), (Z)-11-hexadecenyl acetate (Z11-16:OAc), hexadecanyl acetate (16:OAc), (Z)-9-hexadecenol (Z9-16:OH), Z11-16:OH, and hexadecanol (16:OH) (Cork et al., 1992; Berg and Mustaparta, 1995). However, the H. virescens female glands only produce six aldehydes and alcohols rather than acetates, including tetradecanal (14:Ald), Z9-14:Ald, Z7-16:Ald, Z9-16:Ald, Z11-16:OH, 16:Ald (Tumlinson et al., 1975; Klun et al., 1980; Vetter and Baker, 1983; Ramaswamy et al., 1985; Teal et al., 1986; Groot et al., 2006, 2009, 2013).
Field tests and behavior experiments have shown that binary pheromone blends of Z11-16:Ald and Z9-16:Ald effectively attract H. armigera males (Kehat et al., 1980; Kehat and Dunkelblum, 1990). Z11-16:OH significantly reduced catches but 16: Ald acted in opposite function when mixed with the sex pheromone principal of H. armigera (Wu et al., 1997). In addition, the pheromone component Z9-14:Ald (found in H. armigera but not H. assulta) mixed with binary pheromone blends of Z11-16:Ald and Z9-16:Ald caught more H. armigera males at lower concentrations compared to H. assulta, whereas it significantly inhibited the attraction behavior of H. armigera at higher concentrations (Gothilf et al., 1978; Kehat and Dunkelblum, 1990; Zhang et al., 2012; Wu et al., 2015). In H. assulta, addition of Z9-14:Ald or Z9-16:OH to the principal pheromone blend in certain amounts significantly reduced trap catch of male H. assulta in both field and laboratory experiments (Cork et al., 1992; Park et al., 1994; Boo et al., 1995). However, when Z9-16:OAc and Z11-16:OAc were added to binary pheromone blends of Z9-16:Ald and Z11-16:Ald at a certain ratio the male H. assulta would show attractive and mating behavior (Cork et al., 1992; Park et al., 1994). In H. virescens, males use Z11-16:Ald and Z9-14:Ald as the principal pheromone blend for upwind flight behavior (Vetter and Baker, 1983; Ramaswamy et al., 1985). When 16:Ald was added to pheromone blends of Z11-16:Ald and Z9-14:Ald, close-range sexual behaviors of male moths usually increased (Vetter and Baker, 1983). However, H. virescens does not produce acetates compared to H. armigera and H. assulta (Tumlinson et al., 1975; Klun et al., 1980; Vetter and Baker, 1983; Ramaswamy et al., 1985; Teal et al., 1986; Groot et al., 2006).
In previous studies, electrophysiological responses of sex pheromone have been recorded from a single cell within trichoid sensillum of male antennae in H. armigera, H. assulta, and H. virescens, showing specific neuron responses activated by sex pheromones (Baker et al., 2004; Gould et al., 2010; Wu et al., 2015; Chang et al., 2016; Xu et al., 2016). Genes encoding pheromone receptors (PRs), expressed on the dendritic membrane of specific olfactory receptor neurons (ORNs) in trichoid sensilla of adult male antennae, are vital to the reception of conspecific sex pheromones (Baker, 2009; Wang et al., 2010; Zhang and Löfstedt, 2015). PRs have been identified and characterized by species from genomic databases, cDNA-library screenings, and the antennal transcriptome sequencing, with seven PRs in H. armigera, six in H. assulta, and six in H. virescens (Krieger et al., 2004; Liu et al., 2012; Zhang et al., 2015a). The function and localization of PRs were demonstrated by electrophysiology methods and in situ hybridization studies (Krieger et al., 2004, 2009; Grosse-Wilde et al., 2007; Baker, 2009; Wang et al., 2010, 2016; Liu et al., 2013a; Jiang et al., 2014; Chang et al., 2016; Xu et al., 2016).
To date, several strategies for deorphanizing Lepidoptera PRs have been developed both in vitro and in vivo systems (Supplementary Table S1). The most common method to study insect ORs is in vitro heterologous expression in Xenopus oocytes (de Fouchier et al., 2014; Zhang and Löfstedt, 2015; Cui et al., 2018). Another transgenic fly lines have been used to assay OR function since 2003. The earliest system for studying OR functions was the Drosophila “empty neuron” system (Dobritsa et al., 2003). The advantage of this system is that the target OR gene is expressed in the Drosophila “empty neuron,” offering an actual cellular environment and allowing heterologous OR coupling with endogenous Orco. At the same time, the odorants can be delivered in gaseous form and combined with the Drosophila OBPs, in vivo (Hallem et al., 2004; Carey et al., 2010). However, the “empty neuron” system has some limitations for testing other ORs, such as lepidopteran pheromone receptors (Syed et al., 2010). These limitations likely arise due to some essential factors, for instance, sensory neuron membrane proteins (important for pheromone-evoked neuronal activity) are lacking in the ab3A neuron (Benton et al., 2007). However, some studies have proven that the Or67dGAL4 knock-in system is better for detecting the function of moth pheromone receptors in terms of structural, biochemical, and/or biophysical features of the at1 trichoid sensilla (Syed et al., 2010; Vasquez et al., 2013; Wang et al., 2016).
In this study, we constructed a phylogenetic tree from seven identified Lepidopteran species, and revealed orthology with closely related Noctuidae PRs. According to their evolutionary relationships and functions, we selected three sets of homologous genes, OR6, OR13, and OR16, from H. armigera, H. assulta, and H. virescens, respectively, and predicted highly conserved sequences motifs. Then, we constructed nine transgenic fly lines using the Or67dGAL4 knock-in system for further functional characterization. Specifically, we compare PR functions between the Xenopus oocytes system and the Or67dGAL4 knock-in system, as well as the relationships between PRs and neurons in the peripheral nervous system. Our results summarize the correlations among pheromones, pheromone receptors, and neurons at the periphery of the sensory system from three closely related species in Lepidoptera, as well as provide information to further detect evolutionary relationships of sex pheromones.
Materials and Methods
Insect Rearing
Drosophila stocks were fed cornmeal-agar-molasses medium and maintained under a 12 h light: 12 h dark cycle at 25°C and 60% relative humidity. The medium was changed after 10 days. Three to ten-days adults were used to test.
Fly Strains
Transgenic lines were generated according to standard procedures as described below. The open reading frame encoding OR6/OR13/OR16 genes was cloned into the pVALIUM20 vector (Ni et al., 2011). Independent homozygous UAS-OR lines (with transgene insertions into chromosome II) were generated at the Tsinghua Fly Center (Beijing, China). Driver mutant allele Or67dGAL4 stock was provided by Dr. Barry J. Dickson (Kurtovic et al., 2007). The balancer w-/w-; sp/CyO; TM3/TM6B was used to cross with homozygous driver lines. The driver line in the Or67dGAL4 mutant background was then crossed with the UAS-OR balancer line to establish a final homozygous stock w+/w+; UAS-OR/UAS-OR; Or67dGAL4/ Or67dGAL4 which expressed OR6/OR13/OR16 genes in at1 sensilla neurons. Each OR6/OR13/OR16 insertion was confirmed by sequencing genomic DNA prepared from mutant lines. The final stock was used for electrophysiological experiments.
Sequence Analysis and Phylogenetic Tree Construction
The amino acid sequences of OR6, OR13, and OR16 from H. armigera, H. assulta, and H. virescens, respectively, were aligned using ClustalX software (Version 2.1, European Bioinformatics Institute). Dendrograms were labeled by FigTree software1. The transmembrane domains of PR6, PR13, and PR16 were predicted using TMHMM Server Version 2.02. The phylogenetic tree of PRs genes in different Lepidoptera species was constructed by RaxML version 8 with Jones-Taylor-Thornton amino acid substitution model (JTT) (Stamatakis, 2014). Node support was assessed using a bootstrap method based on 1000 replicates. The PR and Odorant receptor co-receptor (Orco) data set contained 38 PR and seven Orco sequences identified in Lepidoptera [eight from H. armigera (Liu et al., 2012), seven from H. assulta (Zhang et al., 2015a), seven from H. virescens (Wang et al., 2010), eight from B. mori (Nakagawa et al., 2005; Wanner et al., 2007), five from S. exigua (Liu et al., 2013a), five from S. litura (Zhang et al., 2015b), and five from S. littoralis (Montagné et al., 2012; de Fouchier et al., 2015)]. The phylogeny of the seven moth species above was constructed on the basis of cytochrome oxidase subunit I (COI) genes.
Motif-Pattern Analysis
The motif-pattern analysis of proteins was performed broadly using the MEME online server (MEME Suite Version 4.11.2)3. A total of nine PRs from H. armigera, H. assulta, and H. virescens were selected to predict the conserved motif pattern. The parameter settings were as follow: maximum number of motifs was eight, minimum motif width was six, maximum motif width was 15, and Expectation maximization (EM) improvement threshold was 10-5.
Single Sensillum Recordings
Using a transgenic in vivo system, the OR6, OR13, and OR16 genes across three Heliothis/Helicoverpa species were respectively expressed in at1 neurons of Drosophila, and the resulting UAS-OR flies were crossed with a mutant knock-in allele Or67dGAL4 driver line. Extracellular electrophysiological recordings were performed on single at1 sensilla of one to 10 day old flies. The antenna was fixed using standard procedures (de Bruyne et al., 2001; Syed et al., 2006). The reference electrode was placed in the fly eye, under a microscope (LEICA Z16 APO, Germany) at 920 × magnification. Action potentials were recorded by inserting a tungsten wire electrode in the base or in the shaft of a sensillum of the fly antenna. Signals were amplified 10× by a high impedance pre-amplifier (IDAC-4 USB System, Syntech, Kirchzarten, Germany), sent to a PC via an analog-digital converter, and analyzed off-line with AUTOSPIKE v. 3.9 software (Syntech, Kirchzarten, Germany). The filter was set with a 500 Hz low cutoff and a three kHz high cutoff. AC signals were recorded for 10 s, starting 1 s before stimulation. Responses were calculated by counting the number of action potentials 1 s after stimulation (with a delay of 200 ms to allow the odorant to travel down the airstream), and subtracting the number counted in the second before stimulation. Three dimensional bar charts were created in SigmaPlot Version 12.5 (SYSTAT, San Jose, CA, United States). Heatmaps of different PR functions activated by sex pheromone components and analog were generated by Heml 1.0 software (Deng et al., 2014).
Odor Stimulation
In total, nine sex pheromone components and analogs, Z9-14:OAc, Z9-16:OAc, Z11-16:OAc, Z9-14:Ald, Z9-16:Ald, Z11-16:Ald, Z9-14:OH, Z9-16:OH, and Z11-16:OH, were used to screen in vivo functions of all three types ORs across three Heliothis/Helicoverpa species with paraffin oil as a control. Aliquots of sex pheromone components were dissolved in paraffin oil (v/v), and 10 μL of each solution were loaded onto a 5 × 40 mm Whatman filter paper strip, which was placed inside a Pasteur pipette. Paraffin oil alone was tested as a negative control. For dose-response relationships, serial dilutions were made in increasing doses of 0.001, 0.01, 0.1, 1, 10, and 100 μg/μL and loaded on separate filter paper strips. Each preparation was held in a humidified continuous air flow delivered by the Syntech Stimulus controller (CS-55 model, Syntech) at 1.4 L/min. Stimulus pulses were added for 300 ms. During stimulation, the compensatory flow was switched off.
Results and Discussion
Phylogenetic Analysis Reveals Orthology With Closely Related Noctuidae PRs
In some Lepidoptera species (especially in the superfamily Noctuidae), the number of PRs revealed, identified, and characterized by species were four in Spodoptera exigua (OR6, 11, 13, 16), four in Spodoptera litura (OR6, 11, 13, 16), four in Spodoptera littoralis (OR6, 11, 13, 16), seven in H. armigera (OR6, 11, 13, 14, 14b, 15, 16), six PRs in H. assulta (OR6, 11, 13, 14, 14b, 16), and six in H. virescens (OR6, 11, 13, 14, 15, 16). In addition to Noctuidae species, seven PRs (OR1, 3, 4, 5, 6, 7, 9) were identified and characterized in Bombyx mori, belonging to Bombycidae (Nakagawa et al., 2005; Wanner et al., 2007; Wang et al., 2010, 2016; Liu et al., 2012; Montagné et al., 2012; Liu et al., 2013a,b; Jiang et al., 2014; de Fouchier et al., 2015; Zhang et al., 2015a,b; Chang et al., 2016) (Figure 1A and Supplementary Table S1).
FIGURE 1. The phylogeny of pheromone receptors (PRs) from seven moth species. (A) Phylogeny and number of PRs identified by species, including Helicoverpa armigera, H. assulta, Heliothis virescens, Bombyx mori, Spodoptera exigua, S. litura and S. littoralis. (B) A phylogenetic tree of PR genes in different Lepidoptera species. Six clades (I to VI) are shown in this tree representing Orco, OR13, OR11, OR14/15, OR16, and OR6 clades, respectively.
Full-length amino acid sequences of candidate PRs genes were used to construct a phylogenetic tree from seven identified lepidopteran species including B. mori, H. armigera, H. assulta, H. virescens, S. exigua, S. litura, and S. littoralis (Figure 1B). Orthologous genes of the highly conserved co-receptor Orco, were clustered together as Clade I. As expected, sequence identity among them was very high. Another five orthologous clades were shown as noctuids species in Clade II-VI, representing clades OR6, OR11, OR13, OR14/14b/15, and OR16 (Figure 1B). The amino acid sequences of PRs across various noctuids species in OR13 clade are quite conserved, showing functional conservation; the sequences of OR6 or OR16 clade are relatively less conserved, exhibiting functional differentiation (Figure 2) (Wang et al., 2010; Liu et al., 2013b; Jiang et al., 2014; de Fouchier et al., 2015).
FIGURE 2. Motif analysis of pheromone receptors (PRs) identified from three closely related Lepidoptera species, and the alignment of amino acid sequence of three set of PRs. (A) The eight motif-pattern discovered in nine PRs from Helicoverpa armigera, H. assulta, and Heliothis virescens. (B) The locations of each motif-pattern on the predicted protein sequence from N-terminal to C-terminal. Smaller numbers indicate higher conservation. (C) The alignment of amino acid sequence of clade OR6 from H. armigera, H. assulta, and H. virescens. TM1-TM7 indicates seven transmembrane domains. Harm: H. armigera; Hass: H. assulta; Hvir: H. virescens. (D) The alignment of amino acid sequence of clade OR13 from H. armigera, H. assulta, and H. virescens. (E) The alignment of amino acid sequence of clade OR16 from H. armigera, H. assulta, and H. virescens.
Three Sets of Homologous PR Genes Selected and Cloned From Closely Related Species
Evolutionarily, H. armigera, H. assulta, and H. virescens are highly related compared with other Lepidopteran species (Wang et al., 2005; Cho et al., 2008). PRs of these three species could respond to overlapping sex pheromone components (Wang et al., 2010; Liu et al., 2013b; Jiang et al., 2014). Thus, studying evolutionary relationships among PRs in H. armigera and related species will provide valuable information on reproductive isolation.
Based on previous studies of PRs across Heliothis/Helicoverpa species, several pheromone components were used to determine response profiles of all PRs across Heliothis/Helicoverpa species, mainly using an in vitro two-electrode voltage-clamp system (Wang et al., 2010, 2016; Liu et al., 2013b; Jiang et al., 2014; Chang et al., 2016; Xu et al., 2016). We found that none of the OR11 and OR15 PRs across three Heliothis/Helicoverpa species were activated by any pheromone component tested. However, only HvirOR14 of all OR14 PRs across three Heliothis/Helicoverpa species showed response, and was activated by Z11-16:OAc and Z9-14:Ald. Similarly, OR14b from H. virescens was not identified (Supplementary Table S2). Therefore, we selected homologous OR6, OR13, and OR16, which play an important role in mating, for comparing the functions across three species. Three sets of homologous PR genes (total of nine PRs) were cloned from cDNA sequences according to the genomic database and antennal transcriptome sequence (Krieger et al., 2004; Liu et al., 2012; Zhang et al., 2015a). Subsequently, all genes were subcloned into the expression vector of the transgenic fly for further functional screening.
Sequence Analysis of Noctuidae PRs Genes
According to amino acid sequences of orthologous PR genes in the closely related species, H. armigera, H. assulta, and H. virescens, three multiple sequence alignments (OR6, OR13, and OR16) revealed relatively conserved characteristics among orthologous PRs. Each alignment contained seven transmembrane domains (Figures 2C–E), with sequence identities of 89.95, 95.54, and 94.08% corresponding to OR6, OR13, and OR16 alignments, respectively.
Nine PR sequences were used to predict highly conserved motifs. A total of eight motifs composed the most common pattern of sequence “7-6-8-5-2-3-1-4,” which represented traits with three types of ORs in H. armigera, H. assulta, and H. virescens, respectively (Figure 2B). The most typical conserved sequence patterns were located in the conserved C-terminal region as (A/G)-V-Y-(G/L/S)-(V/L)-P-W-(E/D)-(C/Y)-M-D-(T/V)-K-N-R in motif 1, F-H-Q-(A/Y/T)-S-G-C-(L/I)-L-L-L-(E/G)-C-S-Q in motif 2, Q-Q-L-I-Q-(L/I)-S-V-I-F-E-L-(V/L)-G-(S/T) in motif 3, and G-V-(T/Q)-(T/S)-M-(A/T)-(A/S)-I-L-K-T-S-(M/E)-S-Y in motif 4 (Figure 2A). The functions of these motifs were thought to be important in protein-protein interactions (Miller and Tu, 2008), especially in the formation of the OR/Orco heteromeric complex (Benton et al., 2006; Vasquez et al., 2013). In addition, another four motifs, motif 5 (H/N)-(W/C/V)-(I/F/V)-(I/L)-S-Y-(L/T)-C-(S/T/A)-(T/S/C)-(W/Y)-F-C-(M/Y)-(F/Y), motif 6 L-F-N-(L/M/I)-(I/T)-P-(M/F)-Y-(S/N)-(N/C)-(Y/L)-(A/S) -(A/R)-G-(R/M/K), motif 7 K-(I/T)-H-L-F-(Y/H)-(Y/H)-(K/R)-(D/H/N/E)-(R/K)-S-(K/E/D)-(Y/H/Q/A)-A-(M/Y), and motif 8 N-(S/A/T/R)-T-(F/Y)-(E/D)-H-(S/A)-(L/V/M)-(Y/F)-Y-(S/L/P)-(Y/V)-P-F-(D/N), had lower conservation and exhibited more sequence variation. It is possible that some amino acid residues were highly variable, resulting in functional differentiation. However, the reason for evolutionary differences of PRs presented in closely related species remains unclear.
In vivo Functional Assays of Closely Related Noctuidae PRs
In H. armigera and H. assulta, the OR6-expressing neurons in at1 sensilla mainly responded to the sex pheromone component analogs Z9-14:OH and Z9-16:OH, at a dose of 1 mg loaded in the stimulus cartridge, whereas HvirOR6-expressing neurons responded to Z9-14:Ald and analog Z9-14:OH (Figures 3A,D and Supplementary Figure S1). In a dose–response experiment, neurons in at1 sensilla started firing at doses as low as 10 ng, with Z9-14:OH and Z9-16:OH EC50 values of 3.85 × 10-5 and 5.84 × 10-5 g in H. armigera, 9.66 × 10-5 and 5.99 × 10-5 g in H. assulta, and Z9-14:Ald and Z9-14:OH EC50 values of 2.75 × 10-5 and 1.26 × 10-4 g in H. virescens (Figure 4A and Supplementary Figure S1).
FIGURE 3. Odor coding of pheromone receptors (PRs) from three closely related species. (A) Responses of OR6-expressing neurons in at1 sensilla of transgenic flies. (B) Responses of OR13-expressing neurons in at1 sensilla of transgenic flies. (C) Responses of OR16-expressing neurons in at1 sensilla of transgenic flies. (D) Heatmap of response spectra of PR-expressing neurons in at1 sensilla of transgenic flies.
FIGURE 4. Odor coding of pheromone receptors from three closely related species across concentrations. (A) Dose-responses of OR6-expressing neurons in at1 sensilla of transgenic flies. (B) Dose-responses of OR13-expressing neurons in at1 sensilla of transgenic flies. (C) Dose-responses of OR16-expressing neurons in at1 sensilla of transgenic flies.
The function of the OR13 gene was highly conserved (Figure 3D and Supplementary Figure S2). We found that OR13-expressing neurons in at1 sensilla responded specifically to the sex pheromone component Z11-16:Ald at a dose of 1 mg across three Heliothis/Helicoverpa species (Figure 3B). Dose–response results showed neurons in at1 sensilla started to respond to Z11-16:Ald at a threshold of 10 ng, and continued to receive stimulation in a concentration gradient up to 1 mg (Figure 4B and Supplementary Figure S2). The EC50 values of Z11-16:Ald were 2.13 × 10-4, 2.42 × 10-4, and 2.16 × 10-4 g in H. armigera, H. assulta, and H. virescens, respectively.
By comparison, the OR16 gene exhibited functional divergence (Figure 3D and Supplementary Figure S3). In H. armigera, the HarmOR16-expressing neurons in at1 sensilla responded to the sex pheromone components Z9-14:Ald, Z11-16:OH, and Z11-16:OAc (Figure 3C and Supplementary Figure S3). In a dose–response experiment, neurons in at1 sensilla started firing at doses as low as 10 ng, with a Z9-14:Ald EC50 value of 1.26 × 10-3 g and a Z11-16:OH EC50 value of 7.94 × 10-5 g (Figure 4C and Supplementary Figure S3). In H. assulta, the HassOR16-expressing neurons in at1 sensilla responded to the sex pheromone components Z9-16:Ald, Z9-14:Ald, and Z9-16:OH (Figure 3C and Supplementary Figure S3). In addition, neurons in at1 sensilla showed a dose–response, with a Z9-16:Ald EC50 value of 8.56 × 10-5 g (Figure 4C and Supplementary Figure S3). In H. virescens, the sex pheromone components Z11-16:OH, Z9-16:Ald, Z11-16:OAc, and Z9-14:OH activated the HvirOR16-expressing neurons in at1 sensilla (Figure 3C and Supplementary Figure S3). The dose–response experiment showed Z11-16:OH and Z11-16:OAc EC50 values of 3.89 × 10-5 and 8.88 × 10-5 g (Figure 4C and Supplementary Figure S3), respectively.
PR Functional Comparison Test Between Xenopus Oocytes and Or67dGAL4 Knock-In Systems
According to the previous functional identifications of PRs using the Xenopus oocytes system (Wang et al., 2010; Liu et al., 2013b; Jiang et al., 2014; Chang et al., 2016), we summarized the sex pheromone response profiles of PRs across H. armigera, H. assulta, and H. virescens, and the functions of these PRs using the Or67dGAL4 knock-in system (Table 1). Through a comparative analysis of different methods on functional identification, we found that ligand-binding traits of PRs detected by the Xenopus oocytes system are essentially consistent with that of the Or67dGAL4 knock-in system. This was especially true for functionally conserved PR, OR13, where the best binding-ligand of three orthologous OR13s across Heliothis/Helicoverpa species was the sex pheromone component Z11-16:Ald, regardless of which methods we used. In general, OR/Orco expressed in the Xenopus oocytes system was more sensitive to the sex pheromone components. However, the in vivo Or67dGAL4 knock-in system has generally proven to be more accurate and specific (Wang et al., 2016).
By comparison, the function of OR6 was relatively divergent. HvirOR6 was mainly tuned to Z9-14:Ald in both in vivo and in vitro systems. However, HarmOR6/Orco and HassOR6/Orco were all tuned to Z9-14:Ald, Z9-16:Ald, Z9-16:OH, and Z9-14:OH using the Xenopus oocytes system, whereas only Z9-16:OH and Z9-14:OH activated HarmOR6/HassOR6 expressing at1 neurons (Table 1). These results may be explained by additional factors; for instance, the suitability of ligand concentrations, or whether some OR genes were able to work properly in the Or67dGAL4 knock-in system. It is pointed out that Z9-14:OH is not a sex pheromone component in any of these closely related species (Nesbitt et al., 1979; Klun et al., 1980; Cork et al., 1992), but instead activates HarmOR6/ HassOR6/ HvirOR6 expressing at1 neurons. This phenomenon requires further investigation.
The function of OR16 was highly divergent and widely tuned to more than three sex pheromone components or analogs, including Z11-16: OH, Z11-16:OAc, and Z9-14:Ald. The major ligands from HarmOR16, HassOR16, and HvirOR16 using both in vivo and in vitro methods were essentially identical (Table 1).
The Relationship Between PRs and Neurons in the Peripheral Nervous System
Three closely related species use their sensitive olfactory system to specially recognize interspecific-overlapping sex pheromone components. Using previous results from the Xenopus oocytes system (Wang et al., 2010, 2016; Liu et al., 2013a; Jiang et al., 2014; Chang et al., 2016) and our results from the Or67dGAL4 knock-in system combined with in situ hybridization and electrophysiological recordings, functional characterization between neurons and odorant receptors were predicted (Table 1).
In previous studies, electrophysiological responses of peripheral sex pheromone recognition were recorded from a single sensilla within trichoid sensillum of male antennae in H. armigera, H. assulta, and H. virescens (Baker et al., 2004; Gould et al., 2010; Wu et al., 2015; Chang et al., 2016; Xu et al., 2016). A total of three trichoid sensilla subtypes have been identified to perceive sex pheromone components, A-type, B-type (missing in H. assulta) and C-type, each housing two ORNs.
Combined with behavioral results, there may be a correlation between some electrophysiological responses and the functional identification of pheromone receptors. For instance, in H. armigera, Z9-14:Ald was previously found to effectively enhance attractions at lower concentrations, and significantly inhibit attraction behavior at higher concentrations (Gothilf et al., 1978; Kehat and Dunkelblum, 1990; Zhang et al., 2012; Wu et al., 2015), whereas Z11-16: OH was found to be a behavioral inhibitor (Wu et al., 1997). Single sensillum recordings showed an “a-spike” ORN (HarmOR6 or HarmOR14b. The predictions of expressed neurons are given for each ORN) in C-type sensillum was tuned to two behavioral agonists, Z9-14:Ald and Z9-16:Ald, while a “b-spike” ORN (HarmOR16) in C-type sensillum was tuned to three behavioral antagonists, Z9-14:Ald, Z11-16:OH, and Z11-16:OAc (Chang et al., 2016, 2017; Yang et al., 2017) (Table 1). In H. assulta, an “a-spike” ORN (HassOR6 or HassOR14b) in C-type sensillum was tuned to Z9-16:Ald and Z9-14:Ald, while a “b-spike” ORN (HassOR16) in C-type sensillum was tuned to the behavioral antagonist, Z9-14:Ald, and analogs Z9-14:OH and Z9-16:OH (Chang et al., 2016, 2017; Yang et al., 2017). In H. virescens, an “a-spike” ORN (HvirOR14) in C-type sensillum was tuned to Z11-16:OAc, while a “b-spike” ORN (HvirOR16) in C-type sensillum was tuned to Z11-16:OH (interspecific inhibitor) and Z9-14:Ald (Almaas and Mustaparta, 1991; Baker et al., 2004; Wang et al., 2010) (Table 1).
In general, electrophysiological responses showed an “a-spike” ORN (predicting OR13-expressing neuron) in A-type sensillum across all three species was activated by the sex pheromone component Z11-16:Ald, but another “b-spike” ORN (OR11) in A-type sensillum is still uncharacterized (Table 1). In A-type sensillum, the functions of expressed ORs are relatively conserved. In addition, the number of A-type sensilla confers a larger proportion of all trichoid sensilla in H. armigera than in H. assulta, in accordance with the understanding that Z11-16: Ald is major sex pheromone component in H. armigera (Chang et al., 2016).
One “a-spike” ORN (OR14b or OR6) in B-type trichoid sensillum is known to be mainly tuned to the sex pheromone component Z9-14: Ald, whereas none of ligands activate a “b-spike” ORN (OR15) in B-type sensillum.
Overall, we summarized the relationships among sensilla, neurons, and PRs involving sex pheromone recognition in the peripheral-coding olfactory system of three Heliothis/Helicoverpa species (Table 1). It is evident that neuron function in type-A trichoid sensilla completely matched the function of PRs (OR13 and OR11). However, relationships between neurons in type-B or -C trichoid sensilla and PRs did not fully clarified. The Or67dGAL4 knock-in system used to detect the function of moth pheromone receptor is nearly identical to the Xenopus oocytes system (Wang et al., 2016). A few functional differences are observed between PRs and endogenous neurons in moths which may be driven by many factors such as the cell environment, gene expression, lack of accessories, and category and concentration of ligand. In addition, the functions of OR14b or OR6 in H. armigera and H. assulta still exist differences in previous studies (Table 1) (Jiang et al., 2014; Chang et al., 2016; Yang et al., 2017). Therefore, the CRISPR/Cas9 genome editing technique combined with electrophysiological response assays are needed for functional characterization of OR14b (OR6). It is better for elucidation of the molecular and neuronal mechanisms of sex pheromone identification.
The Evolution of Lepidoptera PRs Selectivity
Three Heliothis/Helicoverpa male species can perceive respective sex pheromone components released from their female pheromone blends. A few hypotheses have been proposed on how variation is generated during pheromone evolution of closely related species, such as the “asymmetric tracking” hypothesis and the gene duplication hypothesis (Phelan, 1992; Gould et al., 2010; Heckel, 2010). However, it is still elusive how subtle variations of sex pheromone components are precisely distinguished by males of different species. Certain moth PRs of closely related species are evolutionarily conserved under strong selective pressure, whereas PRs are more functionally divergent if relaxed from evolutionary constraint (Zhang and Löfstedt, 2013, 2015). The latter is broadly tuned to the behavioral antagonists and agonist, which efficiently increased the specificity and selectivity of interspecific pheromone detection (Zhang and Löfstedt, 2015). This is consistent with our finding that OR16 from three closely related species exhibits largely functional divergences. The function of HarmOR16 from H. armigera has been confirmed to be activated by the pheromone antagonist Z11-16:OH, which regulates optimal mating time and influences fecundity (Chang et al., 2017). One study revealed that single mutations in PRs across Asian and European corn borers selectively altered pheromone recognition (Leary et al., 2012). Another study showed that two site mutations of HassOR14b changed ligand selectivity (Yang et al., 2017). Thus, the evolutionary relationship of structure and function of PRs in closely related Lepidoptera species will help reveal the mechanisms underlying reproductive isolation and speciation.
Author Contributions
BW, YL, and G-RW designed the experiments. BW performed the experiments and analyzed the data. YL and G-RW contributed reagents, materials, and gene identification. BW and G-RW wrote and revised the paper.
Funding
This work was funded by the National Natural Science Foundation of China (31725023, 31621064, 31402023, and 31230062) and China Postdoctoral Science Foundation (2014M550905).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. A. Ray for kindly providing the Or67dGAL4 mutant knock-in fly line with Dr. B. Dickson’s permission; Dr. J. Q. Ni for kindly providing the pVALIUM20 vector and for the transformation service; and Dr. J. Shen for kindly providing the balancer w-; sp/CyO; TM3/TM6B strain.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.01188/full#supplementary-material
FIGURE S1 | SSR traces from OR6-expressing neurons in at1 sensilla of Drosophila in response to pheromone compounds.
FIGURE S2 | SSR traces from OR13-expressing neurons in at1 sensilla of Drosophila in response to pheromone compounds.
FIGURE S3 | SSR traces from OR16-expressing neurons in at1 sensilla of Drosophila in response to pheromone compounds.
TABLE S1 | Functional characterizations of PR genes in some Lepidopteran species.
TABLE S2 | Functional characterizations of PRs in three Heliothis/Helicoverpa species in vitro.
Footnotes
- ^ http://tree.bio.ed.ac.uk/software/figtree/
- ^ http://www.cbs.dtu.dk/services/TMHMM/
- ^ http://meme-suite.org/tools/meme
References
Almaas, T. J., and Mustaparta, H. (1991). Heliothis virescens: response characteristics of receptor neurons in Sensilla Trichodea type 1 and type 2. J. Chem. Ecol. 17, 953–972. doi: 10.1007/BF01395602
Baker, T. C. (2009). Nearest neural neighbors: moth sex pheromone receptors HR11 and HR13. Chem. Senses 34, 465–468. doi: 10.1093/chemse/bjp025
Baker, T. C., Ochieng’, S. A., Cosse, A. A., Lee, S. G., Todd, J. L., Quero, C., et al. (2004). A comparison of responses from olfactory receptor neurons of Heliothis subflexa and Heliothis virescens to components of their sex pheromone. J. Comp. Physiol. A 190, 155–165. doi: 10.1007/s00359-003-0483-2
Benton, R., Sachse, S., Michnick, S. W., and Vosshall, L. B. (2006). Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 4:e20. doi: 10.1371/journal.pbio.0040020
Benton, R., Vannice, K. S., and Vosshall, L. B. (2007). An essential role for a CD36- related receptor in pheromone detection in Drosophila. Nature 450, 289–293. doi: 10.1038/nature06328
Berg, B. G., Almaas, T. J., Bjaalie, J. G., and Mustaparta, H. (1998). The macroglomerular complex of the antennal lobe in the tobacco budworm moth Heliothis virescens: specified subdivision in four compartments according to information about biologically significant compounds. J. Comp. Physiol. A 183, 669–682. doi: 10.1007/s003590050290
Berg, B. G., and Mustaparta, H. (1995). The significance of major pheromone components and interspecific signals as expressed by receptor neurons in the oriental tobacco budworm moth, Helicoverpa assulta. J. Comp. Physiol. A 177, 683–694. doi: 10.1007/BF00187627
Berg, B. G., Tumlinson, J. H., and Mustaparta, H. (1995). Chemical communication in Heliothine moths. IV. Receptor neuron responses to pheromone compounds and formate analogs in the male tobacco budworm moth Heliothis virescens. J. Comp. Physiol. A 177, 527–534. doi: 10.1007/BF00207182
Boo, K. S., Park, K. C., Hall, D. R., Cork, A., Berg, B. G., and Mustaparta, H. (1995). (Z)-9-tetradecenal: a potent inhibitor of pheromone-mediated communication in the oriental tobacco budworm moth, Helicoverpa assulta. J. Comp. Physiol. A 177, 695–699. doi: 10.1007/BF00187628
Carey, A. F., Wang, G., Su, C. Y., Zwiebel, L. J., and Carlson, J. R. (2010). Odorant reception in the malaria mosquito Anopheles gambiae. Nature 464, 66–71. doi: 10.1038/nature08834
Chang, H. T., Guo, M. B., Wang, B., Liu, Y., Dong, S. L., and Wang, G. R. (2016). Sensillar expression and responses of olfactory receptors reveal different peripheral coding in two Helicoverpa species using the same pheromone components. Sci. Rep. 6:18742. doi: 10.1038/srep18742
Chang, H. T., Liu, Y., Ai, D., Jiang, X. C., Dong, S. L., and Wang, G. R. (2017). A pheromone antagonist regulates optimal mating time in the moth Helicoverpa armigera. Curr. Biol 27, 1610.–1615. doi: 10.1016/j.cub.2017.04.035
Cho, S., Mitchell, A., Mitter, C., Regier, J., Matthews, M., and Robertson, R. (2008). Molecular phylogenetics of heliothine moths (Lepidoptera: Noctuidae: Heliothinae), with comments on the evolution of host range and pest status. Syst. Entomol. 33, 581–594. doi: 10.1111/j.1365-3113.2008.00427.x
Cork, A., Boo, K. S., Dunkelblum, E., Hall, D. R., Jee-Rajunga, K., Kehat, M., et al. (1992). Female sex pheromone of oriental tobacco budworm, Helicoverpa assulta (Guenee) (Lepidoptera: Noctuidae): identification and field testing. J. Chem. Ecol. 18, 403–418. doi: 10.1007/BF00994240
Cui, W. C., Wang, B., Guo, M. B., Liu, Y., Jacquin-Joly, E., Yan, S. C., et al. (2018). A receptor-neuron correlate for the detection of attractive plant volatiles in Helicoverpa assulta (Lepidoptera: Noctuidae). Insect Biochem. Mol. Biol. 97, 31–39. doi: 10.1016/j.ibmb.2018.04.006
de Bruyne, M., Foster, K., and Carlson, J. R. (2001). Odor coding in the Drosophila antenna. Neuron 30, 537–552. doi: 10.1016/S0896-6273(01)00289-6
de Fouchier, A., Montagné, N., Mirabeau, O., and Jacquin-Joly, E. (2014). “Current views on the function and evolution of olfactory receptors in Lepidoptera,” in Short Views on Insect Biochemistry and Molecular Biology, Vol. 2, eds R. Chandrasekar, B. K. Tyagi, Z. Z. Gui, and G. R. Reeck (Manhattan, KS: Internation Book Mission), 385–408.
de Fouchier, A., Sun, X., Monsempes, C., Mirabeau, O., Jacquin-Joly, E., and Montagné, N. (2015). Evolution of two receptors detecting the same pheromone compound in crop pest moths of the genus Spodoptera. Front. Ecol. Evol. 3:95. doi: 10.3389/fevo.2015.00095
Deng, W. K., Wang, Y. B., Liu, Z. X., Cheng, H., and Xue, Y. (2014). HemI: a toolkit for illustrating Heatmaps. PLoS One 9:e111988. doi: 10.1371/journal.pone.0111988
Dobritsa, A. A., van der Goes van Naters, W., Warr, C. G., Steinbrecht, R. A., and Carlson, J. R. (2003). Integrating the molecular and cellular basis of odor coding in the Drosophila antenna. Neuron 37, 827–841. doi: 10.1016/S0896-6273(03)00094-1
Dunkelblum, E., Gothilf, S., and Kehat, M. (1980). Identification of the sex pheromone of the cotton hollworm, Heliothis armigera, in Israel. Phytoparasitica 8, 209–211. doi: 10.1007/BF03158319
Gothilf, S., Kehat, M., Jacobson, M., and Galun, R. (1978). Sex attractants for male Heliothis armigera (Hbn.). Cell. Mol. Life Sci. 34, 853–854. doi: 10.1007/BF01939662
Gould, F., Estock, M., Hillier, N. K., Powell, B., Groot, A. T., Ward, C. M., et al. (2010). Sexual isolation of male moths explained by a single pheromone response QTL containing four receptor genes. Proc. Natl. Acad. Sci. U.S.A. 107, 8660–8665. doi: 10.1073/pnas.0910945107
Groot, A. T., Estock, M. L., Horovitz, J. L., Santangelo, R. G., Schal, C., and Gould, F. (2009). QTL analysis of sex pheromone blend differences between two closely related moths: insights into divergence in biosynthetic pathways. Insect Biochem. Mol. Biol. 39, 568–577. doi: 10.1016/j.ibmb.2009.05.002
Groot, A. T., Horovitz, J. L., Hamilton, J., Santangelo, R. G., Schal, C., and Gould, F. (2006). Experimental evidence for interspecific directional selection on moth pheromone communication. Proc. Natl. Acad. Sci. U.S.A. 103, 5858–5863. doi: 10.1073/pnas.0508609103
Groot, A. T., Staudacher, H., Barthel, A., Inglis, O., Schofl, G., Santangelo, R. G., et al. (2013). One quantitative trait locus for intra- and interspecific variation in a sex pheromone. Mol. Ecol. 22, 1065–1080. doi: 10.1111/mec.12171
Grosse-Wilde, E., Gohl, T., Bouche, E., Breer, H., and Krieger, J. (2007). Candidate pheromone receptors provide the basis for the response of distinct antennal neurons to pheromonal compounds. Eur. J. Neurosci. 25, 2364–2373. doi: 10.1111/j.1460-9568.2007.05512.x
Grosse-Wilde, E., Svatos, A., and Krieger, J. (2006). A pheromone-binding protein mediates the bombykol-induced activation of a pheromone receptor in vitro. Chem. Senses 31, 547–555. doi: 10.1093/chemse/bjj059
Hallem, E. A., Ho, M. G., and Carlson, J. R. (2004). The molecular basis of odor coding in the Drosophila antenna. Cell 117, 965–979. doi: 10.1016/j.cell.2004.05.012
Hansson, B. S., Almaas, T. J., and Anton, S. (1995). Chemical communication in heliothine moths. V. antennal lobe projection patterns of pheromone-detecting olfactory receptor neurons in the male Heliothis virescens (Lepidoptera, Noctuidae). J. Comp. Physiol. A 177, 535–543. doi: 10.1007/BF00207183
Heckel, D. G. (2010). Smells like a new species: gene duplication at the periphery. Proc. Natl. Acad. Sci. U.S.A. 107, 9481–9482. doi: 10.1073/pnas.1004511107
Jiang, X. J., Guo, H., Di, C., Yu, S., Zhu, L., Huang, L. Q., et al. (2014). Sequence similarity and functional comparisons of pheromone receptor orthologs in two closely related Helicoverpa species. Insect Biochem. Mol. Biol. 48, 63–74. doi: 10.1016/j.ibmb.2014.02.010
Kehat, M., and Dunkelblum, E. (1990). Behavioral responses of male Heliothis armigera (Lepidoptera: Noctuidae) moths in a flight tunnel to combinations of components identified from female sex pheromone glands. J. lnsect Behav. 3, 75–83. doi: 10.1007/BF01049196
Kehat, M., Gothilf, S., Dunkelblum, E., and Greenberg, S. (1980). Field evaluation of female sex pheromone components of the cotton bollworm, Heliothis armigera. Entomol. Exp. Appl. 27, 188–193. doi: 10.1111/j.1570-7458.1980.tb02963.x
Klun, J. A., Bierl-Leonhardt, B. A., Plimmer, J. R., Sparks, A. N., Primiani, M., Chapman, O. L., et al. (1980). Sex pheromone chemistry of the female tobacco budworm moth Heliothis virescens. J. Chem. Ecol. 6, 177–183. doi: 10.1007/BF00987536
Krieger, J., Gondesen, I., Forstner, M., Gohl, T., Dewer, Y., and Breer, H. (2009). HR11 and HR13 receptor-expressing neurons are housed together in pheromone-responsive Sensilla Trichodea of male Heliothis virescens. Chem. Senses 34, 469–477. doi: 10.1093/chemse/bjp012
Krieger, J., Grosse-Wilde, E., Gohl, T., Dewer, Y. M. E., Raming, K., and Breer, H. (2004). Genes encoding candidate pheromone receptors in a moth (Heliothis virescens). Proc. Natl. Acad. Sci. U.S.A. 101, 11845–11850. doi: 10.1073/pnas.0403052101
Kurtovic, A., Widmer, A., and Dickson, B. J. (2007). A single class of olfactory neurons mediates behavioural responses to a Drosophila sex pheromone. Nature 446, 542–546. doi: 10.1038/nature05672
Leary, G. P., Allen, J. E., Bunger, P. L., Luginbill, J. B., Linn, C. E., Macallister, I. E., et al. (2012). Single mutation to a sex pheromone receptor provides adaptive specificity between closely related moth species. Proc. Natl. Acad. Sci. U.S.A. 109, 14081–14086. doi: 10.1073/pnas.1204661109
Liu, C., Liu, Y., Walker, W. B., Dong, S. L., and Wang, G. R. (2013a). Identification and functional characterization of sex pheromone receptors in beet armyworm Spodoptera exigua (Hubner). Insect Biochem. Mol. Biol. 43, 747–754. doi: 10.1016/j.ibmb.2013.05.009
Liu, Y., Liu, C., Lin, K., and Wang, G. (2013b). Functional specificity of sex pheromone receptors in the cotton bollworm Helicoverpa armigera. PLoS One 8:e62094. doi: 10.1371/journal.pone.0062094
Liu, Y., Gu, S. H., Zhang, Y. J., Guo, Y. Y., and Wang, G. R. (2012). Candidate olfaction genes identified within the Helicoverpa armigera antennal transcriptome. PLoS One 7:e48260. doi: 10.1371/journal.pone.0048260
Miller, R., and Tu, Z. J. (2008). Odorant receptor C-terminal motifs in divergent insect species. J. Insect Sci. 8, 1–10. doi: 10.1673/031.008.5301
Montagné, N., Chertemps, T., Brigaud, I., François, A., François, M. C., de Fouchier, A., et al. (2012). Functional characterization of a sex pheromone receptor in the pest moth Spodoptera littoralis by heterologous expression in Drosophila. Eur. J. Neurosci. 36, 2588–2596. doi: 10.1111/j.1460-9568.2012.08183.x
Nakagawa, T., Pellegrino, M., Sato, K., Vosshall, L. B., and Touhara, K. (2012). Amino acid residues contributing to function of the heteromeric insect olfactory receptor complex. PLoS One 7:e32372. doi: 10.1371/journal.pone.0032372
Nakagawa, T., Sakurai, T., Nishioka, T., and Touhara, K. (2005). Insect sex-pheromone signals mediated by specific combinations of olfactory receptors. Science 307, 1638–1642. doi: 10.1126/science.1106267
Nesbitt, B. F., Beevor, P. S., Hall, D. R., and Lester, R. (1979). Female sex pheromone components of the cotton bollworm, Heliothis armigera. J. Insect Physiol. 25, 535–541. doi: 10.1016/S0022-1910(79)80012-8
Ni, J. Q., Zhou, R., Czech, B., Liu, L. P., Holderbaum, L., Yang-Zhou, D., et al. (2011). A genome-scale shRNA resource for transgenic RNAi in Drosophila. Nat. Methods 8, 405–407. doi: 10.1038/nmeth.1592
Park, K. C., Cork, A., Boo, K. S., and Hall, D. R. (1994). Biological activity of female sex pheromone of the oriental tobacco budworm, Helicoverpa assulta (Guenbe) (Lepidoptera: Noctuidae): electroantennography, wind tunnel observations and field trapping. Korean J. Appl. Entomol. 33, 26–32.
Phelan, P. L. (1992). “Evolution of sex pheromones and the role of asymmetric tracking,” in Insect Chemical Ecology: an Evolutionary Approach, eds B. D. Roitberg and M. B. Isman (New York, NY: Chapman & Hall), 265–314.
Ramaswamy, S. B., Randle, S. A., and Ma, W. K. (1985). Field evaluation of the sex-pheromone components of Heliothis virescens (Lepidoptera, Noctuidae) in cone traps. Environ. Entomol. 14, 293–296. doi: 10.1093/ee/14.3.293
Sakurai, T., Nakagawa, T., Mitsuno, H., Mori, H., Endo, Y., Tanoue, S., et al. (2004). Identification and functional characterization of a sex pheromone receptor in the silkmoth Bombyx mori. Proc. Natl. Acad. Sci. U.S.A. 101, 16653–16658. doi: 10.1073/pnas.0407596101
Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and postanalysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033
Syed, Z., Ishida, Y., Taylor, K., Kimbrell, D. A., and Leal, W. S. (2006). Pheromone reception in fruit flies expressing a moth’s odorant receptor. Proc. Natl. Acad. Sci. U.S.A. 103, 16538–16543. doi: 10.1073/pnas.0607874103
Syed, Z., Kopp, A., Kimbrell, D. A., and Leal, W. S. (2010). Bombykol receptors in the silkworm moth and the fruit fly. Proc. Natl. Acad. Sci. U.S.A. 107, 9436–9439. doi: 10.1073/pnas.1003881107
Teal, P. E. A., Tumlinson, J. H., and Heath, R. R. (1986). Chemical and behavioral analyses of volatile sex pheromone components released by calling Heliothis virescens (F.) females (Lepidoptera: Noctuidae). J. Chem. Ecol. 12, 107–126. doi: 10.1007/BF01045595
Tumlinson, J. H., Hendricks, P. E., Mitchell, E. R., Doolittle, R. E., and Brennan, M. M. (1975). Isolation, identification and synthesis of the sex pheromone of the tobacco budworm. J. Chem. Ecol. 1, 203–214. doi: 10.1007/BF00987869
Vasquez, G. M., Fischer, P., Grozinger, C. M., and Gould, F. (2011). Differential expression of odorant receptor genes involved in the sexual isolation of two Heliothis moths. Insect Mol. Biol. 20, 115–124. doi: 10.1111/j.1365-2583.2010.01044.x
Vasquez, G. M., Syed, Z., Estes, P. A., Leal, W. S., and Gould, F. (2013). Specificity of the receptor for the major sex pheromone component in Heliothis virescens. J. Insect Sci. 13:160. doi: 10.1673/031.013.16001
Vetter, R. S., and Baker, T. C. (1983). Behavioral responses of male Heliothis virescens in a sustained flight tunnel to combinations of seven compounds identified from female sex pheromone glands. J. Chem. Ecol. 9, 747–759. doi: 10.1007/BF00988780
Wang, B., Liu, Y., He, K., and Wang, G. R. (2016). Comparison of research methods for functional characterization of insect olfactory receptors. Sci. Rep. 6:32806. doi: 10.1038/srep32806
Wang, G., Vasquez, G. M., Schal, C., Zwiebel, L. J., and Gould, F. (2010). Functional characterization of pheromone receptors in the tobacco budworm Heliothis virescens. Insect Mol. Biol. 20, 125–133. doi: 10.1111/j.1365-2583.2010.01045.x
Wang, H. L., Zhao, C. H., and Wang, C. Z. (2005). Comparative study of sex pheromone composition and biosynthesis in Helicoverpa armigera, H. assulta and their hybrid. Insect Biochem. Mol. Biol. 35, 575–583. doi: 10.1016/j.ibmb.2005.01.018
Wanner, K. W., Anderson, A. R., Trowell, S. C., Theilmann, D. A., Robertson, H. M., and Newcomb, R. D. (2007). Female-biased expression of odourant receptor genes in the adult antennae of the silkworm, Bombyx mori. Insect Mol. Biol. 16, 107–119. doi: 10.1111/j.1365-2583.2007.00708.x
Wu, D. M., Yan, Y. H., and Cui, J. R. (1997). Sex pheromone components of Helicoverpa armigera: chemical analysis and field tests. Insect Sci. 4, 350–356. doi: 10.1111/j.1744-7917.1997.tb00109.x
Wu, H., Xu, M., Hou, C., Huang, L. Q., Dong, J. F., and Wang, C. Z. (2015). Specific olfactory neurons and glomeruli are associated to differences in behavioral responses to pheromone components between two Helicoverpa species. Front. Behav. Neurosci. 9:206. doi: 10.3389/fnbeh.2015.00206
Xu, M., Guo, H., Hou, C., Wu, H., Huang, L. Q., and Wang, C. Z. (2016). Olfactory perception and behavioral effects of sex pheromone gland components in Helicoverpa armigera and Helicoverpa assulta. Sci. Rep. 6:22998. doi: 10.1038/srep22998
Yang, K., Huang, L.-Q., Ning, C., and Wang, C.-Z. (2017). Two single-point mutations shift the ligand selectivity of a pheromone receptor between two closely related moth species. Elife 6:e29100. doi: 10.7554/eLife.29100
Zhang, D. D., and Löfstedt, C. (2013). Functional evolution of a multigene family: orthologous and paralogous pheromone receptor genes in the turnip moth, Agrotis segetum. PLoS One 8:e77345. doi: 10.1371/journal.pone.0077345
Zhang, D. D., and Löfstedt, C. (2015). Moth pheromone receptors: gene sequences, function, and evolution. Front. Ecol. Evol. 3:105. doi: 10.3389/fevo.2015.00105
Zhang, D. D., Wang, H. L., Schultze, A., Froß, H., Francke, W., Krieger, J., et al. (2016). Receptor for detection of a Type II sex pheromone in the winter moth Operophtera brumata. Sci. Rep. 6:18576. doi: 10.1038/srep18576
Zhang, J., Wang, B., Dong, S., Cao, D., Dong, J., Walker, W. B., et al. (2015a). Antennal transcriptome analysis and comparison of chemosensory gene families in two closely related Noctuidae moths, Helicoverpa armigera and H. assulta. PLoS One 10:e0117054. doi: 10.1371/journal.pone.0117054
Zhang, J., Yan, S. W., Liu, Y., Jacquin-Joly, E., Dong, S. L., and Wang, G. R. (2015b). Identification and functional characterization of sex pheromone receptors in the common cutworm (Spodoptera litura). Chem. Senses 40, 7–16. doi: 10.1093/chemse/bju052
Keywords: sex pheromones, Helicoverpa armigera, Helicoverpa assulta, Heliothis virescens, pheromone receptors, transgenic fly
Citation: Wang B, Liu Y and Wang G-R (2018) Proceeding From in vivo Functions of Pheromone Receptors: Peripheral-Coding Perception of Pheromones From Three Closely Related Species, Helicoverpa armigera, H. assulta, and Heliothis virescens. Front. Physiol. 9:1188. doi: 10.3389/fphys.2018.01188
Received: 31 May 2018; Accepted: 07 August 2018;
Published: 30 August 2018.
Edited by:
Peng He, Guizhou University, ChinaReviewed by:
Haonan Zhang, University of California, Riverside, United StatesFeng Liu, Michigan State University, United States
Longwa Zhang, Anhui Agricultural University, China
Copyright © 2018 Wang, Liu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yang Liu, yangliu@ippcaas.cn Gui-Rong Wang, wangguirong@caas.cn; grwang@ippcaas.cn