- Research article
- Open access
- Published:
Pheromone sensing in Drosophila requires support cell-expressed Osiris 8
BMC Biology volume 20, Article number: 230 (2022)
Abstract
Background
The nose of most animals comprises multiple sensory subsystems, which are defined by the expression of different olfactory receptor families. Drosophila melanogaster antennae contain two morphologically and functionally distinct subsystems that express odorant receptors (Ors) or ionotropic receptors (Irs). Although these receptors have been thoroughly characterized in this species, the subsystem-specific expression and roles of other genes are much less well-understood.
Results
Here we generate subsystem-specific transcriptomic datasets to identify hundreds of genes, encoding diverse protein classes, that are selectively enriched in either Or or Ir subsystems. Using single-cell antennal transcriptomic data and RNA in situ hybridization, we find that most neuronal genes—other than sensory receptor genes—are broadly expressed within the subsystems. By contrast, we identify many non-neuronal genes that exhibit highly selective expression, revealing substantial molecular heterogeneity in the non-neuronal cellular components of the olfactory subsystems. We characterize one Or subsystem-specific non-neuronal molecule, Osiris 8 (Osi8), a conserved member of a large, insect-specific family of transmembrane proteins. Osi8 is expressed in the membranes of tormogen support cells of pheromone-sensing trichoid sensilla. Loss of Osi8 does not have obvious impact on trichoid sensillar development or basal neuronal activity, but abolishes high sensitivity responses to pheromone ligands.
Conclusions
This work identifies a new protein required for insect pheromone detection, emphasizes the importance of support cells in neuronal sensory functions, and provides a resource for future characterization of other olfactory subsystem-specific genes.
Background
To fulfil the formidable task of detecting and discriminating many diverse chemical signals in the external world, animal olfactory systems contain tens to thousands of distinct olfactory sensory neuron (OSN) populations. In both vertebrates and insects, each neuronal population is distinguished by the expression of a specific olfactory receptor (or, occasionally, receptors) and their innervation pattern in the brain [1, 2]. In most species, olfactory systems are also characterized by the gross-level organization of different sets of OSNs into structurally and functionally distinct “subsystems.” In rodents, for example, the main olfactory epithelium and vomeronasal organ express different receptor families and have distinct (though not exclusive) roles in detecting environmental odors and pheromones, respectively [3,4,5,6].
Insects also possess olfactory subsystems, which have been best-described in Drosophila melanogaster [7,8,9,10]. The main D. melanogaster olfactory organ is the antenna (Fig. 1A), a head appendage covered with several hundred sensory sensilla [11, 12]. Each sensillum constitutes a porous cuticular hair housing the ciliated dendrites of 1-4 OSNs, whose somas are flanked by non-neuronal support cells [13]. Antennal OSNs can be categorized into two subsystems, defined by their expression of members of distinct olfactory receptor repertoires: the odorant receptors (Ors) [14,15,16,17,18] and the ionotropic receptors (Irs) [10, 19]. Receptors from both families function as odor-gated ion channels, composed of ligand-specific “tuning” receptors and one or more broadly expressed co-receptors: Orco for Ors; Ir8a or Ir25a and Ir76b for Irs [20,21,22,23,24,25,26,27]. Although the Or and Ir subsystems share a similar overall organization and are intermingled within the antenna, they display a number of important differences.
During development, the Or and Ir subsystems are specified by distinct transcription factors—absent multidendritic neurons and olfactory sensilla (Amos) and Atonal (Ato), respectively—which determine the fate of lineages derived from sensory organ precursors in the developing larval antennal imaginal disk [33,34,35]. Amos and Ato induce the expression of downstream genes required to form the corresponding olfactory subsystem, encompassing neurons, non-neuronal support cells, and other structural elements. Morphologically, the sensilla in the Or subsystem comprise three main classes (basiconic, trichoid, and intermediate) while Ir subsystem OSNs are housed in coeloconic sensilla [11, 12]. Some Ir neurons are found in hygrosensory or thermosensory sensilla in the sacculus (an internal multi-chambered pocket) and arista (a long cuticular projection) (Fig. 1A) [30, 36,37,38,39]. Different sensillar classes exhibit several distinct structural features, including hair length and diameter, pore size and number, and the sensory cilia morphology of OSNs and the number of support cells [12, 40,41,42].
The response properties of Or and Ir OSNs are also distinctive, encompassing different odor specificities, sensitivities, and temporal dynamics [10, 43, 44]. Much effort has focused on addressing how Ors or Irs define odor-evoked activity, demonstrating that OSN signaling properties are largely determined by the corresponding tuning receptor [26, 45,46,47,48,49]. However, support cells are also likely to have several important contributions to olfactory detection. During development, the non-neuronal cells have roles in secreting and shaping the cuticular structure of olfactory hairs [13, 50]. In mature antennae, support cells are thought to form an isolated biochemical microenvironment for the OSNs within a given sensillum, through their control of the ionic composition of lymph fluid and the secretion of odorant-degrading enzymes and odorant-binding proteins (Obps) [13, 51]. Obps have various perireceptor roles in modulating olfactory responses [13, 52,53,54,55]. For example, the Obp Lush is secreted by Or subsystem trichoid sensilla support cells and is critical for the responses of Or67d neurons to their cognate pheromone ligand [56, 57]. The differential expression of this and other Obps in the Or and Ir subsystems [54] raises the question of the existence of other types of olfactory subsystem-specific support cell proteins that contribute to olfactory detection properties.
Results
A transcriptomic screen for olfactory subsystem-specific genes
To identify subsystem-specific molecules, encompassing those expressed in neuronal and/or non-neuronal cells, we performed comparative transcriptomics of antennae from animals mutant for either amos or ato, which selectively lack the Or or Ir subsystems, respectively (Fig. 1A) [33,34,35, 58]. This strategy was designed to facilitate sensitive identification of differentially expressed genes between the two subsystems, complementary to previously described independent comparisons of the antennal transcriptomes of these mutants against controls [59, 60]. We isolated RNA from amos or ato antennal olfactory segments (in three biological replicates) and hybridized these samples to D. melanogaster microarrays (Fig. 1A).
Comparison of transcript levels for positive control genes (i.e., Ors and Irs) validated the selectivity and sensitivity of this screen (Fig. 1B and Additional file 4: Table S1): with very rare exceptions, all transcripts for antennal Ors (but not those expressed in other olfactory organs [11], or the exceptional Ir subsystem-expressed Or35a [47]) were detected at significantly lower levels in amos compared to ato antennae (Fig. 1B and Additional file 4: Table S1). Conversely, transcripts for essentially all antennal Irs (but not family members expressed in non-olfactory organs [31, 61,62,63]) were expressed at lower levels in ato antennae compared to amos mutants. We further analyzed members of the “gustatory receptor” (Gr) gene family: most Grs are expressed in contact chemosensory organs [64] and, consistently, are not differentially expressed in these olfactory subsystems (Fig. 1B). However, a few Grs have sensory functions in the antenna: Gr21a and Gr63a encode subunits of a carbon dioxide receptor expressed in basiconic sensilla [65, 66] and their transcripts were enriched, as expected, in the Or subsystem (Fig. 1B). By contrast, transcripts for Gr28b.d (encoding an aristal-expressed thermoreceptor [38]) were expressed predominantly in the Ir subsystem (Fig. 1B). Several other Grs were mildly enriched in one or other subsystem but their endogenous expression and function (if any) is unclear (Fig. 1B) [67].
To characterize other differentially expressed genes of the olfactory subsystems, we first investigated whether those displaying >2-fold higher expression in either subsystem were enriched for particular Gene Ontology (GO) terms (Additional file 9: Data S1) [68]. Almost all significantly enriched terms reflected the differential expression of the olfactory receptor gene families (e.g., Or subsystem: GO:0007608 sensory perception of smell (FDR q-value 1.11E−13); Ir subsystem: GO:0004970 ionotropic glutamate receptor activity (FDR q-value 6.24E−06)), suggesting that there are few other broad categories of genes that distinguish these subsystems. We did note that the Or subsystem was enriched for genes encoding proteins with glutathione transferase activity (GO:0004364, FDR q-value 1.86E−07). Such enzymes have putative roles in odor degradation [69], and this enrichment might reflect differences in the types of chemical ligands detected by each subsystem [10].
We subsequently focused on genes displaying an expression difference between amos and ato antennae of >4-fold. This threshold captured the most enriched 177 and 180 genes of the Or and Ir subsystems, respectively—including all of the known chemosensory receptors—and facilitated manual curation and prioritization (Additional file 4: Table S1). Through protein domain and BLAST analyses, we assigned these genes to distinct categories (Fig. 1C). These classes include small, secreted proteins, notably Obps, such as the Or subsystem-expressed Lush [56] and the Ir subsystem-enriched Obp59a, which functions in hygrosensation in the sacculus [54, 70]. Many differentially expressed genes encode enzymes, including those potentially involved in odorant degradation and intracellular metabolism. A number of genes encode ion channels or transporters, some of which have known roles in olfactory signaling. For example, the Ir subsystem-specific Ammonium transporter (Amt) is a non-canonical olfactory receptor underlying ammonia detection in coeloconic and sacculus neurons [60, 71], while the Or subsystem is enriched for the phospholipid flippase/transporter ATPase 8B, which is required for olfactory sensitivity in several classes of basiconic and trichoid OSNs [72, 73]. Several genes encode proteins with probable roles in sensillum construction and/or maintenance, encompassing OSN cilia morphogenesis (e.g., CG45105, an ortholog of human SDCCAG8, which is implicated in ciliopathies [74]) and cuticle formation (e.g., Vajk3 [75]); some of these may contribute to the distinctive ultrastructural properties of different sensillar classes [12, 40, 41]. Differentially expressed transcription factor genes include unplugged, which is required for correct receptor expression and axon targeting in several Or OSN populations [76]. Diverse additional functional categories were represented, including post-translational gene regulation, neural guidance, immune signaling, and proteolysis, while many genes encode proteins with no obvious known domains. Notably, more than half of the differentially expressed genes are functionally uncharacterized, either in or beyond the olfactory system.
Diverse cell type and breadth of expression of novel olfactory subsystem-enriched genes
To further characterize these differentially expressed genes, we first took advantage of the single-cell RNA-sequencing (scRNA-seq) dataset of the antenna, generated as part of the Fly Cell Atlas [29], to obtain insights into their cellular-level expression patterns. This dataset encompasses the nuclear transcriptomes of cells from all segments of the antenna, which have been clustered and annotated into several dozen neuronal and non-neuronal cell types ([29] and Fig. 1D). We first asked whether the olfactory subsystem-enriched genes were predominantly expressed in neurons or non-neuronal cells by inspection of their expression pattern across cell clusters. We used N-synaptobrevin (N-syb) as a marker of all neuronal cell types (Fig. 1D) and found that other subsystem-enriched genes (excluding chemosensory receptors) could be neuronal, non-neuronal, or both (Fig. 1D and Additional file 4: Table S1). In both Or and Ir subsystems, the majority of detected genes was expressed in non-neuronal cells.
When we used these scRNA-seq data to survey the breadth of expression of each gene, we found patterns ranging from a single cluster to very broad expression across many cell classes of the Ir or Or subsystems (Additional file 4: Table S1). Notably, the vast majority of neuronal genes are broadly expressed, in contrast to the Ors and Irs which are restricted to a single cluster [29]. We confirmed these transcriptomic data by RNA fluorescence in situ hybridization (FISH) for the Or subsystem-enriched Glutathione S transferase E4 (GstE4) and the Ir subsystem-enriched Tetraspanin 47F (Tsp47F), which are co-expressed with the Orco and Ir8a co-receptors, respectively (Fig. 1E). Although we cannot exclude the possibility that lowly expressed neuron subtype-specific genes are not represented in the scRNA-seq datasets, our data reinforce previous conclusions that the olfactory receptor is the main, and perhaps only, determinant of tuning specificity of most individual neuron classes [26, 45,46,47,48,49]. Other neuronal, subsystem-specific genes identified here may have broader functions in defining these neurons’ signaling, metabolic, and/or morphological properties.
We next focused on the non-neuronal genes enriched in either olfactory subsystem. Many non-neuronal cells in the antenna develop independently of the cell lineages of the Or and Ir subsystems, including epithelial cells, muscle cells (which are probably located only within the proximal, non-olfactory antennal segments [77, 78]), and a subset of glia (which migrate from the central brain during development [79]) (Fig. 2A). A distinct set of antennal glia derive from the ato-dependent lineages [79], although it is likely that they contribute to the function of both olfactory subsystems. Non-neuronal cells of the olfactory subsystems—known as support (or auxiliary) cells—derive from the same developmental lineages that produce the OSNs within a given sensillum and comprise several classes: trichogen, thecogen, tormogen (described further below). These cell types are represented by ~20 cell clusters in the antennal atlas, globally demarcated by expression of the transcription factor gene shaven (sv) (Fig. 2A) [29]. The olfactory subsystem-specific genes that we identified display diverse breadth of expression in these cell types (Additional file 4: Table S1), several of which could be validated through RNA FISH (Fig. 2B–E). Some are expressed in multiple cell clusters, such as Juvenile hormone esterase duplication (Jhedup) in the Or subsystem (Fig. 2B) and CG10357, a5 and a10 (all of unknown function) in the Ir subsystem (Fig. 2C). Others are prominently expressed in only one cluster, including the Or subsystem-specific Osiris 8 (Osi8) (Fig. 2B) and the Ir subsystem-specific CG14153, CG13285, and CG34456, which are found in support cells around the sacculus (Fig. 2D). Beyond support cells, we detected glial expression for the Ir subsystem gene defective proventriculus (dve) (Fig. 2E).
We noted that the sacculus-expressed genes were detectable in the same cluster in the Fly Cell Atlas—which was also marked by the previously characterized Obp59a [54, 70] (Fig. 2D)—allowing us to annotate this cluster (arbitrarily named “I11”; Additional file 5: Table S2) as “sacculus support cell” (Fig. 2F,G). We extended this approach by systematically noting the support cell clusters expressing Or and Ir subsystem-enriched genes (Additional file 5: Table S2), which enabled demarcation of essentially all clusters as corresponding to one or other olfactory subsystem (Fig. 2F,G). We were unable to further define support cell types—although I12 and I13 are likely to be thecogen cells (Additional file 5: Table S2)—as many molecular markers are expressed in multiple cell clusters (Fig. 2F,G and Additional file 5: Table S2) and/or multiple sensillum classes (e.g., Obps [54]).
Expression of Osi8 in trichoid sensilla tormogen support cells
We subsequently focused on characterizing the Or subsystem gene Osi8, a member of a family of 24 Osi genes that each encode an N-terminal signal peptide, a domain of unknown function (DUF1676), followed by a presumed transmembrane domain (Fig. 3A). This family appears to be largely insect-specific, with syntenically arranged orthologs of most members present in diverse insect orders and, with rare exceptions [84], no related genes in the genomes of non-insect Arthropoda or other invertebrates [84,85,86]. Notably, of the analyzed insect species, Osi8 is absent in only the human body louse Pediculus humanus [86], which also has a greatly reduced Or repertoire compared to many other insects [87].
Osi8 was the only subsystem-enriched member of this family (Additional file 4: Table S1). To confirm and extend this observation, we quantified all Osi gene transcripts in a bulk RNA-seq dataset of adult antennae from control and amos mutant animals (Fig. 3B). Osi8 is expressed at >10-fold the level of any other Osi gene. Moreover, consistent with our microarray analysis, Osi8 transcripts are almost completely lost in amos mutants (Fig. 3B). Using the Fly Atlas 2.0 [88], we further examined Osi8 expression in RNA-seq datasets across diverse D. melanogaster tissues of both adults and larvae, and found the highest expression levels in adult heads of both males and females (Fig. 3C). This signal is most likely due to the antennal expression of Osi8 because, within the Fly Cell Atlas whole head dataset [29], transcripts of this gene were detected in only a very small subset of sv-positive cells, which presumably correspond to the antennal support cell population (Fig. 3D). Finally, analysis of tissue-specific RNA-seq datasets in the mosquito Aedes aegypti [89] revealed high expression of the Osi8 ortholog (AAEL004275) in the antenna (Fig. 3E). However, transcripts could also be detected in the mosquito brain (Fig. 3E), suggesting additional roles for the gene in this species.
Analysis of the spatial distribution of Osi8 RNA FISH signal revealed broad expression around the latero-distal region of the antenna, where trichoid and intermediate sensilla are located (Figs. 2B and 3F; specificity of signal is validated in Fig. 4B), but not in the region where basiconic sensilla are most abundant [11, 12]. Closer examination revealed two types of Osi8-expressing cells: one population is distally located in the antenna and has an elongated cell morphology; the other is more proximal, exhibiting a rounder shape and has, generally, a weaker Osi8 RNA FISH signal (Fig. 3F). To examine whether Osi8 expression is associated with specific sensillar classes, we visualized Osi8 RNA simultaneously with GFP transgenic reporters for antennal trichoid (at) and antennal intermediate (ai) sensilla (Fig. 3G). For at1 (labeled with Or67d-Gal4-driven mCD8:GFP), the rounder, more weakly labeled Osi8-positive cells were observed adjacent to the OSN somas (Fig. 3G). at4 sensilla distribution (labeled with Or88a-mCD8:GFP) most closely resembled that of the more elongated, strongly expressing Osi8 cells, although in this case the support cell bodies were located distally to the neuronal somas (Fig. 3G). The two intermediate sensillar classes—ai2, labeled with Or83c-mCD8:GFP, and ai3, labeled with Or2a-mCD8:GFP—were intermingled with Osi8-positive cells. While we cannot exclude the possibility that Osi8 is also expressed in support cells in these sensillum types, the lack of a consistent spatial relationship between these sensilla and Osi8-positive cells suggests that Osi8 is expressed in support cells predominantly or exclusively in trichoid sensilla. Consistently, we detected 114 Osi8-positive cells per antenna (±4 SD; n = 5 antennae), which matches well with the total number of at1 and at4 sensilla (68 and 48, respectively; T. O. Auer and L. Abuin, personal communication).
Support cells are named after their developmental roles in construction of sensillar cuticular specializations: trichogen (shaft cell), thecogen (sheath cell), and tormogen (socket cell) (Fig. 3H), although these cells also have functions in mature antennae in secreting high levels of Obps and odorant-degrading enzymes into the sensillum lymph [13, 53]. To determine in which cell type(s) Osi8 is expressed, we performed Osi8 RNA FISH on antennae in which different support cells were labeled with transgenic markers. No overlap of Osi8-positive cells was observed with a marker of thecogen cells, nompA-mCD8:GFP (Fig. 3H). By contrast, both distal and proximal populations of Osi8-positive cells co-localized with the tormogen cell marker ASE5-mCD8:GFP (Fig. 3H and Additional file 1: Figure S1A). Although there is no known marker for the trichogen cell population in the antenna, these observations are consistent with selective Osi8 expression in tormogen cells in trichoid sensilla, most likely corresponding to cluster O3 in the Fly Cell Atlas (Fig. 2B, G).
To examine the localization of Osi8 protein in tormogen cells, we generated a GFP-tagged version of Osi8 and expressed this under the control of an Osi8 promoter-Gal4 driver (Additional file 1: Figure S1B). GFP:Osi8 was detected around the nuclear membrane (presumably the endoplasmic reticulum) and in prominent vesicle-like puncta (Fig. 3I). We also detected GFP within the lumen near the base of the sensillum shaft (Fig. 3I), which may correspond to the occasional protrusions of the tormogen cell and/or extracellular vacuoles that likely derive from this cell type [40]. Confirmation of this localization pattern will require development of specific Osi8 antibodies to detect endogenous protein, but the distribution of GFP:Osi8 in various membranous organelle-like structures is similar to that of other Osi proteins [50, 90, 91].
Osi8 is required for high sensitivity pheromone-evoked responses
To determine the function of Osi8, we generated a null mutant by CRISPR/Cas9-mediated replacement of the Osi8 locus with a DsRed reporter (which was subsequently removed) (Fig. 4A). Homozygous Osi8 mutants are viable and fertile, and the absence of Osi8 transcripts in the antenna was confirmed by RNA FISH (Fig. 4B). Previous work has suggested that Osi proteins have roles in shaping of the cuticle [50, 91]. Notably, in the maxillary palp—a secondary olfactory organ of insects—Osi23 is expressed in developing (but not mature) trichogen cells and is required for the formation of the pores in the sensillar shaft through which odors pass (see “Discussion”) [50]. Scanning electron microscopy of trichoid sensilla of Osi8 mutant antennae did not reveal any overt morphological defects, including in the basal drum (the rounded base of the sensillum secreted by the tormogen cell) (Fig. 4C). Moreover, there were no noticeable defects in the mechanical flexibility of the Osi8 mutant trichoid sensilla as assessed qualitatively during electrophysiological recordings (see below). While we cannot exclude the possibility of more subtle cuticular defects, the different spatio-temporal expression pattern of Osi8 (Fig. 3B) compared to Osi23 (and other family members) suggests that Osi8 plays a role in the mature function, rather than the development, of these sensilla.
To test this hypothesis, we measured the activity of OSNs within trichoid sensilla through single-sensillum electrophysiological recordings. In at4, the responses of Or47b and Or88a neurons to palmitoleic acid and methyl palmitate—their respective cognate ligands [92, 93]—were markedly reduced in Osi8 mutants compared to controls (Fig. 4D–G), although weak responses were observed at the higher concentrations tested. We validated this loss-of-function phenotype in two ways: first, depletion of Osi8 transcripts through RNA interference (RNAi) (Additional file 2: Figure S2A) also led to a decrease in pheromone sensitivity of these neurons (Additional file 2: Figure S2B-E). Second, Osi8-Gal4-driven expression of an Osi8 transgene (but not a negative control CD4:tdTomato transgene) was sufficient to restore pheromone sensitivity in Osi8 mutants to control levels (Fig. 4D–G). By contrast, responses of Or67d neurons to their cognate ligand, (Z)-11-octadecenyl acetate (cis-vaccenyl acetate), were indistinguishable from controls (Fig. 4H,I). Loss of Osi8 did not affect the basal firing rate of any of these neuronal populations (Fig. 4J).
We explored the origin of the electrophysiological phenotype in at4 sensilla by examining the expression of a variety of cellular markers for this sensillum class in Osi8 mutants (Additional file 3: Figure S3). Osi8 cells (labeled by Osi8-Gal4) were indistinguishable in appearance (by confocal microscopy) in the absence of Osi8 to control cells, suggesting that Osi8 is not required for their developmental specification (Additional file 3: Figure S3A). Similarly, a Lush-Gal4 driver (which labels several support cell types in at4 (and at1) sensilla [29, 54, 94]) displayed the same expression in controls and Osi8 mutants (Additional file 3: Figure S3B). Finally, we investigated the expression of Or47b and Or88a using transcriptional reporters, but did not detect any noticeable differences in Osi8 mutants that could explain the attenuated electrophysiological responses of these neurons (Additional file 3: Figure S3C-D).
Together, these results reveal a selective role for Osi8 in promoting high sensitivity of a subset of pheromone-sensing neurons, while not appearing to affect the gross developmental properties of either the neuronal or support cells of these sensilla. The differential requirement for Osi8 in at4 and at1 may reflect differences in the morphology or function of Osi8-expressing cells in these sensillar classes, or another unique property of the at4 sensilla.
Discussion
We have performed a comparative transcriptomic screen to identify genes that are differentially expressed between the D. melanogaster antennal olfactory subsystems, revealing a striking number and diversity of uncharacterized molecules that might contribute to the subsystems’ distinct properties. As RNA samples were collected from adult antennae, it is possible that many of the identified genes have roles in the mature function, rather than the development, of these subsystems. Future transcriptomic analyses of developing ato and amos mutant antennal tissue should be an effective way to identify genes that contribute to the different developmental features of these subsystems.
Our comparative transcriptomic dataset has also been complementary to the Fly Cell Atlas [29] in advancing cell type annotation of support cells in the antenna. Many subsystem-specific genes are expressed in non-neuronal cells, which has allowed us to demarcate those that form part of the Or or Ir subsystem. However, in contrast to OSN populations, we found that very few (or no) unique molecular markers exist for specific support cell types across or within individual sensillar classes. The molecular heterogeneity of support cells suggests a substantial degree of functional overlap, at least within the main types of olfactory sensilla.
We have further exploited these transcriptomic data to demonstrate the selective in situ expression of one of the Or subsystem-specific molecules, Osi8, in tormogen support cells in trichoid sensilla. Importantly, loss-of-function genetic analyses demonstrate a requirement for this protein for pheromone-evoked neuronal responses. Amongst insect olfactory sensory signaling pathways, pheromone detection is well-recognized to require accessory proteins beyond the olfactory receptors [13, 51]. In D. melanogaster, the best characterized are the Obp Lush [56] and two neuronally expressed proteins, the CD36-related Sensory neuron membrane protein 1 (Snmp1) [94, 95] and the DEG/ENaC sodium channel Pickpocket 25 [96]. One important similarity between Osi8 and these proteins is that their loss can be by-passed at high stimulus concentrations [57, 96, 97], indicating that they are not essential components of the signaling cascade but contribute to the high sensitivity and/or selectivity of pheromone detection. Osi8 differs from these other proteins in its function within support cells, and not as a perireceptor protein (like Lush) or in OSNs (like Snmp1 or Pickpocket 25).
While the mechanism of Osi8 is unknown, some insights may be gained from studies on two other Osi proteins, which have been implicated in regulating membrane trafficking in other tissues. Osi21 controls endolysosomal trafficking in photoreceptor neurons, in which it partially localizes to these organelles [90]. Osi23 localizes, at least in part, to endosomes, and the absence of pores in maxillary palp sensillar hairs in Osi23 mutants has been traced back to lack of undulations in the plasma membrane of developing trichogen cells, which may prefigure the porous nature of the secreted cuticle layer [50]. Other Osi proteins appear to be expressed in cuticle-secreting cells in various tissues, and some have been detected in vesicular structures [50, 91]. The biochemical function in any case is, however, unclear. Bearing in mind the caveats of transgenic fusion protein expression, the predominant localization of GFP:Osi8 to the endomembrane system of tormogen cells suggests that Osi8 has an analogous role to other Osi proteins in intracellular membrane trafficking. Future ultrastructural and electrophysiological analysis of Osi8 mutant at4 sensilla will be necessary to determine if this protein contributes to tormogen-specific morphological specializations—such as microvilli and microlamellae bordering the sensillum lymph [40, 41]—and how this might impact the function of pheromone-sensing neurons.
Interestingly, another potential function for Osi8 has emerged from investigations into the molecular basis of octanoic acid resistance of the ecological specialist Drosophila sechellia, which feeds uniquely on octanoic acid-rich “noni” fruit [98]. Osi8 falls within a genomic region identified by unbiased mapping studies as contributing to the tolerance of D. sechellia to this acid [99]. Ubiquitous RNAi of Osi8 in adult or larval D. melanogaster led to a decrease in this species’ resistance to octanoic acid exposure [100, 101]. It is unclear in which tissue Osi8 is required for this role, but we suspect that it is unrelated to its function in the antenna, as octanoic acid does not activate neurons in at4 (or at1) sensilla [49, 102]. Our Osi8 mutant may help validate this RNAi phenotype.
Regardless of its precise function, the conservation of Osi8 across insect taxa—except in P. humanus, which exhibits a drastic loss of Ors—together with the antennal expression of the A. aegypti ortholog, raise the possibility that Osi8 is a conserved regulator of insect pheromone signaling. More generally, the requirement for Osi8 further highlights the important contribution of proteins in non-neuronal cells for the normal function of associated sensory neurons. While there are recent hints of interactions between support cells and neurons in insect sensilla [103], many of the best-characterized examples are found in C. elegans, such as the glial-expressed DEG/ENaC homologs and a chloride channel, which are critical for mechanosensory responses of neighboring nose-touch neurons [104, 105]. Further study of Osi8—as well as the many other non-neuronal genes identified here—should reveal deeper insights into the molecular mechanisms by which specific properties of sensory systems arise from the concerted contributions of neurons and their associated cells.
Conclusions
Together our work provides a resource for characterization of olfactory subsystem-specific genes, emphasizes the importance of support cells in sensory responses of neurons, and identifies a new protein required for insect pheromone detection.
Methods
D. melanogaster culture and strains
Animals were grown on standard Drosophila culture media at 25°C in a 12-h light:12-h dark cycle. Mutant and transgenic lines used are described in Additional file 6: Table S3. Mixed sexes were used in all experiments, except for the electrophysiology, where only male flies were analyzed.
Transcriptomic analyses
Microarrays
Three biological replicates were analyzed for both ato1/Df(3R)p13 and amos3. The mutant lines were back-crossed five times with Oregon-R-P2 prior to the experiment to isogenize the genetic background. Olfactory third antennal segments from ~100–200 adult flies (0–3 days old; mixed sexes (~1:2 males:females—defined by the ratio of ato mutants that successfully eclosed—with no selection for virginity) per biological replicate were harvested via snap-freezing animals in a mini-sieve with liquid nitrogen and agitating the animals to break off the appendages [106]. Collected antennae were homogenized manually with a tissue grinder and total RNA was extracted using a standard TRIzol/chloroform protocol, ethanol precipitated, and dissolved to 45 ng/μl final concentration. Samples were hybridized to Affymetrix Drosophila Genome 2.0 Arrays. Microarray data—which also include an Oregon-R-P2 wild-type genotype not presented in this work—are available at the NCBI Gene Expression Omnibus (GSE183763) [107].
Antennal RNA-seq
Three biological replicates of control (w1118) and amos3 animals were cultured at 22°C. Antennae were collected from 0–4-h-old adults (mixed sexes) by snap-freezing the animals in dry ice and separating the antennae from other tissues by shaking through a 20-μm sieve. Approximately 300 antennae were transferred into a 1.5-ml Eppendorf tube containing 20 μl Trizol. The antennae were homogenized with an RNase-free pestle until no intact cuticle could be detected and the lysate was transferred into 300 μl Trizol. RNA purification was performed following a standard protocol [108] followed by column purification (RNeasy, QIAGEN). RNA-seq libraries were prepared using the Illumina TruSeq Stranded mRNA protocol; library QC was performed using a Fragment Analyzer (Agilent Technologies). Sequence data was processed using the Illumina Pipeline Software v1.82. Library adapters of purity-filtered reads were trimmed with Skewer v0.1.120 [109] and read quality assessed with FastQC v0.10.1 (Babraham Informatics). Reads were aligned against a D. melanogaster reference genome (dmel_r6.02, FlyBase) using TopHat2 v2.0.11 [110]. The numbers of read counts for each Osi gene were summarized with HTSeq-count v0.6.1 [111] using D. melanogaster (GFF:dmel-all-r6.02, FlyBase) gene annotation. RNA-seq data—which also include similar RNA-seq analyses of third instar larval antennal discs not presented in this work—are available at the NCBI Gene Expression Omnibus (GSE190696) [112].
New Drosophila mutant and transgenic lines
Osi8 mutant
The sgRNA expression vector was generated by PCR amplification of three fragments (encoding four different sgRNAs) using the oligonucleotide pairs listed in Additional file 7: Table S4, and cloning these via Gibson assembly into BbsI-digested pCFD5 (Addgene #73914), as described [113]. The donor vector for homologous recombination was generated by amplifying ~1 kb homology arms (HA) flanking the Osi8 coding sequence from genomic DNA of {Act5C-Cas9.P.RFP-}ZH-2A w[118]Lig[169] flies and inserting the products into pHD-DsRed-attP (Addgene #51019) via EcoRI and SacII (HA1) or SapI (HA2) restriction cloning. Injection of the sgRNA vector (150 ng μl−1) and donor vector (400 ng μl−1) into {Act5C-Cas9.P.RFP-}ZH-2A w[118]Lig[169] flies was performed by BestGene Inc. An Osi8 mutant (Osi81) was identified on the basis of DsRed expression and balanced to remove other transgenes; deletion of the Osi8 gene was validated by PCR (Additional file 7: Table S4). The Osi81-DsRed allele was generated by balancing Osi81 with a transgene encoding the Cre recombinase (Additional file 6: Table S3) and selecting for larvae that had lost the DsRed marker.
Osi8-Gal4
3416 bp upstream of the predicted Osi8 transcription start site (i.e., 100 bp 5′ of the start codon) were PCR amplified (Additional file 7: Table S4), cloned in pCRII-TOPO (Thermo Fisher Scientific), sequenced, and subcloned via EcoRI/BamHI sites incorporated into the PCR primers into EcoRI/BglII sites of pGAL4attB [31]. The construct was inserted into attP2 by phiC31-mediated transgenesis (BestGene Inc.).
UAS-Osi8
A genomic region encompassing the Osi8 coding sequence (with 5′ and 3′ UTRs and intron) was PCR amplified (Additional file 7: Table S4) from Canton-S genomic DNA, cloned into pCR2.1-TOPO and sequenced, before subcloning using EcoRI (using the EcoRI site incorporated into the forward primer and, at the 3′ end, the EcoRI site in pCR2.1-TOPO) into pUAST-attB [114]. The construct was inserted into attP40 by phiC31-mediated transgenesis (BestGene Inc.).
UAS-SS:EGFP:Osi8 (GFP:Osi8)
The Osi8 cDNA lacking the first 23 codons (encoding the predicted Osi8 signal sequence) was PCR amplified from Canton-S antennal cDNA (Additional file 7: Table S4), cloned into pCR2.1-TOPO and sequenced, before subcloning using EcoRI and XbaI into pUAST-SS:EGFP attB [26], creating an in-frame fusion with the heterologous signal sequence from calreticulin and the EGFP tag. The construct was inserted into attP40 by phiC31-mediated transgenesis (BestGene Inc.).
Bioinformatic analyses and other RNA-sequence datasets
Gene ontology term enrichment analyses were performed using GOrilla [68], using genes displaying a >2-fold expression difference between ato and amos; the Affymetrix Drosophila Genome 2.0 Array gene list was used as the reference dataset. Subsequently, genes displaying >4-fold expression difference between these genotypes were curated manually, using annotations from FlyBase [115] and additional analyses using SMART [116] and BLAST [117], the latter usually to verify the absence of similarities to protein domains of known function. D. melanogaster antennal and head scRNA-seq data were from the 10× stringent dataset from the Fly Cell Atlas [29]; these were visualized in the HVG tSNE coordinate representation in the SCope interface [118]. Expression data across tissues/life stages for Osi8 were from Fly Atlas 2 [88]. Expression data of the Aedes aegypti Osi8 ortholog (AAEL004275, AaegL3.3 annotation) were from a published dataset [89].
Histology
RNA FISH (based upon the Tyramide Signal Amplification™ method (Perkin Elmer)) and protein immunofluorescence on whole-mount antennae or antennal cryosections were performed essentially as described [106], using 1–10-day-old animals. Templates for RNA FISH probes were amplified by PCR from genomic DNA (Additional file 7: Table S4) and cloned into pGEM-T Easy (Promega) or pCRII-TOPO; the Osi8 probe was synthesized from the same genomic region used to construct the UAS-Osi8 transgene. Antibodies used are listed in Additional file 8: Table S5. Imaging was performed on a Zeiss LSM710 or LSM880 Airyscan confocal microscope using 40× or 63× oil immersion objectives. Confocal stacks were imported into Fiji [119] for processing and analysis and subsequently formatted in Adobe Photoshop 2022. Cell counting was performed manually using the Cell Counter plugin of ImageJ. For each experiment, the expression/phenotype was assessed in antennae from at least 10 animals in at least two independent replicates. Experimenters were not blind to the gene or genotype under assessment.
Electron microscopy
Control (w1118) and Osi81 2-day or 20-day-old animals were fixed in 2.5% glutaraldehyde solution in 0.1 M phosphate buffer pH 7.4 (PB) for 2 h at room temperature (RT). Samples were rinsed 3 × 5 min in PB and post-fixed in a fresh mixture of 1% osmium tetroxide (EMS) with 1.5% potassium ferrocyanide in PB for 2 h at RT. After washing twice in distilled water, samples were dehydrated in acetone solution at graded concentrations (30% (40 min), 50% (40 min), 70% (40 min), 100% (3 × 1 h)). Samples were dried at the critical point (BALTEC CPD 30) and glued on a pin and sputter coated with 10 nm of metal platinum (LEICA EM SCD 500). Micrographs were taken with a scanning electron microscope FEI Quanta FEG 250 (FEI, Eindhoven, the Netherlands) with an Everhart-Thornley Detector (ETD).
Electrophysiology
For all electrophysiological experiments, 7-day-old male flies (housed in groups of 10) were used. An animal was prepared for recordings by wedging it into the narrow end of a truncated 200-μl plastic pipette tip to expose the antennae, one of which was then stabilized between a tapered glass microcapillary tube and a coverslip covered with double-sided tape. Single-unit recordings were performed essentially as described [120]. In brief, the electrical activity of the neurons was recorded extracellularly by inserting a sharp electrode filled with artificial hemolymph solution [121] into either an at4 sensillum (Or47b and Or88a OSN recordings) or an at1 sensillum (Or67d OSN recording). A reference electrode filled with the same solution was inserted into the eye. Odor stimuli: trans-palmitoleic acid (Cayman Chemical, CAS 10030-73-6) and cis-vaccenyl acetate (Cayman Chemical, CAS 6186-98-7) were diluted in ethanol, and methyl palmitate (Sigma-Aldrich, CAS 112-39-0) was diluted in paraffin oil. 4.5 μl aliquots of dilutions were spotted onto filter discs and delivered via a 500-ms air pulse at 250 ml/min directly to the antenna from close range, as described [120]. Ethanol was allowed to evaporate before the experiments. Spike responses were averaged, binned at 50 ms, and smoothed using a binomial algorithm to obtain peristimulus time histograms (PSTHs). For dose-response curves and statistical analysis, responses were quantified by subtracting the pre-stimulus spike rate (1 s) from the peak spike response during odorant stimulation (adjusted peak responses). Basal neuronal activity was measured during 10 s from flies that had not been exposed to odor stimuli. Electrophysiological data are provided in Additional file 10: Data S2.
Availability of data and materials
All data generated or analyzed during this study are included in this published article, its supplementary information files (Additional file 10: Data S2) and publicly available repositories: NCBI Gene Expression Omnibus GSE183763 (comparative antennal microarray data) and GSE190696 (bulk antennal RNA-seq data), SCope (Fly Cell Atlas data), and FlyAtlas 2 (tissue-specific bulk RNA-seq data).
References
Su CY, Menuz K, Carlson JR. Olfactory perception: receptors, cells, and circuits. Cell. 2009;139(1):45–59.
Ache BW, Young JM. Olfaction: diverse species, conserved principles. Neuron. 2005;48(3):417–30.
Bear DM, Lassance JM, Hoekstra HE, Datta SR. The evolving neural and genetic architecture of vertebrate olfaction. Curr Biol. 2016;26(20):R1039–49.
Munger SD, Leinders-Zufall T, Zufall F. Subsystem organization of the mammalian sense of smell. Annu Rev Physiol. 2009;71:115–40.
Breer H, Fleischer J, Strotmann J. The sense of smell: multiple olfactory subsystems. Cell Mol Life Sci. 2006;63(13):1465–75.
Ishii KK, Touhara K. Neural circuits regulating sexual behaviors via the olfactory system in mice. Neurosci Res. 2019;140:59–76.
Vosshall LB, Stocker RF. Molecular architecture of smell and taste in Drosophila. Annu Rev Neurosci. 2007;30:505–33.
Grabe V, Sachse S. Fundamental principles of the olfactory code. Bio Systems. 2018;164:94–101.
Schlegel P, Bates AS, Sturner T, Jagannathan SR, Drummond N, Hsu J, et al. Information flow, cell types and stereotypy in a full olfactory connectome. Elife. 2021;10:e66018.
Silbering AF, Rytz R, Grosjean Y, Abuin L, Ramdya P, Jefferis GS, et al. Complementary function and integrated wiring of the evolutionarily distinct Drosophila olfactory subsystems. J Neurosci. 2011;31(38):13357–75.
Grabe V, Baschwitz A, Dweck HKM, Lavista-Llanos S, Hansson BS, Sachse S. Elucidating the neuronal architecture of olfactory glomeruli in the Drosophila antennal lobe. Cell Rep. 2016;16(12):3401–13.
Shanbhag SR, Muller B, Steinbrecht RA. Atlas of olfactory organs of Drosophila melanogaster. 1. Types, external organization, innervation and distribution of olfactory sensilla. Int J Insect Morphol Embryol. 1999;28(4):377–97.
Schmidt HR, Benton R. Molecular mechanisms of olfactory detection in insects: beyond receptors. Open Biol. 2020;10(10):200252.
Vosshall LB, Amrein H, Morozov PS, Rzhetsky A, Axel R. A spatial map of olfactory receptor expression in the Drosophila antenna. Cell. 1999;96(5):725–36.
Gao Q, Chess A. Identification of candidate Drosophila olfactory receptors from genomic DNA sequence. Genomics. 1999;60(1):31–9.
Clyne PJ, Warr CG, Freeman MR, Lessing D, Kim J, Carlson JR. A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron. 1999;22(2):327–38.
Couto A, Alenius M, Dickson BJ. Molecular, anatomical, and functional organization of the Drosophila olfactory system. Curr Biol. 2005;15(17):1535–47.
Fishilevich E, Vosshall LB. Genetic and functional subdivision of the Drosophila antennal lobe. Curr Biol. 2005;15(17):1548–53.
Benton R, Vannice KS, Gomez-Diaz C, Vosshall LB. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell. 2009;136(1):149–62.
Wicher D, Schafer R, Bauernfeind R, Stensmyr MC, Heller R, Heinemann SH, et al. Drosophila odorant receptors are both ligand-gated and cyclic-nucleotide-activated cation channels. Nature. 2008;452(7190):1007–11.
Sato K, Pellegrino M, Nakagawa T, Nakagawa T, Vosshall LB, Touhara K. Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature. 2008;452(7190):1002–6.
Butterwick JA, Del Marmol J, Kim KH, Kahlson MA, Rogow JA, Walz T, et al. Cryo-EM structure of the insect olfactory receptor Orco. Nature. 2018;560(7719):447–52.
Del Marmol J, Yedlin MA, Ruta V. The structural basis of odorant recognition in insect olfactory receptors. Nature. 2021;597(7874):126–31.
Larsson MC, Domingos AI, Jones WD, Chiappe ME, Amrein H, Vosshall LB. Or83b encodes a broadly expressed odorant receptor essential for Drosophila olfaction. Neuron. 2004;43(5):703–14.
Benton R, Sachse S, Michnick SW, Vosshall LB. Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 2006;4(2):e20.
Abuin L, Bargeton B, Ulbrich MH, Isacoff EY, Kellenberger S, Benton R. Functional architecture of olfactory ionotropic glutamate receptors. Neuron. 2011;69(1):44–60.
Vulpe A, Menuz K. Ir76b is a co-receptor for amine responses in Drosophila olfactory neurons. Front Cell Neurosci. 2021;15:759238.
Vosshall LB, Wong AM, Axel R. An olfactory sensory map in the fly brain. Cell. 2000;102(2):147–59.
Li H, Janssens J, De Waegeneer M, Kolluru SS, Davie K, Gardeux V, et al. Fly Cell Atlas: a single-nucleus transcriptomic atlas of the adult fruit fly. Science. 2022;375(6584):eabk2432.
Knecht ZA, Silbering AF, Cruz J, Yang L, Croset V, Benton R, et al. Ionotropic Receptor-dependent moist and dry cells control hygrosensation in Drosophila. Elife. 2017;6:e26654.
Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, et al. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010;6(8):e1001064.
Task D, Lin C-C, Vulpe A, Afify A, Ballou S, Brbić M, et al. Chemoreceptor co-expression in Drosophila melanogaster olfactory neurons. Elife. 2022;11:e72599.
Gupta BP, Rodrigues V. Atonal is a proneural gene for a subset of olfactory sense organs in Drosophila. Genes Cells. 1997;2(3):225–33.
zur Lage PI, Prentice DR, Holohan EE, Jarman AP. The Drosophila proneural gene amos promotes olfactory sensillum formation and suppresses bristle formation. Development. 2003;130(19):4683–93.
Goulding SE, zur Lage P, Jarman AP. amos, a proneural gene for Drosophila olfactory sense organs that is regulated by lozenge. Neuron. 2000;25(1):69–78.
Frank DD, Enjin A, Jouandet GC, Zaharieva EE, Para A, Stensmyr MC, et al. Early integration of temperature and humidity stimuli in the Drosophila brain. Curr Biol. 2017;27(15):2381–2388 e2384.
Enjin A, Zaharieva EE, Frank DD, Mansourian S, Suh GS, Gallio M, et al. Humidity sensing in Drosophila. Curr Biol. 2016;26(10):1352–8.
Ni L, Bronk P, Chang EC, Lowell AM, Flam JO, Panzano VC, et al. A gustatory receptor paralogue controls rapid warmth avoidance in Drosophila. Nature. 2013;500(7464):580–4.
Knecht ZA, Silbering AF, Ni LN, Klein M, Budelli G, Bell R, et al. Distinct combinations of variant ionotropic glutamate receptors mediate thermosensation and hygrosensation in Drosophila. Elife. 2016;5:e17879.
Nava Gonzales C, McKaughan Q, Bushong EA, Cauwenberghs K, Ng R, Madany M, et al. Systematic morphological and morphometric analysis of identified olfactory receptor neurons in Drosophila melanogaster. Elife. 2021;10:e69896.
Shanbhag SR, Muller B, Steinbrecht RA. Atlas of olfactory organs of Drosophila melanogaster 2. Internal organization and cellular architecture of olfactory sensilla. Arthropod Struct Dev. 2000;29(3):211–29.
Shanbhag SR, Singh K, Singh RN. Fine structure and primary sensory projections of sensilla located in the sacculus of the antenna of Drosophila melanogaster. Cell Tissue Res. 1995;282(2):237–49.
Getahun MN, Wicher D, Hansson BS, Olsson SB. Temporal response dynamics of Drosophila olfactory sensory neurons depends on receptor type and response polarity. Front Cell Neurosci. 2012;6:54.
Cao LH, Jing BY, Yang D, Zeng X, Shen Y, Tu Y, et al. Distinct signaling of Drosophila chemoreceptors in olfactory sensory neurons. Proc Natl Acad Sci U S A. 2016;113(7):E902–11.
Grosjean Y, Rytz R, Farine JP, Abuin L, Cortot J, Jefferis GS, et al. An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature. 2011;478(7368):236–40.
Hallem EA, Ho MG, Carlson JR. The molecular basis of odor coding in the Drosophila antenna. Cell. 2004;117(7):965–79.
Yao CA, Ignell R, Carlson JR. Chemosensory coding by neurons in the coeloconic sensilla of the Drosophila antenna. J Neurosci. 2005;25(37):8359–67.
Dobritsa AA, van der Goes van Naters W, Warr CG, Steinbrecht RA, Carlson JR. Integrating the molecular and cellular basis of odor coding in the Drosophila antenna. Neuron. 2003;37(5):827–41.
Hallem EA, Carlson JR. Coding of odors by a receptor repertoire. Cell. 2006;125(1):143–60.
Ando T, Sekine S, Inagaki S, Misaki K, Badel L, Moriya H, et al. Nanopore formation in the cuticle of an insect olfactory sensillum. Curr Biol. 2019;29(9):1512–1520 e1516.
Leal WS. Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu Rev Entomol. 2013;58:373–91.
Xiao S, Sun JS, Carlson JR. Robust olfactory responses in the absence of odorant binding proteins. Elife. 2019;8:e51040.
Sun JS, Xiao S, Carlson JR. The diverse small proteins called odorant-binding proteins. Open Biol. 2018;8(12):180208.
Larter NK, Sun JS, Carlson JR. Organization and function of Drosophila odorant binding proteins. Elife. 2016;5:e20242.
Scheuermann EA, Smith DP. Odor-specific deactivation defects in a Drosophila odorant-binding protein mutant. Genetics. 2019;213(3):897–909.
Xu P, Atkinson R, Jones DN, Smith DP. Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron. 2005;45(2):193–200.
Gomez-Diaz C, Reina JH, Cambillau C, Benton R. Ligands for pheromone-sensing neurons are not conformationally activated odorant binding proteins. PLoS Biol. 2013;11(4):e1001546.
Okumura M, Kato T, Miura M, Chihara T. Hierarchical axon targeting of Drosophila olfactory receptor neurons specified by the proneural transcription factors Atonal and Amos. Genes Cells. 2016;21(1):53–64.
Mohapatra P, Menuz K. Molecular profiling of the Drosophila antenna reveals conserved genes underlying olfaction in insects. G3 (Bethesda). 2019;9(11):3753–71.
Menuz K, Larter NK, Park J, Carlson JR. An RNA-Seq screen of the Drosophila antenna identifies a transporter necessary for ammonia detection. PLoS Genet. 2014;10(11):e1004810.
Sanchez-Alcaniz JA, Silbering AF, Croset V, Zappia G, Sivasubramaniam AK, Abuin L, et al. An expression atlas of variant ionotropic glutamate receptors identifies a molecular basis of carbonation sensing. Nat Commun. 2018;9(1):4252.
Stewart S, Koh TW, Ghosh AC, Carlson JR. Candidate ionotropic taste receptors in the Drosophila larva. PNAS. 2015;112(14):4195–201.
Koh TW, He Z, Gorur-Shandilya S, Menuz K, Larter NK, Stewart S, et al. The Drosophila IR20a clade of Ionotropic Receptors are candidate taste and pheromone receptors. Neuron. 2014;83(4):850–65.
Chen YD, Dahanukar A. Recent advances in the genetic basis of taste detection in Drosophila. Cell Mol Life Sci. 2020;77(6):1087–101.
Kwon JY, Dahanukar A, Weiss LA, Carlson JR. The molecular basis of CO2 reception in Drosophila. Proc Natl Acad Sci U S A. 2007;104(9):3574–8.
Jones WD, Cayirlioglu P, Grunwald Kadow I, Vosshall LB. Two chemosensory receptors together mediate carbon dioxide detection in Drosophila. Nature. 2007;445:86–90.
Fujii S, Yavuz A, Slone J, Jagge C, Song X, Amrein H. Drosophila sugar receptors in sweet taste perception, olfaction, and internal nutrient sensing. Curr Biol. 2015;25(5):621–7.
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinf. 2009;10:48.
Younus F, Chertemps T, Pearce SL, Pandey G, Bozzolan F, Coppin CW, et al. Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster. Insect Biochem Mol Biol. 2014;53:30–43.
Sun JS, Larter NK, Chahda JS, Rioux D, Gumaste A, Carlson JR. Humidity response depends on the small soluble protein Obp59a in Drosophila. Elife. 2018;7:e39249.
Vulpe A, Kim HS, Ballou S, Wu ST, Grabe V, Nava Gonzales C, et al. An ammonium transporter is a non-canonical olfactory receptor for ammonia. Curr Biol. 2021;31(15):3382–3390 e3387.
Liu YC, Pearce MW, Honda T, Johnson TK, Charlu S, Sharma KR, et al. The Drosophila melanogaster phospholipid flippase dATP8B is required for odorant receptor function. PLoS Genet. 2014;10(3):e1004209.
Ha TS, Xia R, Zhang H, Jin X, Smith DP. Lipid flippase modulates olfactory receptor expression and odorant sensitivity in Drosophila. Proc Natl Acad Sci U S A. 2014;111(21):7831–6.
Schaefer E, Zaloszyc A, Lauer J, Durand M, Stutzmann F, Perdomo-Trujillo Y, et al. Mutations in SDCCAG8/NPHP10 cause Bardet-Biedl syndrome and are associated with penetrant renal disease and absent polydactyly. Mol Syndromol. 2011;1(6):273–81.
Cinege G, Zsamboki J, Vidal-Quadras M, Uv A, Csordas G, Honti V, et al. Genes encoding cuticular proteins are components of the Nimrod gene cluster in Drosophila. Insect Biochem Mol Biol. 2017;87:45–54.
Li H, Li T, Horns F, Li J, Xie Q, Xu C, et al. Single-cell transcriptomes reveal diverse regulatory strategies for olfactory receptor expression and axon targeting. Curr Biol. 2020;30(7):1189–1198 e1185.
Hartenstein V. The muscle pattern of Drosophila. In: Sink H, editor. Muscle Development in Drosophila. New York: Springer New York; 2006. p. 8–27.
Mamiya A, Straw AD, Tomasson E, Dickinson MH. Active and passive antennal movements during visually guided steering in flying Drosophila. J Neurosci. 2011;31(18):6900–14.
Sen A, Kuruvilla D, Pinto L, Sarin A, Rodrigues V. Programmed cell death and context dependent activation of the EGF pathway regulate gliogenesis in the Drosophila olfactory system. Mech Dev. 2004;121(1):65–78.
Pikielny CW, Hasan G, Rouyer F, Rosbash M. Members of a family of Drosophila putative odorant-binding proteins are expressed in different subsets of olfactory hairs. Neuron. 1994;12(1):35–49.
McKenna MP, Hekmat-Scafe DS, Gaines P, Carlson JR. Putative Drosophila pheromone-binding proteins expressed in a subregion of the olfactory system. J Biol Chem. 1994;269(23):16340–7.
Steiner C, Bozzolan F, Montagne N, Maibeche M, Chertemps T. Neofunctionalization of "juvenile hormone esterase duplication" in Drosophila as an odorant-degrading enzyme towards food odorants. Sci Rep. 2017;7(1):12629.
Senthilan PR, Piepenbrock D, Ovezmyradov G, Nadrowski B, Bechstedt S, Pauls S, et al. Drosophila auditory organ genes and genetic hearing defects. Cell. 2012;150(5):1042–54.
Smith CR, Morandin C, Noureddine M, Pant S. Conserved roles of Osiris genes in insect development, polymorphism and protection. J Evol Biol. 2018;31(4):516–29.
Dorer DR, Rudnick JA, Moriyama EN, Christensen AC. A family of genes clustered at the Triplo-lethal locus of Drosophila melanogaster has an unusual evolutionary history and significant synteny with Anopheles gambiae. Genetics. 2003;165(2):613–21.
Shah N, Dorer DR, Moriyama EN, Christensen AC. Evolution of a large, conserved, and syntenic gene family in insects. G3 (Bethesda). 2012;2(2):313–9.
Kirkness EF, Haas BJ, Sun W, Braig HR, Perotti MA, Clark JM, et al. Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle. Proc Natl Acad Sci U S A. 2010;107(27):12168–73.
Krause SA, Overend G, Dow JAT, Leader DP. FlyAtlas 2 in 2022: enhancements to the Drosophila melanogaster expression atlas. Nucleic Acids Res. 2022;50(D1):D1010–5.
Matthews BJ, McBride CS, DeGennaro M, Despo O, Vosshall LB. The neurotranscriptome of the Aedes aegypti mosquito. BMC Genomics. 2016;17:32.
Lee J, Song M, Hong S. Negative regulation of the novel norpA(P24) suppressor, diehard4, in the endo-lysosomal trafficking underlies photoreceptor cell degeneration. PLoS Genet. 2013;9(6):e1003559.
Scholl A, Yang Y, McBride P, Irwin K, Jiang L. Tracheal expression of Osiris gene family in Drosophila. Gene Expr Patterns. 2018;28:87–94.
Dweck HK, Ebrahim SA, Thoma M, Mohamed AA, Keesey IW, Trona F, et al. Pheromones mediating copulation and attraction in Drosophila. Proc Natl Acad Sci U S A. 2015;112(21):E2829–35.
Lin HH, Cao DS, Sethi S, Zeng Z, Chin JS, Chakraborty TS, et al. Hormonal modulation of pheromone detection enhances male courtship success. Neuron. 2016;90(6):1272–85.
Benton R, Vannice KS, Vosshall LB. An essential role for a CD36-related receptor in pheromone detection in Drosophila. Nature. 2007;450(7167):289–93.
Jin X, Ha TS, Smith DP. SNMP is a signaling component required for pheromone sensitivity in Drosophila. PNAS. 2008;105(31):10996–1001.
Ng R, Salem SS, Wu ST, Wu M, Lin HH, Shepherd AK, et al. Amplification of Drosophila olfactory responses by a DEG/ENaC channel. Neuron. 2019;104(5):947–959 e945.
Li Z, Ni JD, Huang J, Montell C. Requirement for Drosophila SNMP1 for rapid activation and termination of pheromone-induced activity. PLoS Genet. 2014;10(9):e1004600.
Auer TO, Shahandeh MP, Benton R. Drosophila sechellia: a genetic model for behavioral evolution and neuroecology. Annu Rev Genet. 2021;55:527–54.
Hungate EA, Earley EJ, Boussy IA, Turissini DA, Ting CT, Moran JR, et al. A locus in Drosophila sechellia affecting tolerance of a host plant toxin. Genetics. 2013;195(3):1063–75.
Andrade Lopez JM, Lanno SM, Auerbach JM, Moskowitz EC, Sligar LA, Wittkopp PJ, et al. Genetic basis of octanoic acid resistance in Drosophila sechellia: functional analysis of a fine-mapped region. Mol Ecol. 2017;26(4):1148–60.
Lanno SM, Shimshak SJ, Peyser RD, Linde SC, Coolon JD. Investigating the role of Osiris genes in Drosophila sechellia larval resistance to a host plant toxin. Ecol Evol. 2019;9(4):1922–33.
van der Goes van Naters W, Carlson JR. Receptors and neurons for fly odors in Drosophila. Curr Biol. 2007;17(7):606–12.
Prelic S, Pal Mahadevan V, Venkateswaran V, Lavista-Llanos S, Hansson BS, Wicher D. Functional interaction between Drosophila olfactory sensory neurons and their support cells. Front Cell Neurosci. 2021;15:789086.
Han L, Wang Y, Sangaletti R, D'Urso G, Lu Y, Shaham S, et al. Two novel DEG/ENaC channel subunits expressed in glia are needed for nose-touch sensitivity in Caenorhabditis elegans. J Neurosci. 2013;33(3):936–49.
Fernandez-Abascal J, Johnson CK, Graziano B, Wang L, Encalada N, Bianchi L. A glial ClC Cl(-) channel mediates nose touch responses in C. elegans. Neuron. 2022;110(3):470–485 e477.
Saina M, Benton R. Visualizing olfactory receptor expression and localization in Drosophila. Methods Mol Biol. 2013;1003:211–28.
Saina M, Benton R: Comparative transcriptomic analysis of the Drosophila melanogaster antennal olfactory subsystems. Genbank 2021, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE183763.
Rio DC, Ares M Jr, Hannon GJ, Nilsen TW. Purification of RNA using TRIzol (TRI reagent). Cold Spring Harb Protoc. 2010;2010(6):pdb.prot5439.
Jiang H, Lei R, Ding SW, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinf. 2014;15:182.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Armida J, Benton J: Transcriptomic analysis of wild-type and amos mutant D. melanogaster third-instar larval antennal discs and adult antennae. Genbank 2021, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE190696.
Port F, Bullock SL. Augmenting CRISPR applications in Drosophila with tRNA-flanked sgRNAs. Nat Methods. 2016;13(10):852–4.
Bischof J, Maeda RK, Hediger M, Karch F, Basler K. An optimized transgenesis system for Drosophila using germ-line-specific phiC31 integrases. Proc Natl Acad Sci U S A. 2007;104(9):3312–7.
Larkin A, Marygold SJ, Antonazzo G, Attrill H, Dos Santos G, Garapati PV, et al. FlyBase: updates to the Drosophila melanogaster knowledge base. Nucleic Acids Res. 2021;49(D1):D899–907.
Letunic I, Khedkar S, Bork P. SMART: recent updates, new developments and status in 2020. Nucleic Acids Res. 2021;49(D1):D458–60.
Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36(Web Server issue):W5–9.
Davie K, Janssens J, Koldere D, De Waegeneer M, Pech U, Kreft L, et al. A single-cell transcriptome atlas of the aging Drosophila brain. Cell. 2018;174(4):982–998 e920.
Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82.
Ng R, Lin HH, Wang JW, Su CY. Electrophysiological recording from Drosophila trichoid sensilla in response to odorants of low volatility. J Vis Exp. 2017;125(125):56147.
Wang JW, Wong AM, Flores J, Vosshall LB, Axel R. Two-photon calcium imaging reveals an odor-evoked map of activity in the fly brain. Cell. 2003;112(2):271–82.
Jarman AP, Grell EH, Ackerman L, Jan LY, Jan YN. Atonal is the proneural gene for Drosophila photoreceptors. Nature. 1994;369(6479):398–400.
Kurtovic A, Widmer A, Dickson BJ. A single class of olfactory neurons mediates behavioural responses to a Drosophila sex pheromone. Nature. 2007;446(7135):542–6.
Chung YD, Zhu J, Han Y, Kernan MJ. nompA encodes a PNS-specific, ZP domain protein required to connect mechanosensory dendrites to sensory structures. Neuron. 2001;29(2):415–28.
Acknowledgements
We are very grateful to Antonio Mucciolo (Electron Microscopy Facility, UNIL) for performing the electron microscopy analysis, Sylvain Pradervand and Karolina Bojkowska (Lausanne Genomic Technologies Facility, UNIL) for microarray and RNA-sequencing analyses and Michele Marconcini for advice on Gene Ontology analysis. We thank John Carlson, Andrew Jarman, the Bloomington Drosophila Stock Center (NIH P40OD018537), the National Institute of Genetics Fly Stock Center (Japan), and the Developmental Studies Hybridoma Bank (NICHD of the NIH, University of Iowa) for reagents. We thank members of the Benton lab for discussions and comments on the manuscript.
Funding
Research in C-Y.S.’s laboratory was supported by NIH grants (R01DC016466, R21DC018912, R21DC020536, and R21AI169343). Research in R.B.’s laboratory was supported by ERC Consolidator and Advanced Grants (615094 and 833548, respectively) and the Swiss National Science Foundation.
Author information
Authors and Affiliations
Contributions
MSc generated molecular reagents, performed histological, bioinformatic, and pilot electrophysiological analyses. RN performed electrophysiological analyses. SC generated molecular reagents and performed histological analyses. MSa performed the microarray screen, generated molecular reagents, and performed pilot bioinformatic and histological analyses. JA performed the bulk antennal RNA-sequencing analysis. CYS supervised electrophysiological experiments. RB supervised the project, performed bioinformatic and image analyses, and wrote the manuscript, with feedback from other authors. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Osi8 tormogen cell expression and validation of Osi8-Gal4. (A) GFP immunofluorescence and Osi8 RNA FISH on a whole-mount antenna in which the tormogen support cells are labeled (UAS-mCD8:GFP;ASE5-Gal4). The white arrowheads point to examples of proximally located cells expressing both Osi8 and GFP. Scale bar, 20 μm. (B) GFP immunofluorescence and Osi8 RNA FISH on antennal cryosections of UAS-mCD8:GFP/+;Osi8-Gal4/+ animals. Scale bar, 20 μm. The white arrowheads point to examples of cells expressing GFP where Osi8 transcripts are not detected. Of 176 cells analyzed in four antennae, 135 (77%) express both Osi8 RNA and GFP, while 41 (23%) express only GFP. The latter category may be due to the lower sensitivity of RNA FISH and/or ectopic Osi8-Gal4 expression. Note also the heterogeneous (and often uncorrelated) expression level of Osi8 RNA and GFP.
Additional file 2: Figure S2.
Phenotypic analysis of Osi8 RNAi. (A) Osi8 RNA FISH on control (Act5C-Gal4,UAS-Dcr-2/+) and Osi8RNAi (UAS-Osi8RNAi/+;Act5C-Gal4,UAS-Dcr-2/+) whole-mount antennae. Scale bar, 20 μm. (B) Representative traces of electrophysiological responses of Or47b OSNs to palmitoleic acid (10-1 v/v) (0.5 s stimulus, gray bar) in Gal4 control (Act5C-Gal4,UAS-Dcr-2/+), UAS control (UAS-Osi8RNAi/+) and Osi8RNAi (UAS-Osi8RNAi/+; Act5C-Gal4,UAS-Dcr-2/+) animals. Raster plots and PSTHs of these responses are shown below each trace. (C) Dose-response curves of Or47b OSN responses to palmitoleic acid; genotypes are color-coded as in (B). Mean responses ±SEM are plotted; n = 12 sensilla; 4-6 flies. Statistical comparisons between genotypes were performed by two-way ANOVA: NS P > 0.05, *** P < 0.001. Full statistical analyses are provided in Additional file 10: Data S2. (D) Representative traces, raster plots and PSTHs of electrophysiological responses of Or88a OSNs to methyl palmitate (10-1 v/v) (0.5 s stimulus, gray bar) in the same genotypes as in (B). (E) Dose-response curves of Or88a OSN responses to methyl palmitate; genotypes are color-coded as in (B). Mean responses ±SEM are plotted; n = 12 sensilla; 4-6 flies. Statistical analyses were performed as in (C).
Additional file 3: Figure S3.
Cellular phenotypic characterization of Osi8 mutants. (A) GFP immunofluorescence labeling Osi8 cells in control (UAS-mCD8:GFP/+;Osi8-Gal4/+) and Osi8 mutant (UAS-mCD8:GFP/+;Osi8-Gal4,Osi81/Osi81) antennae. Scale bar, 20 μm. (B) GFP immunofluorescence labeling Lush cells in control (Lush-Gal4/UAS-mCD8:GFP) and Osi8 mutant (Lush-Gal4/UAS-mCD8:GFP;Osi81) antennae. Scale bar, 20 μm. (C) GFP immunofluorescence labeling Or47b neurons in control (Or47b-Gal4/UAS-mCD8:GFP) and Osi8 mutant (Or47b-Gal4/UAS-mCD8:GFP;Osi81) antennae. Scale bar, 20 μm. (D) GFP immunofluorescence labeling Or88a neurons in control (Or88a-mCD8:GFP/+) and Osi8 mutant (Or88a-mCD8:GFP/+;Osi81) antennae. Scale bar, 20 μm.
Additional file 4: Table S1.
Analysis of Ir and Or subsystem-enriched genes. Annotations of genes enriched in the Ir (magenta) and Or (green) subsystems, as defined as those showing a >4-fold differential expression between amos and ato mutant antennae. This arbitrary cut-off captures all of the known chemosensory receptors expressed in the two subsystems. The dataset was cleaned from the original microarray data (GSE183763) by updating gene names and removing duplicate gene entries and those corresponding to transposons. Neuronal/non-neuronal expression is a qualitative assessment based upon the Fly Cell Atlas antennal dataset [29] and/or RNA FISH data for chemosensory receptors [17,18,19]: “1” = expressed; “0” = not expressed; empty cells indicate robust expression was not detected in either cell type (which may reflect sub-threshold expression level). Breadth of expression provides an approximate categorization of cell-type specific and ubiquitously expressed genes within each olfactory subsystem, based upon the qualitative assessment of the number of clusters in the Fly Cell Atlas antennal dataset in which a gene is detected: “1” = 1 cluster; “2” = 2-5 clusters; “3” = >5 clusters.
Additional file 5: Table S2.
Definition of Ir and Or subsystem support cell populations. Populations were numbered arbitrarily, based upon the Fly Cell Atlas antennal cell cluster representation [29]; it is likely that they can be subclustered further. The “I” and “O” categorization was defined post-hoc, based upon a qualitative assessment of their preferential expression of Ir subsystem- or Or subsystem-enriched genes (from Additional file 4: Table S1; undetected genes were excluded). The I12 and I13 clusters express Obp84a, which is thought to be expressed in thecogen cells based upon exclusive co-expression with a nompA promoter transgenic reporter [54]. Endogenous nompA transcripts are detected consistently (albeit weakly) in I12 but only very sparsely in I13.
Additional file 7: Table S4.
Oligonucleotides.
Additional file 8: Table S5.
Antibodies.
Additional file 9: Data S1.
Gene ontology (GO) term analyses of Or and Ir subsystem-enriched genes. Enriched GO terms for Or subsystem genes (>2-fold more highly expressed in ato versus amos antennae from the microarray analysis) and Ir subsystem genes (>2-fold more highly expressed in amos versus ato antennae), organized by ontology (Process, Function, Component) on different worksheets. Data represent the output from GOrilla [68], including the GO term and gene lists, enrichment terms and statistical analyses (see the bottom of the first worksheet for the definitions of values). Below each table is the graphical visualization of the hierarchy of enriched GO terms (from GOrilla), to emphasize the statistically significant, enriched nodes.
Additional file 10: Data S2.
Electrophysiological quantifications. Raw spike counts, processed data and statistical comparisons are shown for all electrophysiological experiments, organized by figure panel on different worksheets.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Scalzotto, M., Ng, R., Cruchet, S. et al. Pheromone sensing in Drosophila requires support cell-expressed Osiris 8. BMC Biol 20, 230 (2022). https://doi.org/10.1186/s12915-022-01425-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12915-022-01425-w