Silencing the odorant receptor co-receptor impairs olfactory reception in a sensillum-specific manner in the cockroach

Summary Insects detect odors via a large variety of odorant receptors (ORs) expressed in olfactory sensory neurons (OSNs). The insect OR is a heteromeric complex composed of a ligand-specific receptor and the co-receptor (ORco). In this study, we identified the ORco gene of the cockroach, Periplaneta americana (PameORco), and performed RNAi-based functional analysis of PameORco. All OSNs in the basiconic sensilla expressed PameORco and received a large variety of odors including sex pheromones. In trichoid sensilla, a PameORco-positive OSN was consistently paired with a PameORco-negative OSN tuned to acids. In adult cockroaches injected with PameORco dsRNA at the nymphal stage, the expression of PameORco, odor receptions via ORs, and its central processing were strongly suppressed. These results provide new insights into the molecular basis of olfactory reception in the cockroach. The long-lasting and irreversible effects of PameORco RNAi would be an effective method for controlling the household pest.

Whole sequence of ORco in the American cockroach (PameORco) was characterized PameORco expressed in olfactory sensory neurons in a sensillar type-specific manner RNAi chronically and irreversibly suppressed the PameORco expression beyond molts PameORco was essential for receptions of sex pheromones and general odors

INTRODUCTION
Various species of cockroaches including the American cockroach Periplaneta americana are common and widespread pests as vectors of various diseases (Do et al., 2016;Pai et al., 2003). Integrated pest management (IPM), which is a comprehensive strategy for pest control that causes less environmental damage, has received increased attention in recent years. Chemical insecticides and toxic baits have often been used to control cockroaches, but they can also harm non-target organisms. In addition, it is reported in German cockroaches that the excessive use of toxic baits can result in the evolution of pesticide-resistance (Wada-Katsumata et al., 2013). Thus, IPM-based strategies are needed to effectively control the pest cockroaches (Gondhalekar et al., 2021). Because nocturnal cockroaches heavily rely on olfactory cues to search for food and mates, understanding the olfactory system and odor-guided behaviors is essential for optimizing IPM programs for pest cockroaches.
The simple but sophisticated olfactory system in the American cockroach is one of the best models for developing an olfaction-based IPM strategy. To determine how the cockroach encodes vast odors and translates them into odor-guided behaviors, the neural processing of odorants, including pheromones, in P. americana has been studied at the cellular levels ranging from peripheral to higher brain centers (Boeckh and Ernst, 1987;Fujimura et al., 1991;Fusca and Kloppenburg, 2021;Nishino et al., 2018;Sass, 1978Sass, , 1983Takahashi et al., 2019;Tichy and Hellwig, 2018;Tichy et al., 2020;Watanabe et al., 2012Watanabe et al., , 2017. However, the functions of olfactory receptors that confer the response selectivity to a given odor have yet to be identified in the cockroach. In insects, odorants are generally detected by olfactory sensory neurons (OSNs) in the olfactory sensilla, which are widely distributed in the antennae and maxillary palps (Watanabe et al., 2012;2018a). Within the olfactory sensilla, odor molecules diffused into the sensillar lymph are received by olfactory receptors expressed on the dendritic cilia of OSNs. The olfactory receptors in insects are classified into two main types belonging to distinct receptor families: odorant receptors (ORs) belonging to the seven-transmembrane type receptor family, and the ionotropic receptors (IRs) belonging to the ionotropic glutamate receptor family (Wicher and Miazzi, 2021). Both ORs and IRs form heteromeric ligand-gated cation channels composed of ligand-specific receptor proteins (ORx or IRx) and their co-receptor proteins (ORco or IRco), peripheral PB-SNs, central PB-responsive PNs where a large number of PB-SNs converge, and sexual behaviors to PB. This is the first study to reveal the fundamental molecular process of odor transduction in the olfaction of the American cockroach, the animal model for IPM.

Identification of PameORco
A tbalstn search revealed that a single ORco gene is encoded by the genome of P. americana. The coding region of the ORco gene is encoded by seven exons spanning $6.5 kb in the genome (GenBank ID: base 91,777-33,081 of PGRX01009403.1; Figure 1A). Next, we amplified 1444-bp and 1495-bp cDNA fragments encoding two PameORco isoforms (471 a.a. and 472 a.a., respectively) in the cockroach antennae using reverse transcription PCR (RT-PCR)-based cDNA cloning (GenBank ID: LC657818 and LC657819). The Pa-meORco isoforms contain seven transmembrane helixes and putative phosphorylation sites (Guo et al., 2017), a putative calmodulin-binding site (Jain et al., 2021), and a putative N-linked glycosylation site iScience Article (Lundin et al., 2007) conserved in insect ORco homologs ( Figure 1B). The two PameORco isoforms share the N-terminal 456 a.a. and differ at their C-terminal end corresponding to the C-terminal half of the seventh transmembrane helix. The C-terminal tail of PameORco isoform 1 contains multiple hydrophilic residues, which is not predicted as transmembrane region by the TMHMM2.0 algorithm. With these data, we concluded that the PameORco isoform 1 is a missplicing product of the PameORco gene, and does not encode functional PameORco protein. Therefore, we focused on the PameORco isoform 2 in this study. A molecular phylogenetic analysis of insect ORco homologs confirmed that PameORco is closely related to the ORco homologs of other hemimetabolous insect species ( Figure S1).

Anti-PameORco antiserum selectively labeled OSNs in specific olfactory sensilla
To examine the distribution of PameORco-expressing OSNs in the cockroach antennae, we generated anti-PameORco antiserum in a guinea pig and conducted immunohistochemistry. The antennae of adult male cockroaches were double-immunolabeled with the anti-PameORco antiserum and the polyclonal anti-HRP antibody generated in a rabbit. The anti-HRP antibody selectively labels sensory neurons in cockroach antennae . The types of olfactory sensilla were unambiguously identified by differential interference contrast (DIC) optics (Figures 2A1-2D1, Videos S1-S4). We investigated the distributions of PameORco-expressing OSNs in three types of olfactory sensilla; perforated basiconic, trichoid (Figures 2A,2B,and 2C), and grooved basiconic sensilla ( Figure 2D). Perforated basiconic sensilla are further categorized into two subtypes on the basis of their lengths; shorter single-walled A (sw-A) sensilla and sex pheromone-responsive longer sw-B sensilla (Schaller, 1978;Watanabe et al., 2018a). Anti-HRP antibody revealed two and four OSNs in single sw-A and sw-B sensilla, respectively, and these OSNs were also labeled by anti-PameORco antiserum ( Figure 2A). However, OSNs in the grooved basiconic sensilla were not labeled by anti-PameORco antiserum regardless of the subtypes (three OSNs in double-walled A1 (dw-A1) and four OSNs in dw-A2 sensilla in Figure 2D). The longer and shorter trichoid sensilla were termed sw-C1 and sw-C2 sensilla, respectively (Watanabe et al., 2012). Anti-HRP antibody revealed that single trichoid sensilla had two OSNs. Regardless of subtype, a PameORco-positive OSN was consistently paired with a PameORco-negative OSN in single trichoid sensilla ( Figures 2B and 2C). Mechanosensory and contact chemosensory neurons in the chaetic sensilla, and hygrosensory and thermosensory neurons in the capitular sensilla were not labeled by anti-PameORco antiserum ( Figure S2). The distribution patterns of anti-Pa-meORco immunoreactive OSNs are summarized in schemes ( Figure 2E).

PameORco RNAi and its effectiveness
To validate the effectiveness of PameORco RNAi, we conducted RT-qPCR analysis to examine the expression level of PameORco in the antennae at different time points after dsRNA injection ( Figure 3). Nymphal and adult male cockroaches were injected with 4 mg of dsRNA of PameORco or Escherichia coli b-lactamase (b-lac; negative control), and the expression level of PameORco in both antennae was quantified at arbitrarily selected days after dsRNA injection ( Figure 3A). In the RT-qPCR, the expression levels of PameORco were normalized with that of elongation factor 1-alpha in P. americana (PameEf1a; GenBank ID: LC657820). There was a significant reduction of the expression levels of PameORco in PameORco dsRNA-injected cockroaches compared to Naive or b-lac dsRNA-injected cockroaches ( Figures 3B and 3C). Because expression levels of PameORco were not different between b-lac dsRNA-injected and Naive cockroaches ( Figure 3B), the injection procedure did not affect the expression of PameORco. Compared to the Naive cockroaches, the expression levels of PameORco were reduced to 2.1 G 1.5% at 3 days after PameORco dsRNA injection and the strong suppression of PameORco expression was observed across different post-dsRNA injection periods ( Figure 3C). The effects of PameORco RNAi were maintained beyond developmental stages; adult cockroaches injected with PameORco dsRNA at the last instar stage also exhibited the prominent RNAi effect ( Figure 3C).

Silencing of the PameORco expression impaired sex pheromone reception
Because all OSNs in sex pheromone-responsive sw-B sensilla expressed PameORco, we investigated the role of PameORco in the perception of sex pheromones in the American cockroach. At arbitrary selected days after dsRNA injection, we performed single sensillum recordings (SSRs) from single sw-B sensilla and recorded responses to periplanone-A (PA) and -B (PB), which are two components of P. americana sex pheromone ( Figure 4). A sw-B sensillum has one PA-SN and one PB-SN in addition to two OSNs (OSN1 and OSN2), which respond to various general odors ( Figure 4A; Sass, 1978;Boeckh and Ernst, 1987;Fujimura et al., 1991 Figure 2. Sensillar type-specific distribution patterns of PameORco-expressing OSNs (A-D) Anti-PameORco-immunopositive and immunonegative OSNs in all morphological types of antennal olfactory sensilla. Six morphological types of antennal olfactory sensilla are unambiguously identified via differential interference contrast (DIC) observations and indicated by different colored arrowheads (A1-D1). All antennal sensory neurons labeled by anti-HRP antibody are green-colored (A2-D2 and A4-B4), whereas anti-PameORco-immunopositive OSNs are magenta-colored (A3-D3 and A4-B4). Sensilla and corresponding OSNs are denoted as arrowheads and arrows of the same color, respectively. Two OSNs in single trichoid sensilla (sw-C1 and sw-C2 sensilla) are indicated by thin arrows (B2-B4 and C2-C4). In single trichoid sensilla, a PameORco-positive OSN was consistently paired with a PameORco-negative OSN (thin arrows in B3,4 and C3,4). Scale bar = 20 mm. See also Figure S2 and Videos S1, S2, S3, and S4. (E) A schematic overview of sensillar type-specific distribution patterns of PameORco-expressing OSNs.
Because of the difficulty of discriminating small-amplitude spikes from PA-SN and PB-SN according to their spike shapes ( Figures S3A and S3B), PA-evoked responses of PA-SNs and PB-evoked responses of PB-SNs were characterized by performing SSRs with a cross-adaptation odor stimuli system ( Figures S3C and S3D). In brief, PA-SN responded exclusively to PA, whereas the PB-SN responded to both PB and PA ( Figure 4A referred to from Boeckh and Ernst, 1987). Therefore, we recorded the PA response of PA-SN after eliminating the PA response of PB-SN by repeating stimulations with PB ( Figure S3D). Our SSRs revealed that A B C Developmental processes of OSNs and sexual maturation are indicated by green bars (Silverman, 1977;Watanabe et al., 2018a). (B) PameORco expressions in Naive, b-lac, and PameORco dsRNA-injected cockroaches. In dsRNA-injected cockroaches, RT-qPCRs were performed at 7 days after dsRNA injection. The expression levels of PameORco mRNA were normalized with that of PameEf1a mRNA. (C) PameORco expressions across the different post-injection periods of dsRNA. The expression levels of PameORco were normalized to the average of those of Naive cockroaches (baseline expression level). In (B and C), the sample number is noted in parentheses and black bars indicate means. The different letters above each plot indicate significant differences (ANOVA post-hoc Tukey-Kramer test; p < 0.05). Samples from eight Naive cockroaches (black dots in B and C) were specifically prepared for the RT-qPCR analysis, but those obtained from RNAi cockroaches were also used in the following electrophysiological experiments. Therefore, throughout following figures, results obtained from the same individuals are represented as symbols of the same color and shape.   Figures 4C, 4E, and 4F). The PameORco expression was significantly suppressed from 3 days after PameORco dsRNA injection ( Figure 3C), but responses to PA and PB at the time point was still stronger than those at >7 days after injection. Adult cockroaches injected with PameORco dsRNA at the nymphal stage responded to neither PA nor PB. The effects of RNAi were maintained at least 88 days after Pa-meORco dsRNA injection into the 10th instar nymph.
Next, we characterized the sensitivities of PA-SNs and PB-SNs in PameORco RNAi cockroaches (Figures 4G and 4H). In the control cockroaches (Naive or b-lac in Figure 3A), responses of PB-SNs and PA-SNs plateaued at 0.5 ng of PB and 0.2 ng of PA, respectively (Figures 4G and 4H). At 3 days after dsRNA injection, PB-SNs exhibited sufficiently strong responses to high concentrations of PB, but the sensitivity to PB was lower than that in control cockroaches. Adult cockroaches injected with PameORco dsRNA at nymphal stages did not respond to high concentrations of sex pheromones (LI or 10th in Figures 4G and 4H). Our results revealed that sex pheromones in the cockroach are received by functional ORs, and sensitivities and response intensities to sex pheromones depended on the post-dsRNA injection periods. Especially, the effects of PameORco RNAi were long-lasting and most prominent in adult cockroaches injected with PameORco dsRNA at nymphal stages.
Silencing of the PameORco expression impaired general odor reception in a sensillum-specific manner Next, we investigated the role of PameORco in the perception of general odors. We focused on OSNs in the sw-A, sw-C2, and dw-A2 sensilla whose ligands have been identified (Fujimura et al., 1991). Because the spikes were not easily distinguishable in many recorded olfactory sensilla, we characterized recorded sensilla in terms of the total number of spikes per sensillum. Based on the results of the effects of Pa-meORco RNAi on sex pheromone reception, we compared sensillar responses to general odors among the three groups; Naive and b-lac dsRNA-injected cockroaches (Control), adult cockroaches at >7days after PameORco dsRNA injection (Adult injection), and adult cockroaches injected with PameORco dsRNA at the last instar stage (LI injection).
Single sw-A sensilla had two PameORco-positive OSNs (Figures 2A and 2E). Consistent with the previous studies, OSNs in sw-A sensilla in the Control group were broadly tuned to alcohols and terpenes ( Figure 5A; Sass, 1983;Fujimura et al., 1991;Tateishi et al., 2020). In the Control group, the response spectra to tested odors varied among different sw-A sensilla, suggesting that various ORxs are expressed in OSNs (Figure 5A). In both the Adult and LI injection groups, both the spontaneous activities and olfactory responses of OSNs in tested sw-A sensilla were significantly decreased (Figures 5A-5C). Consistent with the results of the sw-B sensilla, activities of OSNs were hardly observed in the LI injection group. These results indicated that OSNs in sw-A sensilla perceived general odors via functional ORs.
In single trichoid sensilla, a PameORco-positive OSN was consistently paired with a PameORco-negative OSN regardless of subtype; sw-C1 and shorter sw-C2 sensilla ( Figures 2C and 2E). SSRs from single sw-C2 sensilla revealed that there were OSNs with large and small spontaneous spikes, and OSN with small spikes consistently exhibited phasic-tonic excitatory responses to C6-acid and/or C7-acid RNAi cockroaches. Among ten tested odors, OSNs in the dw-A2 sensilla are selectively activated by three acids (See also Figure S4A). Responses to a given odorant are color-coded according to the response intensity. Olfactory responses of 6 sensilla from 5 Control cockroaches and of 13 sensilla from 10 PameORco RNAi cockroaches are summarized.
(E and F) Effects of PameORco RNAi on activities of dw-A2 sensilla. Both spontaneous activities (E) and responses to the best ligand (F: heptanoic acid; C7-acid) of dw-A2 sensilla were not affected by PameORco RNAi. Black bars in each plot indicate means. Throughout figures, symbols of the same color and shape correspond to recordings from the same individual. Significant differences are denoted above each plot (ANOVA post-hoc Tukey-Kramer test; *** = p < 0.001, n.s. = p > 0.05).
Effects of PameORco RNAi on sex pheromone processing in the central brain To evaluate the effect of PameORco RNAi on brain neurons, we focused on the PB-responsive PNs, which receive sensory inputs from a large number of PB-SNs in the B-glomerulus. Among 12-15 PB-responsive PNs in the B-glomerulus, L1-PN is morphologically homologous to other uniglomerular PNs that process ordinary odorants (Nishino et al., 2018). It sends axon terminals to specific regions of the mushroom body calyx and lateral horn via the medial antennal lobe tract (mALT) ( Figure 7A). We intracellularly recorded and stained L1-PNs and compared PB responses and dendritic arborizations in the B-glomerulus between Naive and PameORco RNAi cockroaches.
Silencing of the PB responses of PB-SNs neither affected the volume of the SN terminals ( Figure S5) nor dendritic arborizations of the L1-PN dendrites in PameORco RNAi cockroaches ( Figure 7B). It suggests silencing the Pa-meORco expression did not affect the development of both peripheral sensory neurons and central brain neurons. Although PB responses of PB-SNs were significantly impaired in Adult injection group (Figure 4), L1-PNs iScience Article exhibited moderately strong responses to PB, and response intensities and latencies of L1-PNs did not differ from those of Naive cockroaches ( Figure 7B-7D). Thus, our results suggest that brain neurons can compensate for the impaired olfactory responses of OSNs to some extent, and it is corresponding to the result of behavioral experiment using the Adult injection group cockroaches; they displayed in an interest in PB, but the ratio of individuals which exhibited the sexual behavior (wing-raising) is significantly lower than that of Naive cockroaches ( Figure S6).
In LI injection group, L1-PNs exhibited weak responses to PB and later response latencies ( Figures 7B-7D). The PameORco dsRNA injected at the LI stage effectively impairs olfactory responses in both peripheral and central neurons in the cockroach. iScience Article DISCUSSION In this study, we functionally analyzed the olfactory co-receptor in the cosmopolitan pest P. americana (Pa-meORco) and revealed fundamental mechanisms underlying the olfactory receptions in the insect. Here, whole sequence of PameORco and its sensillar type-specific expression patterns were characterized. Silencing of PameORco via systemic RNAi impaired receptions of sex pheromones and a large variety of general odors, especially alcohols and terpenes, through OSNs in the antennal perforated basiconic and trichoid sensilla. In addition, PameORco was essential for evoking the spontaneous activities of OSNs as has been reported in other insect species (Larsson et al., 2004;Stengl and Funk, 2013;Li et al., 2016;Butterwick et al., 2018). These results suggest that PameORco not only forms a functional heteromeric complex with the ligand-specific ORx but also forms a homotetramer on its own that controls the spontaneous activities of OSNs (Nakagawa et al., 2012;Butterwick et al., 2018;Wicher and Miazzi, 2021). In this study, we also revealed RNAi can effectively suppress the PameORco expression during long-lasting period beyond several molts. Our study is the first to characterize the fundamental molecular process of odor transduction in the cockroach, and provide new insights into further understanding of the molecular mechanisms of peripheral and central olfactory processing in the cockroach.
We characterized ORco from the genome of P. americana. Among species in the order Blattodea, whole alignment of ORco has been conducted with sequences from the German cockroach B. germanica and two species of termites Reticulitermes chinensis and Odontotermes formosanus (Robertson et al., 2018;Gao et al., 2020;He et al., 2020). Consistent with the result of a recent nuclear phylogenomic study of the order Blattodea (Evangelista et al., 2019), PameORco exhibits higher homology with ORco of termites than that of B. germanica ( Figure S1; Robertson et al., 2018). Because termites might have acquired eusociality by developing the chemical and pheromonal communication, analysis of P. americana olfactory receptors could provide new insights into the evolutionary transition from cockroaches to eusocial termites (Harrison et al., 2018). PameORco exhibits high homology with ORco in hemimetabolous insects but the low homology with ORco in holometabolous insects, such as ORco in D. melanogaster (DmelORco). However, the C-terminal region of PameORco is conserved across insect ORco, which suggests that this region might mediate the functional interactions of ORx and ORco proteins ( Figure 1B, Carraher et al., 2015). In addition, the putative phosphorylation and calmodulin-binding-sites-related channel activities that have been identified in DmelORco are conserved in PameORco (Kumar et al., 2013;Carraher et al., 2015;Guo et al., 2017;Jain et al., 2021). Thus, the odor transduction mechanisms for DmelORco may be similar to those for PameORco.
We established effective experimental condition for PameORco RNAi in the American cockroach. PameORco expression was almost completely suppressed at 3 days after PameORco dsRNA injection, but PameORco-expressing OSNs exhibited sufficient responses to ligand odors at that time point. The reduction of OSN responses plateaued at >7 days after PameORco dsRNA injection (Figure 4), suggesting that the turnover rate of OR expression in OSNs is approximately one week in the cockroach. The effect of PameORco RNAi was the most prominent in cockroaches injected with dsRNA at nymphal stage. The strong suppression of Pa-meORco expression was maintained for at least 88 days, and we did not observe recovery from the effects of RNAi in all tested cockroaches. Long-lasting and irreversible RNAi effects have also been observed when RNAi was applied to silencing genes related to visual and mechanosensory reception in the cockroach (French et al., 2015;Hennenfent et al., 2020). In other hemimetabolous insects, such as the cricket Gryllus bimaculatus and the firebrat Thermobia domestica, long-term effects of RNAi on circadian rhythm-related genes have also been reported (Uryu et al., 2013). These results suggest that RNAi may have long-lasting effects on hemimetabolous insects regardless of the gene targeted, and the RNAi technique utilize in any experiments.
The long-lasting and irreversible PameORco RNAi effects are very surprising because it has been believed that RNAi-mediated breakdown of target mRNA is transient, and the knockdown effect disappears after the depletion of short interfering RNA. In several insect species, mutagenesis of ORco by CRISPER/Cas9 induced the loss or smaller volumes of targeted glomeruli of OSNs expressing ORs (Trible et al., 2017;Fandino et al., 2019). However, our morphological analysis of a target glomerulus of PB-SNs indicated that ORco RNAi did not affect developmental process of OSNs ( Figure S5). Therefore, we hypothesize that the possible mechanism of the long-lasting effect is piwi-mediated transcriptional silencing (Onishi et al., 2021). Piwi is a paralog of Argonaute proteins, RNA-binding proteins required for RNAi machinery, and predominantly expressed in germline cells to suppress transposon activity. Interestingly, we found that the piwi homolog of P. americana is expressed in the antennae by conducting a blast search on the ll OPEN ACCESS iScience 25, 104272, May 20, 2022 iScience Article transcriptome database of the cockroach antennae (GenBank ID: SRR3089536 and SRR3089537; Chen et al., 2016). Further studies are needed to reveal the molecular/cellular mechanisms underlying dsRNAmediated long-lasting gene knockdown in long-living hemimetabolous insects.
Both anti-PameORco antiserum staining and SSRs combined with PameORco RNAi unambiguously revealed that the cockroaches detect the sex pheromones, PA and PB, via OSNs expressing ORx/ORco complex in sw-B sensilla (Figures 2 and 4). Transcriptome analysis has revealed that the two OR genes, PameOR1 and PameOR2, are predominantly expressed in adult male antennae and are thought to be ligand-specific ORxs of sex pheromone receptors . To determine the sex pheromone receptors in the cockroach, functional analysis of PameOR1 and PameOR2 using the methods established in this study is needed. Additional studies must reveal the sex pheromone receptors in the cockroach.
Sex pheromones in the cockroach elicit the sexual behavior in adult males but not nymphs (Tobin et al., 1981). However, nymphs also have both peripheral and central neurons to process sex pheromones (Schaller-Selzer, 1984;Nishino et al., 2009;Tateishi et al., 2020). In nymphal cockroaches, PB-SNs are housed in sw-A sensilla (Tateishi et al., 2020), and at final molt, the sex pheromone-responsive sw-A sensilla elongate and transform into sw-B sensilla (Schaller, 1978;Watanabe et al., 2018a). The results of a previous study and our study indicate that the sensitivities of PB-SNs to PB dramatically increased according to the sensillar elongation ( Figure 4G, Tateishi et al., 2020). The reported dose-response curve of PB-SNs in nymphal cockroaches is consistent with that in adult cockroaches at 3 days after PameORco dsRNA injection, whose ORs might be partially metabolized (Figures 3 and 4). Because the dendrites of PB-SNs expand in according to the sensillar elongation (Schaller, 1978), a greater number of receptors may be expressed on the dendrites at final molt. Thus, postembryonic changes of sex pheromone receptor expression may be linked with the sensitivities of sex pheromone-responsive OSNs, and trigger sexual behaviors in adult cockroaches.
Intracellular recordings from brain neurons using PameORco RNAi cockroaches revealed relationships between the activities of pre-and post-synaptic olfactory neurons. In adult cockroaches at 7 days after injection of Pa-meORco dsRNA, PB-SNs exhibited impaired responses to PB, whereas post-synaptic L1-PNs exhibited sufficient responses to 1 ng of PB, which were comparable to those in Naive cockroaches, and exhibited impaired behavioral responses to the PB. Responses of PB-SNs to 1 ng of PB in the cockroaches at 7 days after injection of Pa-meORco dsRNA were equivalent to the responses to 20 pg of PB in Naive cockroaches ( Figure 4G), suggesting that a low concentration of sex pheromone effectively activates brain neurons in Naive cockroach. In fact, 0.1 pg of PB can elicit sufficient excitatory responses in L1-PN of Naive cockroaches and elicited sexual behaviors (Boeckh and Selsam, 1984;Domae et al., 2019). The high sensitivity of the L1-PN is partly explained by the high sensory convergence of PB-SNs; approximately 37,000 PB-SNs send axons into the B-glomerulus and synapse onto a L1-PN (Boeckh and Ernst, 1987). Nevertheless, L1-PNs hardly responded to PB in the cockroach injected with PameORco dsRNA at the LI stage, suggesting that PameORco RNAi performed in the nymphal stage resulted in the near-complete suppression of sex pheromone processing.
The single-walled trichoid sensilla are common in many insect species, but their functions vary; these sensilla are specifically tuned sex pheromones in many moth species (Ng et al., 2020) and broadly tuned to a large variety of odors in D. melanogaster (Munch and Galizia, 2016) and the locust Locusta migratoria Yang et al., 2012;You et al., 2016). Regardless of the functional differences among insect species, OSNs in the trichoid sensilla only express ORs, with the exception of the cockroach as shown in this study. In the cockroach, trichoid sensilla are classified into sw-C1 and sw-C2 sensilla based on differences in their lengths and target glomeruli (Watanabe et al., 2012). In single sw-C2 sensilla, a Pa-meORco-positive OSN always pairs with a PameORco-negative OSN that is tuned to acids (Figures 2  and 6). Acids are specific ligands for IRs (Ai et al., 2010;Raji et al., 2019), suggesting that an OR-OSN and an IR-OSN are co-housed in a sensillum. Except for a subtype of double-walled coeloconic sensilla (ac4B sensilla) in D. melanogaster, single sensilla with OSNs expressing ORs and IRs have not been reported in other insects to date (Benton et al., 2009;Munch and Galizia, 2016). To confirm the specific features of the sensilla, the distributions and functions of IRs in the cockroach should be analyzed in the future. Although we were not able to record from sw-C1 sensilla in this study, single sw-C1 sensilla are equipped with two unique OSNs termed ON-cell and OFF-cell, which are excited by the onset and offset of a given odor stimulus, respectively (Burgstaller and Tichy, 2011;Tichy and Hellwig, 2018). Our anti-PameORco antiserum staining suggests that ON-OFF information of a given odor stimulus may be received by different types of olfactory receptors.

