Original Research ARTICLE
Characterization and Comparative Analysis of Olfactory Receptor Co-Receptor Orco Orthologs Among Five Mirid Bug Species
- 1State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
- 2Department of Plant Protection, College of Horticulture and Plant Protection, Yangzhou University, Yangzhou, China
- 3DanDong Entry-Exit Inspection and Quarantine Bureau, Dandong, China
- 4Key Laboratory of Tea Quality and Safety Control, Ministry of Agriculture, Tea Research Institute, Chinese Academy of Agricultural Sciences, Hangzhou, China
The phytophagous mirid bugs of Apolygus lucorum, Lygus pratensis as well as three Adelphocoris spp., including Adelphocoris lineolatus, A. suturalis, and A. fasciaticollis are major pests of multiple agricultural crops in China, which have distinct geographical distribution and occurrence ranges. Like many insect species, these bugs heavily rely on olfactory cues to search preferred host plants, thereby investigation on functional co-evolution and divergence of olfactory genes seems to be necessary and is of great interest. In the odorant detection pathway, olfactory receptor co-receptor (Orco) plays critical role in the perception of odors. In this study, we identified the full-length cDNA sequences encoding three putative Orcos (AsutOrco, AfasOrco, and LpraOrco) in bug species of A. suturalis, A. fasciaticollis, and L. pratensis based on homology cloning method. Next, sequence alignment, membrane topology and gene structure analysis showed that these three Orco orthologs together with previously reported AlinOrco and AlucOrco shared high amino acid identities and similar topology structure, but had different gene structure especially at the length and insertion sites of introns. Furthermore, the evolutional estimation on the ratios of non-synonymous to synonymous (Ka/Ks) revealed that Orco genes were under strong purifying selection, but the degrees of variation were significant different between genera. The results of quantitative real-time PCR experiments showed that these five Orco genes had a similar antennae-biased tissue expression pattern. Taking these data together, it is thought that Orco genes in the mirid species could share conserved olfaction roles but had different evolution rates. These findings would lay a foundation to further investigate the molecular mechanisms of evolutionary interactions between mirid bugs and their host plants, which might in turn contribute to the development of pest management strategy for mirid bugs.
Due to long-term adoption of transgenic Bt (Bacillus thuringiensis) cotton and the associated reduction in broad-spectrum insecticide used for controlling Helicoverpa spp. (Wu et al., 2008), several species of the mirid bugs (Hemiptera: Miridae) including Apolygus lucorum, Lygus pratensis as well as three Adelphocoris spp., including Adelphocoris lineolatus, A. suturalis and A. fasciaticollis have become most important pest species in cotton fields of northern China (Lu et al., 2010). Besides cotton, these polyphagous mirid species cause severe destructions to many other important crops including vegetables, fruits trees and tea plants (Lu and Wu, 2008). It was reported that these five mirid species are significantly different in geographic distribution and seasonal abundance in China (Lu et al., 2008a). The A. lucorum is widely distributed in whole China, while three Adelphocoris species and L. pratensis occur mainly in Yangtze River region and the northern parts of Yellow River region, and in the colder region of northwest China, respectively (Lu and Wu, 2008). The screening of overwintering and early season host plant ranges suggested that mirid bugs from different regions employed distinct host plant ranges to survive winter and early spring, and these differences are significantly linked to their reliance on local plants (Lu et al., 2011). Consequently, the interactions between mirid species and local host plants should play crucial roles in determining ecological landscape-level especially their different geographic distribution and seasonal abundance. A better understanding of the underlying species-preferential host plants tracking would help to define co-evolution between different mirid species and their host plants, and ultimately facilitate the development of regional forecasting and pest management strategies.
Insect olfaction plays important roles in locating host plant. Several classes of molecules including odorant binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors (ORs), sensory neuron membrane proteins (SNMPs) and odorant degradation enzymes (ODEs) play important roles in odorant signal transduction pathway (Leal, 2013). ORs located in the dendrite membrane of olfactory sensory neurons (OSNs) and are considered to play a central role in identifying the distinct odorants and activating the OSNs (Clyne et al., 1999; Hallem et al., 2004). Compared with mammal ORs, insect ORs have seven transmembrane domains (TMDs) but employ a “reversed” topology with their N-terminus inside the cell and the C-terminus exposed to the external environment (Benton et al., 2006; Lundin et al., 2007; Hull et al., 2012). To detect the odorants, ORs could interact with a conserved olfactory receptor co-receptor (Orco) and then form ligand-gated ion channels (Sato et al., 2008; Wicher et al., 2008).
Orco is previously referred to as OR83b in Drosophila melanogaster, OR2 in Bombyx mori, and OR7 in mosquito species (Vosshall and Hansson, 2011). Conventional ORs demonstrate low sequence identity, whereas Orco is strikingly well conserved across insect species. It was reported that Orco has no direct relation with odor binding or discrimination (Nichols and Luetje, 2010; Nichols et al., 2011), but is essential for ion channel formation and olfactory cues transduction. In fact, Orco could interact with conventional ORs to form heterodimeric complexes, whereas conventional ORs were responsible for specifically binding to structurally diverse odorants (Larsson et al., 2004; Benton et al., 2006). Also, Orco was confirmed to be activated by VUAA1 as a functional ion channel in homomeric complex, even in the absence of conventional olfactory receptors (Jones et al., 2011). However, VU0183254, one of the analogs of VUAA1, showed the ability to “lock” hemomeric and homomeric ion channels in a non-competitive way due to its affinity to Orco (Jones et al., 2012). Coincidentally, these functional hemomeric and homomeric channels can be also blocked by amiloride derivatives when they were activated (Pask et al., 2013). Disruption in the transcript expression of Orco could significantly impair olfactory behavior responses in all the tested insect species, including D. melanogaster (Larsson et al., 2004), Acyrthosiphon pisum (Zhang et al., 2017), Locust amigratoria (Li et al., 2016), Spodoptera litura (Dong et al., 2013), Lymantria dispar (Lin et al., 2015), Aedes aegypti (DeGennaro et al., 2013), Microplitis mediator (Li et al., 2012), A. lucorum (Zhou et al., 2014), and Bactrocera dorsalis (Zheng et al., 2012). Due to the crucial role in olfactory perception, Orco is known as an excellent target for investigating co-evolution across sibling insect species (Lu et al., 2009).
The plant mirid species of Lygus spp., Adelphocoris spp., and other species strongly rely on olfactory cues to regulate their chemical perception behaviors. Series of studies on chemoreception of plant mirids were reported such as antennal morphological and electrophysiological characteristic (Chinta et al., 1997; Sun et al., 2014a), putative odorants (Koczor et al., 2012; Sun et al., 2013), physiological functions of OBPs (Gu et al., 2011; Hull et al., 2014; Sun et al., 2014b) and conventional ORs (Yan et al., 2015; An et al., 2016; Xiao et al., 2016; Zhang et al., 2016). In the current study, we focused on the evolutionary divergence of Orco orthologs among plant bug species from distinct geographic regions of China. Three Orco genes from A. suturalis, A. fasciaticollis and L. pratensis are were newly identified. Gene structures, substitution rates and tissues-biased expression of Orco orthologs from five bug species were investigated to further figure out the evolutionary divergence in different mirid bugs.
Materials and Methods
Insect Collection and Rearing
Five mirid bug species including A. lucorum, L. pratensis, A. lineolatus, A. suturalis and A. fasciaticollis were collected from cotton fields at Langfang (Latitude 39.53°N, Longitude 116.70°E) or Kuerle (Latitude 41.45°N, Longitude 85.48°E) experimental station of the Chinese Academy of Agricultural Sciences. The laboratory colony was kept in 20 × 10 × 6 cm rearing containers and was reared on green beans (Phaseolus vulgaris L.) and a 10% sucrose solution (Lu et al., 2008b). Green beans also served as the oviposition substrate and were changed every other day. Beans containing eggs were subsequently placed in rearing containers lined with filter paper. After the emergence of the nymphs, the individuals were transferred to identical containers that were covered with nylon organdy mesh to allow air circulation. The nymphs were provided with fresh food every 2 d until the emergence of adults. Each container housed approximately 100 nymphs or 60 adults. The laboratory colony was maintained at 29 ± 1°C, 60 ± 5% relative humidity (RH), and 14 h:10 h light: dark (L: D) photoperiod.
RNA Extraction and cDNA Synthesis
Antennae from newly eclosion adults were excised and immediately frozen in liquid nitrogen, then stored at −80°C until use. Total RNA was isolated by Trizol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions. The RNA quantity and integrity were checked using 1.2% agarose gel electrophoresis and a NanoDrop 2000 spectrophotometer (NanoDrop, Wilmington, DE, USA). Total RNA was treated with RQ1 RNase-Free DNase (Promega, Madison, USA) at 37°C for 30 min to remove residual DNA. The cDNAs were synthesized using the Superscript III Reverse Transcriptase system (Invitrogen, Carlsbad, CA).
Gene Cloning and Sequence Analysis
AsutOrco, AfasOrco, and LpraOrco genes were cloned using degenerate primers (Table S1). Each reaction contained 300 ng antennal cDNA and 0.5 units of Ex Taq DNA Polymerase (TaKaRa, Dalian, China). The cycling parameters were: 95°C for 2 min followed by 35 cycles at 94°C for 30 s, 55°C for 30 s, 72°C for 60 s, and final extension at 72°C for 10 min. The PCR product was gel-purified and sub-cloned into the pEASY-T3 vector (TransGen, Beijing, China) and then sequencing validation was performed. The 5′ and 3′ regions of Orco genes were amplified using SMARTer™ RACE cDNA amplification kit (Clontech, Mountain View, CA, USA) using gene-specific primers (GSP) (Table S1). Touchdown PCR was performed as follows: 95°C for 2 min followed by 5 cycles at 94°C for 30 s, 72°C for 2 min; 5 cycles at 94°C for 30 s, 70°C for 30 s, and 72°C for 90 s, 30 cycles at 94°C for 30 s, 68°C for 30 s, and 72°C for 90 s; and a final 10 min incubation at 72°C. The RACE PCR products were sub-cloned into the pEASY-T3 vector (Transgene, Beijing, China) and then sequenced. The full-length Orco genes were confirmed with LA Taq DNA polymerase (Takara, Dalian, China) by PCR using gene-specific primers (Table S1).
The full length Orco sequences were aligned by ClustalX 2.1 and edited by GeneDoc 2.7.0 software. TOPCONS (http://topcons.cbr.su.se/) (Tsirigos et al., 2015) was used to identify the number and location of predicted transmembrane domains. The topology diagrams were constructed using TOPO2 Transmembrane Protein Display by the server at http://www.sacs.ucsf.edu/TOPO2/ (SJ)1.
Gene Structure and Selective Pressure Analysis
Genomic DNAs from antennae were extracted using TIANamp genomic DNA kit (TIANGEN, Beijing, China) followed the manufacturer's instruction. Introns of Orco genes were amplified using specific primers (Tables S2–S4).The neighbor joining tree of Orco gene from various insect species were constructed using MEGA7.0 program with a p-distance model and a pairwise deletion of gaps. Bootstrapping was performed by the re-sampling amino acid positions of 1000 replicates, the synonymous and non-synonymous divergence was analyzed using modified Nei-Gojobori (Jukes-Cantor) (assumed transition/transversion bias = 1.21) method in MEGA 7.0 (Jukes and Cantor, 1969; Zhang et al., 1998; Kumar et al., 2016).
Quantitative Real-Time PCR (qPCR) Measurement
The expressions profiles of Orco gene in different tissues of both genders were evaluated by using qPCR measurement on an ABI Prism 7,500 Fast Detection System (Applied Biosystems, Carlsbad, CA, USA).The reference genes β-actin (GenBank accession number: GQ477013, KU230353, KF921006, KU188517, and MG397129, separately) were used as the endogenous control to normalize the target gene expression and correct for any sample-to-sample variation. The primers (Table S5) of the target and reference genes were designed by BEACON DESIGNER 7 (PREMIER Biosoft International). The specificity of each primer set was validated by melt-curve analysis, and the efficiency was calculated by analyzing standard curves with a five-fold cDNA dilution series. Each qPCR reaction was conducted in 20 μL mixture containing 10 μL of 2 × Super-Real PreMix Plus (TIANGEN, Beijing, China), 0.6 μL of each primer (10 μM), 0.4 μL of 50 × Rox Reference Dye, 1 μL of sample cDNA and 7.4 μL of sterilized H2O. The qPCR cycling parameters consisted of 95°C for 15 min, followed by 40 cycles of 95°C for 10 s and 62°C for 30 s, and melt curve stages at 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s. The experiments for the test samples, endogenous control and negative control were performed in triplicate to ensure reproducibility. The comparative 2−ΔΔCT method was used to calculate the relative transcript levels in each tissue samples (Livak and Schmittgen, 2001). All of the data were normalized to endogenous β-actin levels from the same tissue samples.
Results and Discussion
Cloning and Sequence Analysis of Orcos
Among the five plant bug species, two Orcos, AlinOrco from A. lineolatus and AlucOrco from A. lucorum were identified in our previous work (Zhou et al., 2014; Xiao et al., 2016). Here, we focused on other three Orco genes from A. suturalis, A. fasciaticollis, and L. pratensis. The rest three Orco genes were obtained by homology-based cloning (Hull et al. 2012) using degenerate primers (Table S1). A 400 bp fragment encoding putative Orco was amplified from A. fasciaticollis, A. suturalis, and L. pratensis, respectively. The remaining 5′and 3′ end sequences were further obtained using RACE PCR using gene specific primers. Finally, three full length sequences encoded AfasOrco, AsutOrco, and LpraOrco were assembled and deposited in GenBank with the accession numbers MF153393, MF153394, and MF153395, separately. The open reading frames (ORFs) of AsutOrco, AfasOrco, and LpraOrco were 1416, 1416, and 1422 bp, respectively, which resembled the full length of previously reported Orco genes (Hull et al., 2012; An et al., 2016; Xiao et al., 2016).
Results of sequence alignment indicated that all five Orcos including AfasOrco, AsutOrco, LpraOrco, AlinOrco and AlucOrco were rather conserved across the species (Figure 1). The amino acid identity among species of genus Adelphocoris and even across the genera of Adelphocoris, Lygus, and Apolygus was up to 99.6 and 96.8 %, respectively (Table S6). Unlike highly conventional ORs (Clyne et al., 1999; Gao and Chess, 1999), alignment of 200 Orco amino acid sequences (Table S7) from 8 orders showed a 62.6% identity (data not shown). These findings coincide with the previous point of view that Orco is highly conserved (Krieger et al., 2003; Melo et al., 2004; Briguad et al., 2009; Zhao et al., 2013).
Figure 1. Sequence alignment of Orcos from five mirid bug species. Amino acid sequences are aligned by ClustalX 2.1 and edited by GeneDoc 2.7.0 software; the predicted positions of seven putative transmembrane domains (TM1-7) are marked with red transverse line.
Generally, different regions in the gene may play different roles. A predicted algorithm based on TOPCONS revealed these five Orco shared a similar atypical seven trans-membrane topology with their N-terminus inside the cell and the C-terminus exposed to the external environment (Figure 1 and Figure S1). Consequently, the full Orco sequences can be divided into 15 regions, including the intracellular N terminal region, the seven transmembrane regions, the three intracellular loops, the three extracellular loops, and the C terminal region. These data were also consistent with the previous reports (Carraher et al., 2012; Missbach et al., 2014). The amino acid variation among different regions was significantly different with the highest variable level observed at transmembrane regions TM3 and intracellular loop 2 (IL2) that could be involved in ligands binding (Chao et al., 1999; Capendeguy et al., 2006). While no variation was found at intracellular loop 3 (IL3), TM7 and C terminus (Figures S1, S2). It was reported that IL3 participates in the channel activation interaction between conventional ORs and Orco in D. melanogaster and (Benton et al., 2006; Turner et al., 2014). As a key residue, the conserved aspartic acid in TM7 could influence the responses of Orco hemomeric and homomeric ion channels to agonist VUAA1 and odors (Kumar et al., 2013).
Gene Structures of Orcos From Five Bug Species
Introns in Orco genes from different bug species were distinct and sequence identity of at the same position across five bug species was extremely low (about 41 %) in comparison to rather conserved amino acids (data not shown) (Figure 2). Orco within genus Adelphocoris shared similar seven exons, six introns and their insertion loci suggesting the most closely relationships among the three bug species. AlucOrco also had seven exons and six introns, but the length of each intron was significantly larger than that of corresponding introns from genus Adelphocoris. Moreover, the insertion sites of third and fourth introns were also different from Orco in genus Adelphocoris (Figures 2A,B). Notably, only six exons and five introns were found in LpraOrco gene, the last intron of which was lost and the third intron was located between Glu244 and Leu245.
Figure 2. Gene structure and intron insertion loci of five Orcos. (A) Location of extrons (orange rectangles) and introns (blue line) in different Orco genes. (B) Insertion loci labeled using black triangle of different introns in Orco sequence.
Generally, the more intron number and larger intron length indicate a higher phylogenetic level (Nixon et al., 2002; Koonin, 2006; Wu et al., 2013; Park et al., 2014). Adult A. lucorum displays the most extensive distribution in China, whereas L. pratensis mainly occurred in Xinjiang Uygur Autonomous Region (Jiang et al., 2015). The host range is consistence with phylogenetic level among the five mirid bug species; A. lucorum has the widest host range including 54 families, however, L. pratensis merely owns 21 families (Jiang et al., 2015). Additionally, adult A. lucorum prefers to track better host plant food during different seasons than that of other four bug species (Pan et al., 2015; Wang et al., 2017). Likewise, olfaction especially the OR family is believed to play essential roles in the host selection for mirid bugs plant (Yan et al., 2015; Zhang et al., 2016). Therefore, our analyses indicate there might be a potential association between Orco evolution rate and the ecological adaption among these five mirid species, which could contribute to clarify the molecular mechanisms of evolutionary interactions between mirid bugs and their host plants. However, this speculation still needs to be proved by more evidences.
Evolution Analysis of Orco Orthologs
There was a clear conserved orthologous relationship among AfasOrco, AsutOrco, LpraOrco, and other four bug Orcos (AlucOrco, AlinOrco, LlinOrco, LhesOrco). Phylogenetic relationship was largely consistent with the species tree constructed from the alignment of species-specific cytochrome oxidase subunit I (COI) (Figures 3A,B). So, we suggested that Orco was significantly conserved and could function as a molecular marker of evolution across bug species. Also relatedness analysis of these seven Orcos to the other 193 Orco sequences from eight insect orders indicated that Orco was highly conserved within insect order. Orco sequences of the same order were strictly clustered together with strong bootstrapping support (Figure 3C), indicating this phylogenetic clade was highly conserved and may fulfill conserved function.
Figure 3. Neighbor joining tree of Orcos from different insect species. (A) Phylogenetic tree of Orcos from seven bug species. (B) Phylogenetic relationships among seven species constructed using species-specific cytochrome oxidase subunit I (COI). (C) Phylogenetic tree of insect Orcos from different orders. Yellow triangle, Coleoptera; Red triangle, Hymenoptera; Dark blue solid diamond, Diptera; Dark green solid square, Lepidoptera; Black solid circle, Hemiptera; Light blue hollow circle, Orthoptera; Light blue hollow diamond, Blattaria; Light red hollow circle, Anoplura.
The ratios of non-synonymous to synonymous substitutions estimated for 14 Orco genes from 5 orders were listed in Table 1. All the ratios were far less than 1.0 indicating that Orco genes are under strong purifying selection pressure. The strong purifying selection pressure suggested a functional conservation, which had been proven by substantial documents. The lack of Orco leading to a similar reduction of olfaction indicates the consistent roles in odor perceptions, suggesting the interspecific conservation of Orco indirectly (Zhou et al., 2014; Liu et al., 2017; Trible et al., 2017). Furthermore, the interspecific functional conservation has been confirmed directly by transgenic rescue experiment. The defects of olfaction in DmelOrco mutant could be rescued by transgenic expression of DmelOrco, CcapOrco, AgamOrco and HzeaOrco, respectively (Jones et al., 2005). It was also demonstrated that Orco, as an obligatory part of ligand-gated ion channel, played conservative functions in ligand binding, and was activated by the agonist VUAA1 dutifully (Benton et al., 2006; Sato et al., 2008; Jones et al., 2011).
Orcos are under strongly purifying selection pressure and exhibit potential conserved olfaction roles. However, our estimation on the ratios of non-synonymous to synonymous substitutions (Ka/Ks) revealed that their levels of the purifying selection pressure significantly varied in the genera and species. Generally, the values of Ka/Ks were similar among species within same genus, but were different from species across genera. As shown in Table 1, when used DmelOrco, AgamOrco or Orco genes from other model species as outgroup, the range of Ka/Ks values of Orco genes from three Adelphocoris species were evaluated as (0.288–0.454), which was significantly different to that of AlucOrco (0.365–0.496) from Apolygus genus, or LpraOrco, LlinOrco and LhesOrco (0.361–0.497) from Lygus genus. These findings indicted there might be a strong constraint on functional variation within Orco from same genus, as illustrated above. In addition, these results were faultlessly correlated to the phylogenetic analyses (Figure 3). Three Orco genes from Adelphocoris species fall into the same clade, AlucOrco and three Orco from Lygus species clustered in another clade. Because of the evolutionary synchronization between Orco genes and their mirid species (Figures 3A,B), we proposed that the degrees of variation (suggested by Ka/Ks values) on Orco protein coding regions could reflected the phylogenetic levels of mirid bug species, and our data would lay a foundation on the further studies on the molecular mechanisms of speciation of mirid bugs.
Expression Profiles of Five Orcos
In general, target gene with different tissue expressions would play different physiological function. To figure out the potential roles of Orco in mirid bugs species, qPCR measurement was conducted to assess their tissue-specific expressions (Figure 4). The results demonstrated that these five Orco genes share similar antennae-biased expression profiles, which were similar to that in L. hesperus (Hull et al., 2012). So, we suspected that Orco in different mirid bugs could be associated with clear olfactory roles. It was reported that silencing in A. lucorum of the olfactory co-receptor Orco gene by RNA interference could induce EAG response declining to two putative semiochemicals (Zhou et al., 2014). However, some Orco could be also expressed in non-olfactory organs such as proboscis and legs, suggesting that Orco might be involved in the contact chemosensory perception and could help to search hosts in close distance and perceive the status of hosts (Lu et al., 2009; Hull et al., 2012). In this study, faint transcript levels of these five Orcos were detected in stylets, legs, head and other non-olfactory organs (Figure 4) suggesting the potential roles of Orco in taste recognition of bugs. Besides in mirid bug species, Orco of B. dorsalis could fulfill a role involved in the perception of Rhodojaponin-III, a non-volatile compound (Yi et al., 2014).
Figure 4. Orco expressions in different tissues of five mirid bug species. The error bars represent standard error, and different letters above each bar denote significant differences (P < 0.05).
Y-JZ and LS conceived and designed the experimental plan. QiW, Y-LZ, and H-HC performed the experiments. QiW, QianW, SS, YX, KD, and AK analyzed the data. QiW and QianW drafted the manuscript. LS and Y-JZ refined and approved the final manuscript.
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.
Dr. Xianhui Wang (Institute of Zoology, Chinese Academy of Sciences, Beijing, China) and Dr. Bo Zhang (Institute of Plant Protection, Chinese Academy of Agricultural Sciences) kindly helped analyse the data on evolution. This work was supported by the China National Basic Research Program (2012CB114104), the National Natural Science Foundation of China (31471778, 31501652,31621064 and 31772176), Central public-interest Scientific Institution Basal Research Fund (No. 1610212016015), Research Foundation of State Key Laboratory for Biology of Plant Diseases and Insect Pests (SKLOF201514 and SKLOF201719), and The Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2015-TRICAAS).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00158/full#supplementary-material
Figure S1. Predicted transmembrane topologies of Orco with variable sites highlighted in red. The double line indicates the membrane region with extracellular and cytoplasmic sides labeled.
Figure S2. Ratio of the relative amino acid differences per domain averaged for Orco.
Table S1. Primers used in identification of Orco genes from mirid bug A. suturalis, A. fasciaticollis, and L. pratensis.
Table S2. Primers used in amplification of AlucOrco gene introns.
Table S3. Primers used in amplification of AlinOrco, AfasOrco and AsutOrco gene introns.
Table S4. Primers used in amplification of LpraOrco gene introns.
Table S5. Primers used in qPCR measurement.
Table S6. Identities of amino acid sequences among Orco genes from five mirid bugs.
Table S7. Orcos used in phyologenetic construction and sequence analysis.
An, X. K., Sun, L., Liu, H. W., Liu, D. F., Ding, Y.-X., Li, L. M., et al. (2016). Identification and expression analysis of an olfactory receptor gene family in green plant bug Apolygus lucorum (Meyer-Dür). Sci. Rep. 6:37870. doi: 10.1038/srep37870
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
Briguad, I., Montagné, N., Monsempes, C., François, M., and Jacquin-Joly, E. (2009). Identification of an atypical insect olfactory receptor subtype highly conserved with noctuids. FEBS J. 276, 6537–6547. doi: 10.1111/j.1742-4658.2009.07351.x
Capendeguy, O., Chodanowski, P., Michielin, O., and Horisberger, H. D. (2006). Access of extracellular cations to their binding sites in Na, K-ATPase: role of the second extracellular loop of the α subunit. J. Gen. Physiol. 127, 341–352. doi: 10.1085/jgp.200509418
Carraher, C., Authier, A., Steinwender, B., and Newcomb, R. D. (2012). Sequence comparisons of odorant receptors among Tortricid moths reveal different rates of molecular evolution among family members. PLoS ONE 7:e38391. doi: 10.1371/journal.pone.0038391
Chao, T. H., Ember, J. A., Wang, M., Bayon, Y., Hugli, T. E., and Ye, R. D. (1999). Role of the second extracellular loop of human C3a receptor in agonist binding and receptor function. J. Biol. Chem. 274, 9721–9728. doi: 10.1074/jbc.274.14.9721
Chinta, S., Dickens, J. C., and Baker, G. T. (1997). Morphology and distribution of antennal sensilla of the tarnished plant bug, Lygus lineolaris (Palisot de beauvois) (Hemiptera: Miridae). Int. J. Insect Morphol. 26, 21–26. doi: 10.1016/S0020-7322(96)00022-0
Clyne, P. J., Warr, C. G., Freeman, M. R., Lessing, D., Kim, J., and Carlson, J. R. (1999). A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron 22, 327–338. doi: 10.1016/S0896-6273(00)81093-4
DeGennaro, M., McBride, C. S., Seeholzer, L., Nakagawa, T., Dennis, E. J., Goldman, C., et al. (2013). Orco mutant mosquitoes lose strong preference for humans and are not repelled by volatile DEET. Nature 498, 487–491. doi: 10.1038/nature12206
Dong, X., Zhong, G., Hu, M., Yi, X., Zhao, H., and Wang, H. (2013). Molecular cloning and functional identification of an insect odorant receptor gene in Spodoptera litura (F.) for the botanical insecticide rhodojaponin III. J. Insect Physiol. 59, 26–32. doi: 10.1016/j.jinsphys.2012.11.004
Gu, S. H., Wang, W. X., Wang, G. R., Zhang, X. Y., Guo, Y. Y., Zhang, Z., et al. (2011). Functional characterization and immunolocalization of odorant binding protein 1 in the lucerne plant bug, Adelphocoris lineolatus (Goeze). Arch. Insect Biochem. Physiol. 77, 81–99. doi: 10.1002/arch.20427
Hull, J. J., Hoffmann, E. J., Perera, O. P., and Snodgrass, G. L. (2012). Identification of the western tarnished plant bug (Lygus hesperus) olfactory co-receptor orco: expression profile and confirmation of atypical membrane topology. Arch. Insect Biochem. Physiol. 81, 179–198. doi: 10.1002/arch.21042
Hull, J. J., Perera, O. P., and Snodgrass, G. L. (2014). Cloning and expression profiling of odorant-binding proteins in the tarnished plant bug, Lygus lineolaris. Insect Mol. Biol. 23, 78–97. doi: 10.1111/imb.12064
Jones, P. L., Pask, G. M., Rinker, D. C., and Zwiebel, L. J. (2011). Functional agonism of insect odorant receptor ion channels. Proc. Natl. Acad. Sci. U.S.A. 108, 8821–8825. doi: 10.1073/pnas.1102425108
Jones, P. L., Pask, G. M., Romaine, I. M., Taylor, R. W., Reid, P. R., Waterson, A. G., et al. (2012). Allosteric antagonism of insect odorant receptor ion channels. PLoS ONE 7:e30304. doi: 10.1371/journal.pone.0030304
Jones, W. D., Nguyen, T. T., Kloss, B., Lee, K. J., and Vosshall, L. B. (2005). Functional conservation of an insect odorant receptor gene across 250 million years of evolution. Curr. Biol. 15, R119–R121. doi: 10.1016/j.cub.2005.02.007
Koczor, S., Vuts, J., and Tóth, M. (2012). Attraction of Lygus rugulipennis and Adelphocoris lineolatus to synthetic floral odour compounds in field experiments in Hungary. J. Pest Sci. 85, 239–245. doi: 10.1007/s10340-012-0422-5
Krieger, J., Klink, O., Mohl, C., Raming, K., and Breer, H. (2003). A candidate olfactory receptor subtype highly conserved across different insect orders. J. Comp. Physiol. A Neuroethol. Sens. Neural. Behav. Physiol. 189, 519–526. doi: 10.1007/s00359-003-0427-x
Kumar, B. N., Taylor, R. W., Pask, G. M., Zwiebel, L. J., Newcomb, R. D., and Christie, D. L. (2013). A conserved aspartic acid is important for agonist (VUAA1) and odorant/tunning receptor-dependent activation of the insect odorant co-receptor (Orco). PLoS ONE 8:e70218. doi: 10.1371/journal.pone.0070218
Larsson, M. C., Domingos, A. I., Jones, W. D., Chiappe, M. E., Amrein, H., and Vosshall, L. B. (2004). Or83b encodes a broadly expressed odorant receptor essential for Drosophila olfaction. Neuron 43, 703–714. doi: 10.1016/j.neuron.2004.08.019
Li, K. M., Ren, L. Y., Zhang, Y. J., Wu, K. M., and Guo, Y. Y. (2012). Knockdown of Microplitis mediator odorant receptor involved in the sensitive detection of two chemicals. J. Chem. Ecol. 38, 287–294. doi: 10.1007/s10886-012-0085-y
Li, Y., Zhang, J., Chen, D., Yang, P., Jiang, F., Wang, X., et al. (2016). CRISPR/Cas9 in locusts: successful establishment of an olfactory deficiency line by targeting the mutagenesis of an odorant receptor co-receptor (Orco). Insect Biochem. Mol. Biol. 79, 27–35. doi: 10.1016/j.ibmb.2016.10.003
Lin, W., Yu, Y., Zhou, P., Zhang, J., Dou, L., Hao, Q., et al. (2015). Identification and knockdown of the olfactory receptor (Orco) in gypsy moth, Lymantria dispar. Int. J. Biol. Sci. 11, 772–780. doi: 10.7150/ijbs.11898
Liu, Q., Liu, W., Zeng, B., Wang, G., Hao, D., and Huang, Y. (2017). Deletion of the Bombyx mori odorant receptor co-receptor (BmOrco) impairs olfactory sensitivity in silkworms. Insect Biochem. Mol. Biol. 86, 58–67. doi: 10.1016/j.ibmb.2017.05.007
Lu, B., Wang, N., Xiao, J., Xu, Y., Murphy, R. W., and Huang, D. (2009). Expression and evolutionary divergence of the non-conventional olfactory receptor in four species of fig wasp associated with one species of fig. BMC Evol. Biol. 9:43. doi: 10.1186/1471-2148-9-43
Lu, Y. H., Jiao, Z. B., Li, G. P., Wyckhuys, K. A. G., and Wu, K. M. (2011). Comparative overwintering host range of three Adelphocoris species (Hemiptera: Miridae) in northern China. Crop Prot. 30, 1455–1460. doi: 10.1016/j.cropro.2011.07.010
Lu, Y. H., Qiu, F., Feng, H. Q., Li, H. B., Yang, Z. C., Wyckhuys, et al. (2008a). Species composition and seasonal abundance of pestiferous plant bugs (Hemiptera: Miridae) on Bt cotton in China. Crop Prot. 27, 465–472. doi: 10.1016/j.cropro.2007.07.017
Lu, Y., Wu, K., Jiang, Y., Xia, B., Li, P., Feng, H., et al. (2010). Mirid bug outbreaks in multiple crops correlated with wide-scale adoption of Bt cotton in China. Science 328, 1151–1154. doi: 10.1126/science.1187881
Lundin, C., Käll, L., Kreher, S. A., Kapp, K., Sonnhammer, E. L., Carlson, J. R., et al. (2007). Membrane topology of the Drosophila OR83b odorant receptor. FEBS Lett. 581, 5601–5604. doi: 10.1016/j.febslet.2007.11.007
Melo, A. C. A., Rützler, M., Pitts, R. J., and Zwiebel, L. J. (2004). Identification of a chemosensory receptor from the yellow fever mosquito, Aedes aegypti, that is highly conserved and expressed in olfactory and gustatory organs. Chem. Senses 29, 403–410. doi: 10.1093/chemse/bjh041
Nichols, A. S., and Luetje, C. W. (2010). Transmembrane segment 3 of Drosophila melanogaster odorant receptor subunit 85b contributes to ligand-receptor interactions. J. Biol. Chem. 285, 11854–11862. doi: 10.1074/jbc.M109.058321
Nichols, A. S., Chen, S., and Luetje, C. W. (2011). Subunit contributions to insect olfactory receptor function: channel block and odorant recognition. Chem. Senses 36, 781–790. doi: 10.1093/chemse/bjr053
Nixon, J. E., Wang, A., Morrison, H. G., McArthur, A. G., Sogin, M. L., Loftus, B. J., et al. (2002). A spiliceosomal intron in Giardialamblia. Proc. Natl. Acad. Sci. U.S.A. 99, 3701–3705. doi: 10.1073/pnas.042700299
Pan, H., Liu, B., Lu, Y., and Wyckhuys, K. A. (2015). Seasonal alterations in host range and fidelity in the polyphagous mirid bug, Apolygus lucorum (Heteroptera: Miridae). PLoS ONE 10:e0117153. doi: 10.1371/journal.pone.0117153
Park, S. G., Hannenhalli, S., and Choi, S. S. (2014). Conservation in first intron is positively associated with the number of exons within genes and the presence of regulatory epigenetic signals. BMC Genomics 15:526. doi: 10.1186/1471-2164-15-526
Pask, G. M., Bobkov, Y. V., Corey, E. A., Ache, B. W., and Zwiebel, L. J. (2013). Blockade of insect odorant receptor currents by amiloride derivatives. Chem. Senses 38, 221–229. doi: 10.1093/chemse/bjs100
Sato, K., Pellegrino, M., Nakagawa, T., Nakagawa, T., Vosshall, L. B., and Touhara, K. (2008). Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature 452, 1002–1006. doi: 10.1038/nature06850
Sun, L., Gu, S. H., Xiao, H. J., Zhou, J. J., Guo, Y. Y., Liu, Z. W., et al. (2013). The preferential binding of a sensory organ specific odorant binding protein of the alfalfa plant bug Adelphocoris lineolatus AlinOBP10 to biologically active host plant volatiles. J. Chem. Ecol. 39, 1221–1231. doi: 10.1007/s10886-013-0333-9
Sun, L., Xiao, H. J., Gu, S. H., Guo, Y. Y., Liu, Z. W., and Zhang, Y. J. (2014a). Perception of potential sex pheromones and host-associated volatiles in the cotton plant bug, Adelphocoris fasciaticollis (Hemiptera: Miridae): morphology and electrophysiology. Appl. Entomol. Zool. 49, 43–57. doi: 10.1007/s13355-013-0223-1
Sun, L., Xiao, H. J., Gu, S. H., Zhou, J. J., Guo, Y. Y., Liu, Z. W., et al. (2014b). The antenna-specific odorant-binding protein AlinOBP13 of the alfalfa plant bug Adelphocoris lineolatus is expressed specifically in basiconic sensilla and has high binding affinity to terpenoids. Insect Mol. Biol. 23, 417–434. doi: 10.1111/imb.12089
Trible, W., Olivos-Cisneros, L., McKenzie, S. K., Saragosti, J., Chang, N. C., Matthews, B. J., et al. (2017). Orco mutagenesis causes loss of antennal lobe glomeruli and impaired social behavior in ants. Cell 170, 727–735. doi: 10.1016/j.cell.2017.07.001
Tsirigos, K. D., Peters, C., Shu, N., Käll, L., and Elofsson, A. (2015). The TOPCONS web server for consensus prediction of membrane protein topology and signal peptides. Nucleic Acids Res. 43, W401–W407. doi: 10.1093/nar/gkv485
Turner, R. M., Derryberry, S. L., Kumar, B. N., Brittain, T., Zwiebel, L. J., Newcomb, R. D., et al. (2014). Mutational analysis of cysteine residues of the insect odorant co-receptor (Orco) from Drosophila melanogaster reveals differential effects on agonist- and odorant-tuning receptor-dependent activation. J. Biol. Chem. 289, 31837–31845. doi: 10.1074/jbc.M114.603993
Wang, Q., Bao, W. F., Yang, F., Xu, B., and Yang, Y. Z. (2017). The specific host plant DNA detection suggests a potential migration of Apolygus lucorum from cotton to mungbean fields. PLoS ONE 12:e0177789. doi: 10.1371/journal.pone.0177789
Wicher, D., Schäfer, R., Bauernfeind, R., Stensmyr, M. C., Heller, R., Heinemann, S. H., et al. (2008). Drosophila odorant receptors are both ligand-gated and cyclic-nucleotide-activated cation channels. Nature 452, 1007–1011. doi: 10.1038/nature06861
Wu, J., Xiao, J., Wang, L., Zhong, J., Yin, H., Wu, S., et al. (2013). Systematic analysis of intron size and abundance parameters in diverse lineages. Sci. China Life Sci. 56, 968–974. doi: 10.1007/s11427-013-4540-y
Wu, K., Lu, Y., Feng, H., Jiang, Y., and Zhao, J. Z. (2008). Suppression of cotton bollworm in multiple crops in China in areas with Bt toxin–containing cotton. Science 321, 1676–1678. doi: 10.1126/science.1160550
Xiao, Y., Sun, L., Ma, X. Y., Dong, K., Liu, H. W., Wang, Q., et al. (2016). Identification and characterization of the distinct expression profiles of candidate chemosensory membrane proteins in the antennal transcriptome of Adelphocoris lineolatus (Goeze). Insect Mol. Biol. 26, 74–91. doi: 10.1111/imb.12272
Yan, S. W., Zhang, J., Liu, Y., Li, G. Q., and Wang, G. R. (2015). An olfactory receptor from Apolygus lucorum (Meyer-Dür) mainly tuned to volatiles from flowering host plants. J. Insect Physiol. 79, 36–41. doi: 10.1016/j.jinsphys.2015.06.002
Yi, X., Zhao, H., Wang, P., Hu, M., and Zhong, G. (2014). BdorOrco is important for oviposition-deterring behavior induced by both the volatile and non-volatile repellents in Bactrocera dorsalis (Diptera: Tephritidae). J. Insect Physiol. 65, 51–56. doi: 10.1016/j.jinsphys.2014.05.007
Zhang, J., Rosenberg, H. F., and Nei, M. (1998). Positive Darwinian selection after gene duplication in primate ribonuclease gene. Proc. Natl. Acad. Sci. U.S.A. 95, 3708–3713. doi: 10.1073/pnas.95.7.3708
Zhang, Z., Zhang, M., Yan, S., Wang, G., and Liu, Y. (2016). A female-biased odorant receptor from Apolygus lucorum (Meyer-Dür) tuned to some plant odors. Int. J. Mol. Sci. 17:1165. doi: 10.3390/ijms17081165
Zhao, H., Gao, P., Zhang, C., Ma, W., and Jiang, Y. (2013). Molecular identification and expressive characterization of an olfactory co-receptor gene in the Asian honeybee, Apis cerana cerana. J. Insect Sci. 13, 1–14. doi: 10.1673/031.013.8001
Zheng, W., Zhu, C., Peng, T., and Zhang, H. (2012). Odorant receptor co-receptor Orco is upregulated by methyl eugenol in male Bactrocera dorsalis (Diptera: Tephritidae). J. Insect Physiol. 58, 1122–1127. doi: 10.1016/j.jinsphys.2012.05.011
Zhou, Y. L., Zhu, X. Q., Gu, S. H., Cui, H. H., Guo, Y. Y., Zhou, J. J., et al. (2014). Silencing in Apolygus lucorum of the olfactory coreceptor Orco gene by RNA interference induces EAG response declining to two putative semiochemicals. J. Insect Physiol. 60, 31–39. doi: 10.1016/j.jinsphys.2013.10.006
Keywords: olfactory receptor co-receptor, mirid bugs, gene structure, sequence analysis, evolution analysis
Citation: Wang Q, Wang Q, Zhou Y-L, Shan S, Cui H-H, Xiao Y, Dong K, Khashaveh A, Sun L and Zhang Y-J (2018) Characterization and Comparative Analysis of Olfactory Receptor Co-Receptor Orco Orthologs Among Five Mirid Bug Species. Front. Physiol. 9:158. doi: 10.3389/fphys.2018.00158
Received: 14 November 2017; Accepted: 16 February 2018;
Published: 05 March 2018.
Edited by:Peng He, Guizhou University, China
Reviewed by:Ya-Nan Zhang, Huaibei Normal University, China
Yihan Xia, Lund University, Sweden
Joe Hull, Agricultural Research Service (USDA), United States
Copyright © 2018 Wang, Wang, Zhou, Shan, Cui, Xiao, Dong, Khashaveh, Sun and Zhang. 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 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.