OPEN ACCESS
iScience 25, 104272, May 20, 2022 13 iScience Article PameORco-negative and acid-tuned OSNs in grooved basiconic sensilla strongly suggest that they also express IRs. In the cockroach brain, general olfactory signals detected by perforated basiconic sensilla are independently processed from those processed by trichoid and grooved basiconic sensilla through the two parallel olfactory pathways; OSNs in two different classes of sensilla project to different groups of glomeruli in the antennal lobe (AD glomeruli and PV glomeruli in Figure S5A) and synapse onto two different classes of PNs, which terminate in different areas of the mushroom body, and each area has dendrites of different classes of Kenyon cells (Watanabe et al., , 2012(Watanabe et al., , 2017Takahashi et al., 2019). Therefore, with the exception of PameORco-positive OSNs in trichoid sensilla, the OR-and IR-related olfactory circuits in the cockroach brain are ambiguously segregated. In the moth (Fandino et al., 2019) and fruit fly (Grabe and Sachse, 2018), odor information received by ORs and IRs elicits different odor-guided behaviors, suggesting that OR-and IR-related circuits have different functions. The simple but sophisticated olfactory processing from peripheral to higher brain centers in the cockroach is one of the best systems for studying how animals encode various odors and translate them into odor-guided behaviors via OR-and IRrelated parallel circuits.

Limitations of the study
Our study revealed the distribution of ORco and its function in odor reception in the antennae of the American cockroach, P. americana. ORco-impaired cockroaches exhibited partial anosmia, especially to sex pheromones, in a sensillum-specific manner. Electrophysiological experiment and RT-qPCR suggest systemic RNAi of PameORco chronically and irreversibly suppressed odor reception via odorant receptors (ORs). Because we can design dsRNA to target species-specific gene, RNAi is considered a robust approach for assessing the physiological actions of repellent or attractants in IPM (Gondhalekar et al., 2021). In this study, adult cockroaches at >7 days after PameORco dsRNA injection displayed in an interest in the sex pheromones, but the ratio of individuals which exhibited sexual behaviors is significantly lower than that of Naive cockroaches. It is highly corresponding to the results of phylogenetically related German cockroach; behavioral responses to the sex pheromone and food odors were partially disrupted in adult cockroaches at 3 days after ORco dsRNA injection (He et al., 2020). Our results in the American cockroach suggest that the injection of a small amount of PameORco dsRNA to nymphal cockroaches more effectively impairs odor and sex pheromone receptions in adults, and the effect is long-lasting and irreversible. It needs the further behavioral studies using the adult cockroaches which received dsRNA injection at nymphal stage. In addition, it needed to determine the most efficient way to deliver dsRNA to nymphal cockroaches (e.g., baits including the dsRNA) under field conditions.

STAR+METHODS
Detailed methods are provided in the online version of this paper and include the following:

Lead contact
Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Hidehiro Watanabe, Ph. D (nabehide@fukuoka-u.ac.jp).

Material availability
Antiserum against to PameORco generated in this study is available from the lead contact upon request.

Data and code availability
Section 1: All raw data reported in this paper will be shared by the lead contact upon request.
Section 2: This paper does not report original code.
Section 3: Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

EXPERIMENTAL MODEL AND SUBJECT DETAILS Animals
Tenth instar, last instar (LI: 11th instar), and adult male cockroaches P. americana with intact antennae were obtained from laboratory colonies maintained at 27 G 1 C under a 12:12 light-dark cycle at Fukuoka University. The nymphal stages were identified according to the lengths of the body and hind tibia (Gier, 1947).

REAGENT or RESOURCE SOURCE IDENTIFIER
Primer with T7 promoter sequence (underlined) used to amplify the template DNAs for in vitro transcription for E.coli b-lac,

Molecular cloning of PameORco
We cloned an ORco homolog expressed in the cockroach. First, we conducted a tblastn search in the genomic DNA sequence database of P. americana (GenBank ID: PGRX00000000.1) using the deduced amino acid sequence of Drosophila melanogaster ORco protein (GenBank ID: Q9VNB5). Next, we performed RT-PCR to amplify cDNA fragments containing the full-length open reading frame (ORF) of the gene using primers designed at the predicted 5 0 and -3 00 -untranslated regions of PameORco. Complementary DNA (cDNA) was prepared according to Watanabe et al. (2010b). Nucleotide sequences of the primers are listed in in key resource table. The nucleotide sequences of the obtained cDNAs of PameORco isoform 1 and 2 were registered in GenBank (GenBank ID: LC657818 and LC657819).

Structural and phylogenetic analyses of insect ORco homologs
The deduced amino acid sequence of PameORco isoform 2 was aligned with those of ORco homologs of other insect species using the MAFFT algorithm (Katoh and Standley, 2013) and refined by manual inspection using the Geneious Prime program (v2021.

Immunohistochemistry
The polyclonal anti-PameORco antiserum was generated in a guinea pig against the 14-a.a. synthetic peptide (NH 2 -C+VSANTNNPAATGDD-COOH) located within the second intracellular loop of the protein (Figure 1; residues 263 to 276) (Eurofins genomics, Tokyo, Japan). We also used a polyclonal antibody generated in rabbit against peroxidase from horseradish (anti-HRP; Sigma-Aldrich Cat# P7899, RRID: AB_261181). We previously showed that the anti-HRP antibody selectively labeled sensory neurons in cockroach antennae .
Immunohistochemistry protocols followed those used for the cockroach antennae with modifications (Watanabe et al., 2018a). The two antennae of a day-old adult cockroaches were removed, and each antenna was divided into several fragments. Each fragment was manually sliced along the longitudinal axis with a razor blade. Sliced fragments were fixed in 4% paraformaldehyde with 0.1 M phosphate buffer (PB, pH 7.4) at 4 C for a day. After washing in PB saline (PBS: 0.14 M NaCl in 0.01 M PB, pH 7.4) with 1% Triton X-100 (PBST), the fragments were blocked for an hour using 5% normal goat serum (NGS; 143-06561, Fujifilm Wako, Osaka, Japan) in PBST with 0.25% bovine serum albumin (PBST-BSA; A2153, Sigma-Aldrich, St. Louis, MO). Antennal fragments were then incubated with a mixture solution of polyclonal anti-HRP antibody (rabbit) and anti-PameORco antiserum (guinea pig) at 4 C for 4 days. The anti-HRP antibody and anti-PameORco antiserum were diluted 1:5000 and 1:1000 in PBST-BSA with 2% NGS, respectively. On the following day, the fragments were washed again in PBST-BSA and then incubated with the mixture solution of Alexa488-labeled secondary anti-rabbit IgG antibody (diluted 1:1000; A32731, Thermo Fisher Scientific; RRID: AB_2633280) and Alexa 555-labeled secondary anti-guinea pig IgG antibody (diluted 1:1000; A21435, Thermo Fisher Scientific; RRID: AB_2535856) in PBST-BSA with 2% NGS at 4 C for 4 days. The fragments were then finally washed in PBST-BSA, dehydrated with an ascending ethanol series (from 50% to 100%), and cleared in methyl salicylate. No signal was detected when specimens were labeled with only secondary antibodies (data not shown).

Preparation and injection of double-stranded RNA (dsRNA)
DsRNA for PameORco and Escherichia coli b-lactamase (b-lac) genes were synthesized by T7 RNA polymerase-based in vitro transcription. A partial cDNA fragment of each target gene with T7 promoter sequences at both ends was used as a template for in vitro transcription. The details for each template are provided below.

OPEN ACCESS
iScience Article PameORco A 424-bp partial cDNA fragment of PameORco encoding the C-terminal region (including the sixth and seventh transmembrane helixes) of the protein was used as a template for in vitro transcription. The nucleotide sequence of PameORco dsRNA is specific for P. americana, and the PameORco dsRNA targets both mRNA isoforms. The cDNA was amplified by PCR from a plasmid containing a partial cDNA clone of Pa-meORco. PCR was carried out with a pair of gene-specific primers (GSPs) with the T7 promoter sequence at the 5 0 end (Key resource table).

E. coli b-lac
A 369-bp partial cDNA fragment of E. coli b-lac was used as a template for in vitro transcription. The cDNA fragment was amplified by PCR from the pGEM-T easy vector (Promega, Madison, WI, USA) with a pair of GSPs with the T7 promoter sequence at the 5 0 end (Key resource table).
All PCRs were carried out using the Q5 High-Fidelity DNA polymerase (New England Biolabs, Tokyo, Japan), and the PCR products were purified using the QIAquick PCR Purification Kit (Qiagen, Tokyo, Japan). Nucleotide sequences of the primers were listed in key resource table. In vitro transcription was carried out using the HiScribe T7 High Yield RNA Synthesis Kit (New England Biolabs, Tokyo, Japan) per the manufacturer's instructions. Approximately 1 mg of template DNA was added to a 20 mL reaction, and incubated at 37 C for $16 hr. The dsRNAs were purified by lithium chloride precipitation. The synthesized dsRNAs were diluted into the nuclease-free phosphate-buffered saline to adjusted the dsRNA concentration to 0.8 mg/mL and were stored at À30 C.
For systemic RNAi of PameORco, we prepared 4 days-old adult males or 7 days-old nymphs ( Figure 3A). After anesthetizing the cockroaches on ice, 5 mL of dsRNA solution (4 mg of dsRNA in total) was manually injected into the dorsal region of the head capsule using a 27 G needle attached to a 10-mL glass syringe. The cockroaches injected with dsRNA were separately reared in isolated glass tubes, and provided food and water ad libitum. Adult cockroaches were used for experiments at 3, 7 or 10 days after dsRNA injection. The nymphal cockroaches injected with dsRNA were reared until they became adults and used at 11 days after their final molt ( Figure 3A). In control experiments, 11 days-old virgin males were used as Naive cockroaches.

Quantitative reverse transcription PCR (RT-qPCR)
The RT-qPCR analysis was conducted following the methods of Watanabe et al. (2018b) with modifications. A pair of antennae was homogenized in 500 mL of the TRIzol Reagent (Life Technologies, CA, USA) using a beads homogenizer (mT-01, TAITEC, Saitama, Japan) at 4600 r/min for 90 s at room temperature. The total RNA was extracted from each sample according to the manufacturer's instructions. After RNA extraction, genomic DNA was digested with Recombinant DNase I (5 U/reaction; TaKaRa, Shiga, Japan). After heatinactivation of the enzyme (75 C, 10 min), each RNA sample was reverse transcribed using the High-Capacity cDNA Reverse Transcription Kit (Life Technologies) according to the manufacturer's instruction. A 1.5 mL of total RNA solution was added to the 10 mL reaction.
The relative quantification of target genes was carried out using the KAPA SYBR Fast qPCR Kit (Kapa Biosystems, MA, USA) and the Thermal Cycler Dice Real Time System (TaKaRa, Shiga, Japan). Gene expression levels were measured using the standard curve method. The expression levels of PameORco mRNA were normalized to that of PameEf1a (GenBank ID: LC657820) using the 2 ÀDDCt method. The nucleotide sequences of the primers are listed in key resource table.
Single sensillum recording from OSNs in single olfactory sensilla An ice-anesthetized cockroach was immobilized ventral-side-up on an acrylic plate. To prevent movement, the body and legs were gently mounted using low-melting wax, the neck was immobilized with small acrylic plates, and the antennae were fixed using the wax. The antennae were observed through the light microscope at 5003 magnification (AZ100, Nikon, Tokyo, Japan).
We performed single sensillum recordings (SSRs) from arbitrarily selected olfactory sensilla on the ventral surface of the flagellum.  , 1991;Tateishi et al., 2020) and two purified synthesized sex pheromones (PA and PB; Mori, 1990a, 1990b). Both PA and PB were diluted in hexane at a concentration of 0.1 ng/mL as stock solutions. Each sex pheromone solution was placed on an aluminum plate (15 3 5 mm) and was dried for more than 1 min to evaporate the solvent. Therefore, their concentrations were denoted as dry weights. General odors were diluted in paraffin oil at a concentration of 10 mM, and 20 mL of the odorant solution was added to a piece of filter paper (15 3 5 mm). Immediately before recordings, the aluminum plates and the filter papers loaded with odorants were separately inserted into glass pipettes.
Fresh air was cleaned and dried with charcoal and silica-gel filters. The air stream was maintained at 1 L/min using a flowmeter. The main tube was tandemly connected to two electric-driven three-way solenoid valves ( Figure S3C; Valve 1 and 2) which were independently operated by Powerlab. During the inter-stimulus period, the constant air which passed through the two valves and the blank glass pipette flowed over the antenna (top panel of Figure S3C). During the general odor stimulation period, the constant air stream was stopped and the air stream from the other outlet of Valve 1 passed through the glass pipette containing a given odorant. The tip of the glass pipette was positioned approximately 1 cm apart from the recording sensillum. The general odor stimulus was presented for 0.5-s period and each odorant was presented two to five times with > 40-s intervals. After stimulation, a new glass pipette containing another odorant was attached to the outlet of Valve 1.
To discriminate the PA response of PA-SN and the PB response of PB-SN, cross-adaptation experiments were performed. A glass pipette containing PB and one containing PA were connected to the outlets of Valve 1 and Valve 2, respectively. In the experiment, PA was presented at a concentration equal to PB. The tips of two glass pipettes were positioned approximately 1 cm apart from the recording sensillum. During the sex pheromone stimuli period, the recording sensillum received three successive PB stimuli and then one PA stimulus with inter-stimulus intervals of 0.5 s ( Figure S3D). Each sex pheromone stimulus was presented for 0.5-s periods. The cross-adaptation experiments were performed with > 2-min intervals. PA-SN and PB-SN easily adapt to PA and PB, respectively. Although low to high concentrations of sex pheromones were presented, each concentration was presented to the recording sensillum only once during each recording period. The air around the preparation was continuously exhausted through a duct behind the recording electrode. We regarded the timing of the solenoid valve gating as the onset of odor stimulation.

Data analysis of SSRs
In each SSR, different shapes of spikes from different OSNs were concurrently recorded. Although the spike amplitudes were distinguishable in several traces, the spikes were not easily distinguishable in many SSRs. We grouped spikes with different shapes in each SSR and calculated the increase in total spike frequency ll OPEN ACCESS iScience 25, 104272, May 20, 2022 iScience Article from the spontaneous frequency as follows; RÀR0, where R0 and R were the total numbers of spikes during the 0.5-s period before and after the onset of odor stimulation, respectively. The average response intensity was defined as the response intensity to a given odorant. Spontaneous spikes were counted during the 10-s period (in perforated basiconic and grooved basiconic sensilla) or 40-s period (in trichoid sensilla) before the onset of all given odor stimuli, and the spontaneous firing rate was determined by converting the averaged numbers to Hz.
Each of spikes derived from OSNs in single sw-B sensilla and sw-C2 sensilla was sorted using the spike sorting function of Spike 2 ver. 8.08 (CED, Cambridge, UK) (Figures 4, 6, and S3). We extracted each of spikes with the duration of 2 msec-period based on its spike peak. And then, principal component analysis was performed in all recorded spikes according to the spike shapes. Each spike was plotted in a three-dimensional space using the first three principal components (PC1 -PC3), and clustered ( Figure S3B). Because spikes from PB-SN and PA-SN was difficult to segregate by their shapes, spikes pre-adapted to PB stimulus and those post-adapted to PA stimulus were regarded as PB-elicited spikes of PB-SN and PA-elicited spikes of PA-SN, respectively.
In electrophysiological experiments, statistical differences were determined by one-way ANOVAs and post-hoc Tukey-Kramer test using BellCurve for Excel (Social Survey Research Information, Tokyo, Japan).

Intracellular recording and staining of L1-PN
A glass microelectrode was filled with 8% Lucifer Yellow (Sigma-Aldrich) in 1 M LiCl 2 . The input resistances were set 30-50 MU. A cockroach, briefly anesthetized by carbon dioxide, was fixed onto a handmade acrylic chamber with wax. The head was positioned in a shallow chamber and fixed with low melting point wax. A small rectangular window was cut in the head cuticle between the compound eyes, muscles and trachea were removed to expose the brain surface. The brain was immersed in cockroach saline to form an electrical bridge with an indifferent silver rod (diameter, 300 mm) that was inserted into a cavity between the circumesophageal connectives. For intracellular recording and subsequent staining, the electrode was inserted into the axon just below the point where two calyces meet ( Figure 7A). Electrical signals were amplified by a DC amplifier (MEX-8301, Nihon Kohden, Japan), displayed in an oscilloscope, and then digitized. L1-PNs were readily identifiable based on their spontaneous activities and responsiveness to PB. The 1 ng of PB stimulus was presented to antennae for 0.3-s periods 5-8 times with intervals of >30 s. Statistical differences were determined by one-way ANOVAs and post-hoc Tukey-Kramer test.
After recording olfactory responses, the L1-PN was filled with Lucifer Yellow by injecting a hyperpolarizing current. After staining, the electrode was removed from the brain, and the head capsule was fixed on a wax plate. Anterograde staining of antennal afferents was then conducted as described the above section. The double-stained brain was fixed, dehydrated, and cleared.

LSM observations and three-dimensional reconstruction of the B-glomerulus
The cleared specimens were examined with a confocal laser scanning microscope (LSM-510; Carl Zeiss, Jena, Germany) equipped with Argon and Helium-Neon lasers. Neurons labeled by Alexa488 or Lucifer Yellow were visualized using an Argon laser with a band-pass emission filter (505-530 nm), whereas those labeled Alexa555 or micro-ruby were visualized using a Helium-Neon laser with a long-pass emission filter (>560 nm). Antennal sensilla were visualized via the differential interference contrast (DIC) imaging function of the LSM. Images were obtained using three different objectives: a Plan-Apochromat 203/0.8 objectives for low magnification images, and an oil-immersion Plan-Neofluar 403/1.3 and Plan-Apochromat 633/1.4 objectives were used for high magnification images. Optical sections were captured at intervals of 0.7-1.2 mm thickness with 1024 3 1024 pixels. The contrast and brightness of all images were adjusted using Adobe Photoshop CS3 and Illustrator CS3. The orientations and positions of the brain are shown relative to the body axis.
For the three-dimensional reconstruction of the B-glomerulus, a series of TIFF-formatted optical sections of macroglomeruli were processed using image processing software ( Figure S5B; Amira 6.0, Thermo Fisher Scientific). The boundary of B-glomerulus was manually traced with the aid of the anterogradely stained OSN afferents. The volume of the 3-D reconstructed B-glomerulus was calculated by the function attached in Amira 6.0. Volumes of B-glomeruli in three RNAi groups were statistically compared by Kruskal-wallis test. We compared sexual behaviors to PB between Naive and PameORco RNAi cockroaches. We collected 4 days-old adult virgin males and divided them into Naive and PameORco RNAi (Adult injection) groups.
In the Naive group, two or three cockroaches were reared together in a plastic cup (4 = 8 cm (bottom), height = 8 cm) spread a filter paper, and provided a paper shelter, food and water for 7 days. In the Adult injection group, two or three cockroaches injected 4 mg of PameORco dsRNA were reared together in a plastic cup for 7days. After then, we removed the paper shelter, food and water from the plastic cup, and performed behavioral experiment at three hours later. We observed cockroach behaviors to PB using an infrared camera (FDR-AX60, Sony, Tokyo, Japan) under the dark situation. At 4 min after the onset of recording, 1 ng of PB on a filter paper (15 3 5 mm) was gently placed on the plastic cup, and recorded behaviors to PB during 4-min period (Videos S5 and S6). After recording, tested cockroaches were again placed a new plastic cup, and were received the same behavioral experiment at 7 days later.
We evaluated the behavior of individual cockroach during the 2-min period after presentation of PB. We identified levels of sexual behavior of P. americana according to the previous studies as follows: 0: no response, 1: waving of antennae, 2: waving of antennae and oriented locomotion, 3: waving of antennae, oriented locomotion, partial wing-raising and wing-fluttering response, 4: waving of antennae, oriented locomotion, wing-raising, abdominal extension ( Figure 6S; Silverman 1977;Yang et al., 1998). We statistically compared the ratio of cockroaches which exhibit wing-raising behavior (levels 3 and 4) between Naive and PameORco RNAi groups using the Fisher's exact test.

QUANTIFICATION AND STATISTICAL ANALYSIS
To validate the effectiveness of PameORco RNAi, we statistically compared experimental data obtained from individuals at different timepoints after dsRNA injection. We used the one-way ANOVAs and posthoc Tukey-Kramer test to assess the results of RT-qPCR and electrophysiological experiments, the Kruskal-wallis test for morphological analysis of B-glomerulus, and the Fisher's exact test for behavioral experiments. Statistical significance was calculated using BellCurve for Microsoft Excel (Social Survey Research Information, Tokyo, Japan). Black bars within the scatterplots indicate means. Error bars in Figures 4G and 4H, and shaded areas in Figure 7D represent as mean G SEM The asterisks indicate significant differences between groups (n.s > P = 0.05, * = P < 0.05, ** = P < 0.01, *** = P <0.001). The different letters above each of plots indicate significant differences according to multiple comparisons with Tukey-Kramer test (P < 0.05).