Olfactory receptor accessory proteins play crucial roles in receptor function and gene choice
Abstract
Each of the olfactory sensory neurons (OSNs) chooses to express a single G protein-coupled olfactory receptor (OR) from a pool of hundreds. Here, we show the receptor transporting protein (RTP) family members play a dual role in both normal OR trafficking and determining OR gene choice probabilities. Rtp1 and Rtp2 double knockout mice (RTP1,2DKO) show OR trafficking defects and decreased OSN activation. Surprisingly, we discovered a small subset of the ORs are expressed in larger numbers of OSNs despite the presence of fewer total OSNs in RTP1,2DKO. Unlike typical ORs, some overrepresented ORs show robust cell surface expression in heterologous cells without the co-expression of RTPs. We present a model in which developing OSNs exhibit unstable OR expression until they choose to express an OR that exits the ER or undergo cell death. Our study sheds light on the new link between OR protein trafficking and OR transcriptional regulation.
https://doi.org/10.7554/eLife.21895.001eLife digest
Olfaction, or the sense of smell, is perhaps the most complicated and least understood of the five basic senses. Olfactory neurons in the nose can detect and distinguish between tens of thousands of different odor producing substances. They do so by using hundreds of unique sensors called olfactory receptors, each of which responds to a specific type of odor. During development, each olfactory neuron “chooses” to produce only one type of olfactory receptor. Once the neuron recognizes that the functional receptor is being generated and transported to the cell surface, it will stop making all the other olfactory receptors.
Chaperone proteins are responsible for transporting many olfactory receptors to the cell surface. To investigate how the loss of these chaperones affects how the olfactory system develops, Sharma et al. studied mice that were unable to produce the olfactory chaperone proteins. In these mice, developing neurons that chose to produce a type of olfactory receptor that depends on chaperone protein transport could not fully shut off other olfactory receptor genes. This led either to the neuron attempting to produce another type of receptor, or the death of the neuron. As a result, more neurons than usual produced receptors that do not require chaperone proteins to transport them to the cell surface. The olfactory neurons therefore produced only a fraction of all possible olfactory receptors, which decreased the ability of the mice to respond to odors.
In the future, it will be important to understand what determines whether an olfactory receptor can be transported to the cell membrane in the absence of chaperone proteins. Olfactory receptors are G protein-coupled receptors (GPCRs), which are the largest molecular class of drug targets for cancer and diseases that affect the brain and heart. Thus, results presented by Sharma et al. will also be relevant to researchers who study how GPCR malfunction causes diseases.
https://doi.org/10.7554/eLife.21895.002Introduction
Seven transmembrane G-protein coupled receptors (GPCRs), are diverse and the largest superfamily of receptors. Their roles are well established in sensing various stimuli including odorants, tastants, light, hormones, neurotransmitters and proteins. Some GPCRs require the presence of specific accessory proteins such as chaperones, vesicular targeting molecules and co-receptors for their cell surface expression (Lu et al., 2003; Salahpour et al., 2004; Dey and Matsunami, 2011; Wu et al., 2003). Mammalian olfactory receptors (ORs), which are GPCRs (Buck and Axel, 1991), are retained in the ER when expressed in non-olfactory cells. RTP1 (Receptor Transporting Protein 1) and RTP2 (Receptor Transporting Protein 2), both single transmembrane proteins strongly and exclusively expressed in the peripheral olfactory organs (Lu et al., 2003; Saito et al., 2004; Zhuang and Matsunami, 2008; Gimelbrant et al., 1999), greatly enhance the trafficking of ORs to the cell surface of heterologous cells. However, the role played by the RTPs in vivo remains unclear.
The mouse genome encodes over one thousand intact OR genes (Niimura et al., 2014), which are expressed in a singular and monoallelic fashion in each olfactory sensory neuron (OSN) (Shykind et al., 2004; Chess et al., 1994; Serizawa et al., 2000; Malnic et al., 1999). Each OR is not chosen at random; rather, OSNs express different ORs with dramatically varying probabilities (Khan et al., 2011; Ibarra-Soria et al., 2014).
OSNs in the olfactory epithelium (OE) are organized in overlapping zones defined by the expression of each OR (Ibarra-Soria et al., 2014; Ressler et al., 1993; Vassar et al., 1993; Miyamichi et al., 2005; Kanageswaran et al., 2015; Saraiva et al., 2015) as well as in a pseudostratified manner with progenitor cells forming the basal layer and mature neurons forming the upper layers. Mature OSN dendrites project into the nasal cavity forming a dendritic knob at the surface of the OE where they express ORs to interact with odorant molecules. Mature OSN axons expressing the same OR project to the olfactory bulb (OB) to converge onto specific glomeruli (Mombaerts et al., 1996; Ressler et al., 1994; Vassar et al., 1994; Hayar et al., 2004; Aungst et al., 2003; Gire et al., 2012). When the β2 Adrenergic Receptor (β2AR) is expressed instead of an OR, the β2AR –expressing OSNs target their axons to the OB and form glomeruli (Feinstein et al., 2004a, 2004b; Omura et al., 2014; Nakashima et al., 2013). Hence, the development of the peripheral olfactory system is dependent on functional GPCRs.
The mechanisms by which an OSN makes an OR choice have not been fully elucidated. Locus control-region like enhancers scattered on the genome and relative location of ORs from these elements have important roles in determining the probabilities of OR gene choice (Khan et al., 2011; Serizawa et al., 2003; Markenscoff-Papadimitriou et al., 2014). Various epigenetic mechanisms, for example histone modification by lysine demethylase (LSD1), allow the escape of an OR gene from repression (Magklara et al., 2011; Lyons et al., 2013, 2014; Armelin-Correa et al., 2014; Kilinc et al., 2016). OR expression is unstable until one OR is functionally expressed which then represses the expression of other OR alleles via negative feedback signaling through the unfolded protein response (UPR) and G proteins (Serizawa et al., 2003; Dalton et al., 2013; Wang et al., 2012; Li and Matsunami, 2013; Lewcock and Reed, 2004; Ferreira et al., 2014; Nguyen et al., 2007; Abdus-Saboor et al., 2016).
Here, we generated Rtp1 and Rtp2 double knockout mice (RTP1,2DKO) to investigate their role in the functioning and development of the olfactory system in vivo. We show that the RTP1,2DKO have OR trafficking defects, a substantial reduction in the number of mature OSNs, and an overall diminished olfactory capacity. Unexpectedly, we found that some ORs are overrepresented (referred to as oORs) while others are underrepresented (referred to as uORs) in RTP1,2DKO. Cells expressing a uOR lack stable gene choice in the mutant compared to wild-types while cells expressing an oOR do not show this instability, a result that links OR protein trafficking and OR transcriptional regulation.
Results
Generation of RTP1,2DKO mice
In order to study the role played by RTP1 and RTP2 in regulating OR expression and trafficking in vivo, we consecutively knocked out these genes while the intervening ~500 kb genomic region was not disrupted in ES cells (Figure 1A). Following chimeric mice production and germline transmission, we established mouse lines with Rtp1 and Rtp2 double knock out alleles. We found no phenotypic difference between Rtp1(+/+);Rtp2(+/+) (wild-type) and Rtp1(+/−);Rtp2(+/−) (het) mice. The Rtp1(−/−);Rtp2(−/−) homozygous mutants (RTP1,2DKO) showed no gross defects outside the olfactory system. Heterozygous crosses gave rise to wild-type, heterozygous and homozygous adults in roughly 1:2:1 ratio (Rtp1(+/+);Rtp2(+/+) 19, Rtp1(+/−);Rtp2(+/−) 34, Rtp1(−/−);Rtp2(−/−) 15, n = 10 mating pairs), suggesting no embryonic or postnatal lethality. We validated the absence of RTP1 and RTP2 transcripts in RTP1,2DKO by performing RNA in situ hybridization (Figure 1B).
RTP1,2DKO OSNs show defects in OR trafficking
It has been previously shown that RTP1 and RTP2 promote cell surface expression of ORs in the heterologous expression assays (Saito et al., 2004; Zhuang and Matsunami, 2007). Therefore, we used M71-IRES-tauGFP mice in which Olfr151 (also known as M71 and MOR171-2) expressing OSNs co-express tauGFP to examine the OSNs for OR trafficking defects (Feinstein et al., 2004a) (Figure 1C). In the RTP1,2DKO;M71-IRES-tauGFP OE, GFP staining was observed in the dendrites of Olfr151 positive OSNs (Figure 1D), indicating that the morphology of their OSNs remains unchanged. In contrast, immunostaining against Olfr151 (Barnea et al., 2004) was restricted to the cell body, indicating these OSNs are unable to traffic the OR to the dendrite (Figure 1D). Altogether, the data suggest that RTP1 and RTP2 are essential for OR trafficking.
RTP1,2DKO mice have fewer mature sensory neurons
Upon examination of the OE, we found that its thickness was significantly reduced in RTP1,2DKO mice. (p=0.02 paired student t test) (Figure 2A). We therefore examined the expression of various OSN developmental markers and signaling molecules in the OE to evaluate areas occupied by mature and immature OSNs in RTP1,2DKO. We compared OMP and adenylate cyclase 3 (ACIII), markers for mature neurons (Carter et al., 2004; Rogers et al., 1987), in 21 day old RTP1,2DKO mice and their littermates (Figure 2B, Figure 2—figure supplement 1). We measured the area occupied by RNA in situ hybridization signal against OMP and found that mice showed an average of 22% reduction in RTP1,2DKO when compared to the wild-type (p<0.0001, paired student t test, wild-type mean area 71% ± 5 (SD), RTP1, 2DKO mean area 49% ± 4 (SD)) (Figure 2C, See methods for details). Comparison of the OMP positive layer from wild-type and RTP1,2DKO OE collected at 1-day-old, 21-day-old and 6-month-old mice showed a significant reduction in OMP expression at 1 day and 21 days (p=0.0003, Mann Whitney U test, p=0.0003, Mann Whitney U test) but not at 6 months (Figure 2D). Immunohistochemical analysis of expression of adenylate cyclase 3 (ACIII), a signaling molecule expressed in mature OSNs (Wei et al., 1998; Wong et al., 2000; Col et al., 2007) showed a 17% decrease in the area occupied by the staining in 21 day old RTP1,2DKO OE (p=0.0001, Mann Whitney U test, wild-type mean area 44% ± 9 (SD), RTP1,2DKO mean area 27%, ± 3 (SD)) (Figure 2B). Consistent with OMP expression, we observed a significant difference in ACIII expression at 1 day old (p=0.0079, Mann Whitney U test) but not at 6 months (Figure 2E).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.21895.005
GAP43 is a marker for immature olfactory neurons in the OE (Meiri et al., 1988; Verhaagen et al., 1990; Treloar et al., 1999) and area occupied by it shows a 7% increase in the area of the OE it occupies in 21-day-old RTP1,2DKO (p=0.03 Mann Whitney U test, wild-type mean area 20%, ± 5 (SD), RTP1,2DKO mean area 27%, ± 6 (SD)) (Figure 2B). No significant difference in the GAP43 positive layer is observed between RTP1,2DKO and their littermates at 1 day nor at 6 months (Figure 2F).
Odorant evoked electrophysiological responses in RTP1,2DKO mice are diminished
Upon observation of fewer OSNs in RTP1,2DKO mice and lack of OR trafficking to the cilia (Figure 1D), we sought to test the olfactory ability by electroolfactogram (EOG). We tested a diverse set of 7 odorants, in both wild-type and RTP1,2DKO littermates. Wild-type mice show robust EOG responses to all odorants at concentrations as low as 0.01% (Figure 3A). In contrast, RTP1,2DKO mice showed striking deficits in their response. Responses to most odors were identical to the blank stimulus (air only), although some sensitivity was maintained for a subset of odorants (2-heptanone, amyl acetate, isomenthone) compared to the wild-types (Figure 3B–C). Thus, the reduction of mature OSNs and the loss of surface OR expression corresponds to a dramatic loss of odorant sensitivity in RTP1,2DKO.
OR expression is biased in RTP1,2DKO
To obtain a comprehensive view of gene expression changes in RTP1,2DKO, we performed an RNA-Seq on isolated whole olfactory mucosa including the OE and surrounding tissues. Differential expression analysis comparing RTP1,2DKO to wild-type littermates revealed that 3.8% of all genes (926/24,661) were differentially expressed between the two genotypes, among which 805 were downregulated and 121 were upregulated in RTP1,2DKO (FDR corrected p<0.05, see Experimental Procedures for details). Canonical signaling molecules known to be expressed in mature OSNs including Gnal (Gαolf), Adcy3 (ACIII), and Cnga2 were less abundant in the RTP1,2DKO consistent with a reduced number of mature OSNs in absence of RTP1 and RTP2. We found no significant difference in the expression levels of housekeeping genes like Gapdh and β actin (Supplementary file 1) (Kouadjo et al., 2007), neither did we see any compensatory increase in other RTP family members Rtp3 or Rtp4.
We then asked whether the loss of RTP1 and RTP2 equally affected all ORs. In a comparison between wild-type and RTP1,2DKO we found that 62% of intact ORs (678/1088) were significantly affected by the loss of RTP1 and RTP2 (Figure 4A). Close to half of the annotated intact ORs (562/1088) were downregulated in RTP1,2DKO (FDR corrected p<0.05), consistent with fewer OSNs in the mutant. Unexpectedly however, a small subset of OR transcripts (116/1088) were upregulated in RTP1,2DKO mice (FDR corrected p<0.05) (Figure 4B, Figure 4—figure supplement 1(A-B)).
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.21895.009
The disparity in the abundance of transcripts for these ORs raised the possibility of a difference in probabilities of OSNs expressing each OR in RTP1,2DKO. To remove any possible confounding variables from non-OSN cells, we normalized our read counts using only reads mapped on intact ORs and found that 531/1088 were underrepresented and 202/1088 were overrepresented (FDR corrected p<0.05) (Figure 4C,F, Figure 4—figure supplement 1C).
To further validate changes in numbers of OSNs expressing individual ORs in RTP1,2DKO mice, we carried out RNA in situ hybridization with probes against either underrepresented ORs (uORs) (Figure 4C) or overrepresented ORs (oORs) (Figure 4F). For all uORs tested, fewer OSNs were positive in RTP1,2DKO (Figure 4D–E) (p<0.05 Mann-Whitney U test). In stark contrast, the frequencies for the tested oORs were greater in RTP1,2DKO (p<0.05, Mann-Whitney U test) (Figure 4G–H). These results demonstrate that OR gene choice is biased in RTP1,2DKO towards a specific subset of receptor types.
Curiously, we found that oORs as a group are more abundantly expressed than uORs in the wild-type. The OR genes that were not classified as either underrepresented nor overrepresented (NS, not significant) exhibited a wide range of changes in expression levels between the wild-type and RTP1,2DKO, but are expressed at significantly lower abundance levels than both oORs and uORs (Figure 4I) (p<0.0001 one-way ANOVA, Tukey’s post hoc test).
The proportion of OSNs expressing oORs increases in older RTP1,2DKO mice
We wondered what happens to the proportion of OSNs expressing uORs and oORs in RTP1,2DKO mice at different ages. We performed RNA in situ hybridization with a (1) a probe mix containing 11 uORs and (2) a probe mix containing 25 oORs, all expressed in the dorsal region of the OE on 1-day, 21-day and 6-month-old OE. In the case of the uOR mix, wild-type showed an increase for OSNs expressing the uORs we tested both at 21 days and 6 months (Figure 5A–B) consistent to the increasing proportion of mature OSNs in the OE indicated by larger OMP and ACIII layers (Figure 2D–E). However, RTP1,2DKO showed no obvious increase in the fraction of cells expressing these ORs with age, while on the other hand, in RTP1,2DKO the number of neurons expressing oORs showed a dramatic increase both from 1 day old to 21 days and from 21 days to 6 months (p<0.0001 one-way ANOVA, Tukey’s post hoc test) (Figure 5C–D) demonstrating that the RTP1,2DKO OE is progressively populated by oORs.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.21895.012
Expression bias of an OR in RTP1,2DKO depends upon the protein sequence and not on the genomic location
We wondered whether the OR expression bias arose due to the effect of RTP on the regulatory elements of an OR’s gene locus or the protein. In our initial investigation, we did not find an obvious pattern or clustering for the genomic locations nor did we find obvious conserved residues or motifs amongst uORs or oORs (Figure 6A, Supplementary file 2). In an attempt to causally identify the basis of the bias, we used a mouse expressing β2AR-IRES-LacZ from the Olfr151 locus (Feinstein et al., 2004b) (Figure 6B) and asked whether the numbers of OSNs expressing Olfr151 or β2AR are similarly affected in RTP1,2DKO. We chose Olfr151 as it is a uOR and β2AR because it is a non-OR GPCR, capable of reaching the cell surface without the RTPs in heterologous culture and can replace a functional OR in native OSNs (Omura et al., 2014). We saw that fewer Olfr151 expressing OSNs were present in RTP1,2DKO (p=0.0028, Mann-Whitney U test), as expected for a uOR. Strikingly, more β2AR positive OSNs were present in RTP1,2DKO compared with wild-type (p<0.001, Mann-Whitney U test) (Figure 6C–D), suggesting that it is the protein sequence and not the locus of the OR that determines whether a given OR is underrepresented or overrepresented.
-
Figure 6—source data 1
- https://doi.org/10.7554/eLife.21895.014
Given that β2AR is known to be efficiently trafficked to the cell surface when heterologously expressed in the absence of the RTPs, we asked whether uORs and oORs show differential capabilities in cell surface trafficking in heterologous cells. We carried out live cell surface staining of HEK293T cells transfected with either Rho tagged uORs or oORs in the absence of RTP1 and RTP2. In order to quantify the surface staining, we carried out FACS to measure the surface OR levels. oORs as a group showed more OR surface expression than uORs (p<0.05, Mann-Whitney U test) (Figure 6E–G). Notably, the ORs that show most robust cell surface expression were all oORs. Even though trafficking mechanisms between OSNs and HEK293T cells are likely to be different, our data are consistent with the idea that RTP-independent trafficking of ORs may be related to increased frequencies of OSNs expressing oORs.
OSNs expressing oORs mature and function
Our results thus far suggest that OSNs expressing oORs are able to function despite the loss of the RTPs. In order to test this, we examined whether OSNs expressing uORs or oORs co-express OMP (Figure 7A). We found that the number of immature OMP-negative OSNs expressing uORs are similar in RTP1,2DKO and het controls, whereas the number of OMP-positive OSNs expressing uORs show a 69% decrease in RTP1,2DKO. (p=0.024 Mann Whitney U test) (Figure 7B–D).
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.21895.016
-
Figure 7—source data 2
- https://doi.org/10.7554/eLife.21895.017
To evaluate the function of OSNs expressing uORs or oORs, we chose a uOR and an oORs that have been previously deorphanized. Olfr1395 is an oOR found to respond to 2,5-dihydro-2,4,5-trimethylthiazoline (TMT) and Olfr923, a uOR, to acetophenone in vivo (Jiang et al., 2015). Using the induction of phospho ribosomal protein S6 (pS6), a marker for neuronal activation (Knight et al., 2012), we found that OSNs expressing Olfr1395 in both het and RTP1,2DKO were activated to its cognate ligand TMT (p<0.0001 one-way ANOVA, Tukey’s post hoc test) (Figure 7E,G). In contrast OSNs expressing Olfr923 in het, but not RTP1,2DKO were activated by their cognate ligand acetophenone (Figure 7E–F). These data show that OSNs expressing oORs mature and function in the RTP1,2DKO.
Immunostaining against active caspase 3, a cell death marker, suggested that RTP1,2DKO mice have an increased number of OSNs undergoing cell death (p<0.01, Mann-Whitney U test) (Cowan et al., 2001) (Figure 7H–I). Active caspase three staining in conjunction with OMP suggested that more OSNs in immature and mature layers both undergo cell death in RTP1,2DKO (Figure 7H). Notably, immature OSNs in the OMP negative layer rarely undergo cell death in wild-type (Figure 7J) (Jia et al., 2010). Together, these observations are consistent with the idea that OSNs expressing uORs are more likely to undergo cell death in RTP1,2DKO.
Persistent expression of nATF5 is observed in RTP1,2DKO OSNs expressing uORs
The above findings reveal the expression of a biased OR repertoire in RTP1,2DKO, with about a half of the ORs being underrepresented but a small subset of ORs overrepresented. Previous studies have shown that the unfolded protein response (UPR) plays an important role in OR gene choice mechanisms. During the initial phase of OR expression, the UPR pathway gives rise to an increased translation of nuclear activating transcription factor (nATF5) over an upstream inhibitory ORF via eIF2α signaling (Godin et al., 2016). Once the OSN has matured, OR gene choice is stabilized and the UPR is relieved (Dalton et al., 2013). We hypothesized that the lack of RTP1 and RTP2 causes persistent UPR in OSNs expressing uORs, leading to unstable OR gene choice in these OSNs, which in turn contributes to the skewed OR repertoire in mice lacking RTP1 and RTP2. We first asked whether the expression of nATF5, a marker for UPR, is different in RTP1,2DKO. Indeed, we observed that there were more nATF5 positive OSNs in RTP1,2DKO mice (p<0.001, Mann-Whitney U Test) and some of these OSNs were located closer to the apical surface, a phenomenon that was not observed in the wild-types, suggesting that nATF5 expression persists during the OSN development in RTP1,2DKO animals (Figure 8A–B). Similarly, expanded expression was observed for LSD1, a histone demethylase whose expression depends on nATF5 (Figure 8C). Next, we asked whether the ectopic expression of nATF5 in RTP1,2DKO is due to their protein sequence or the gene locus. We first compared the co-expression of nATF5 and Olfr151 between the wild-type and RTP1,2DKO using M71-IRES-tauGFP mice, and then asked whether the same co-localization occurs for β2AR expressed from the Olfr151 gene locus using β2AR-IRES-LacZ mice. We observed that the number of OSNs co-expressing Olfr151 and nATF5 was significantly higher in RTP1,2DKO than in the wild-type (p<0.05, Fisher’s exact test). In contrast, the number of OSNs co-expressing β2AR and nATF5 in the β2AR:IRES:tauLacZ mouse was not different from the wild-type (p=1, Fisher’s exact test)(Figure 8D–F). As Olfr151 but not β2ARs require RTPs for surface expression, these data suggest that delivery of ORs to the membrane plays a role in terminating the UPR.
-
Figure 8—source data 1
- https://doi.org/10.7554/eLife.21895.019
We next asked whether the expression of nATF5 in RTP1,2DKO mice was different in OSNs expressing uORs or oORs. In order to answer this question, we carried out fluorescent in situ hybridization against 7 uORs and 6 oORs along with immunohistochemistry for nATF5 and quantified colocalization of the OR signal with nATF5 staining. As expected, higher numbers of OSNs expressing uORs colocalized with nATF5 in RTP1,2DKO as a group (p=0.007, Mann-Whitney U test), whereas no significant difference was observed when oORs were tested (p=0.937, Mann-Whitney U test) (Figure 8G–H). These observations suggest that OSNs expressing uORs contribute to the expanded expression of nATF5 in RTP1,2DKO and are further consistent with the idea that surface trafficking of ORs is linked to turning off UPR.
OR gene expression is unstable in OSNs expressing Olfr151 in RTP1,2DKO
Increased nATF5 levels in uOR-expressing OSNs of RTP1,2DKO mice suggests a lack of stable OR gene choice in these neurons. This led us to hypothesize that OSNs that initially express a uOR may later turn off the OR and stabilize the expression of another OR. To directly address this, we used a lineage tracing strategy to study the stability of OR gene choice in RTP1,2DKO (Shykind et al., 2004; Dalton et al., 2013; Abdus-Saboor et al., 2016). In our assay (Figure 9A), we crossed a mouse that had one allele expressing Cre recombinase under the Olfr151 promoter, M71-IRES-Cre, to a mouse that had the Cre inducible fluorescent reporter Rosa26-lox-stop-lox-tdTomato (Madisen et al., 2010). In this mouse, any OSN expressing the M71-IRES-Cre allele at any point in time, would give rise to the permanent expression of tdTomato even if the OSN went on to express a different OR gene. We counted the number of tdTomato positive neurons that were also Olfr151 positive (double positive), which reflect the OSNs that initially chose as well as stably express Olfr151. OSNs that were tdTomato positive but not Olfr151 positive are the ones that switched their initial gene choice. We found that only 17% of tdTomato positive cells from RTP1,2DKO OE (n = 24/140 neurons) also express Olfr151 in comparison to 31% OSNs from the wild-type (n = 79/258 neurons) (p<0.05, Fisher’s exact test) and 38% (n = 66/172 neurons) in the heterozygotes (p<0.05, Fisher’s exact test) (Figure 9B–C). This suggests that the loss of RTPs leads to the frequent termination of Olfr151 gene expression.
To determine whether the Olfr151 gene switched to an oOR within the same locus we carried out a co-localization analysis between Olfr143, an oOR within the Olfr151 locus and tdTomato under the control of M71-cre (Figure 9—figure supplement 1). We found no Olfr143 and tdTomato double positive OSNs in both RTP1,2DKO or their wild-type littermates, indicating that there is no higher likelihood of the gene switching between these ORs.
Olfr151 glomeruli are not formed in RTP1,2DKO
Lastly, we investigated OSN axon targeting to the olfactory bulb (OB) in RTP1,2DKO. Olfactory axons entered the OB and innervate to the glomerular layer based on our OMP immunostaining in RTP1,2DKO (Figure 10A). However, using M71-IRES-tau GFP mice we found that Olfr151 expressing OSNs were unable to converge in the OB in RTP1,2DKO while their wild-type littermates had two Olfr151 glomeruli in each of their OBs as expected (Figure 10B–C). To investigate whether the axon targeting defect was ubiquitous to all receptors, we used β2AR-IRES-tauLacZ (Feinstein et al., 2004b) (Figure 6B). Both the wild-type and RTP1,2DKO mice formed glomeruli, however the mutant mice had ectopic glomeruli for OSNs expressing this GPCR (Figure 10B–C). These data suggest the RTP1,2DKO mice do not have a complete set of glomeruli, but retain the ability to form them. We observed tdTomato-positive axons forming small glomeruli in RTP1,2DKO;M71-IRES-cre; Rosa26-lox-stop-lox-tdTomato mice (2 out of 3 mice examined had 1 glomerulus each)(Figure 10B–C). This may indicate that OSNs initially expressing Olfr151 switch and/or stabilize expression of a specific OR, presumably an oOR, and axons from these OSNs can converge in the OB.
Discussion
In the current study, we investigated the in vivo role of RTP1 and RTP2 in the mouse olfactory system by generating and analyzing RTP1,2DKO mice. Our results demonstrate a surprising link between receptor trafficking of ORs, the UPR and OR gene choice.
Differential control of OR representations by RTPs
The absence of RTP1 and RTP2 leads to the underrepresentation of nearly half of the ORs, while about 10% of the ORs are significantly overrepresented. This translates into a change in the numbers of OSNs choosing each OR. How do the RTPs regulate the probability of OR gene choice?
Our data indicate that protein sequences of ORs differentially influence OR gene choice, which is linked to trafficking of ORs to the cell surface. Our attempts to identify protein motifs, domains, or features specific to either uORs or oORs were so far unsuccessful. This fits with recent reports where large-scale mutational analysis of Olfr151 in heterologous cells failed to identify any specific amino acids or domains that regulate its cell surface expression (Hague et al., 2004; Jamet et al., 2015). Nevertheless, this study provides a large set of sequence information of ORs that will allow us to conduct future structure-function studies by testing uOR/oOR chimeras and/or searching for hidden features within uORs or oORs, which may in turn give us clues as to why nearly half of the ORs are retained in the ER in the absence of the RTPs.
Prolonged UPRs in OSNs expressing uORs in RTP1,2DKO
Previous studies have suggested the UPR protein nATF5 is induced once an OSN starts actively expressing an OR and is lost when OR expression is stabilized (Dalton et al., 2013). Our results show an expanded expression of nATF5 in RTP1,2DKO OSNs, suggesting persistent UPR during OR gene choice. Importantly, the frequency of co-localization between nATF5 and uORs increases in RTP1,2DKO whereas this increase is not observed for oORs. This suggests that the persistent UPR phenotype observed in RTP1,2DKO is due to the OSNs that express uORs resulting in unstable OR gene choice for these OSNs. This is reinforced by our observation that OSNs initially expressing Olfr151, a uOR, are more likely to terminate their expression in RTP1,2DKO mice. We present a model for the role played by RTP1 and RTP2 in the gene choice made by an OSN (Figure 11A). In our model, the RTPs suppress UPR response by allowing ORs to exit the ER and be transported to the plasma membrane. uORs are not trafficked to the cell surface in the absence of the RTPs, giving rise to persistent UPR in the OSNs that express them; as a consequence, destabilization of the initial OR gene choice leads to cell death or to the stabilization of another OR. In contrast, oORs are trafficked to the cell surface as functional proteins in the absence of the RTPs, allowing these OSNs to terminate UPR and stabilize gene expression. Although we cannot rule out that RTPs themselves play a role in the elimination of the UPR response, the lack of increase in nATF5 observed in OSNs expressing oORs in RTP1,2DKO suggests interaction between the RTPs and the UPR pathway through ORs. Even for OSNs expressing oORs and β2AR, UPR is likely to be induced in the initial stage (Dalton et al., 2013). Recent reports utilizing single cell RNA-Seq suggest that immature OSNs express as many as 12 ORs at low levels but mature OSNs only show the expression of one dominant OR (Saraiva et al., 2015; Tan et al., 2015; Hanchate et al., 2015). Low-level expression of uORs precedes the expression of the RTPs and may be sufficient to trigger UPR in the developing OSNs. Alternatively, both oORs and uORs induce UPR at the initial stage.
RNA-Seq data in the wild-type suggests that oORs as a group tend to be more frequently chosen. It could be that initial expression of RTP1 and RTP2 in the developing OSNs is not stable or abundant enough, causing oORs to be stabilized. Alternatively oORs tend to ‘win the competition’ (Nguyen et al., 2007; Abdus-Saboor et al., 2016) even in the presence of RTPs probably because of its ability to suppress UPR. Both change in probability of OR gene choice and biased cell death could alter OR population representation both in the level of transcripts and in the number of OSNs. The relative contributions of cell death and gene switching to the differential representation of ORs can be further clarified using Bax knockout mice where cell death in developing neurons is suppressed (Robinson et al., 2003).
OR genes that did not significantly change in RTP1,2DKO mice showed lower abundance in wild-type mice suggesting that these ORs are chosen less frequently. Deeper sequencing and/or increased sample sizes will help classify these ORs as underrepresented, overrepresented, or not changed.
In our lineage tracing experiment, we were unable to distinguish tdTomato-positive, Olfr151-negative OSNs between those express Olfr151 as one of the low-level ORs, and the others express Olfr151 as the dominant OR later switch to express another OR. Irrespective of this our data show greater gene instability for Olfr151 expression without the RTPs. Consistent with multiple OR gene expression in developing OSNs, our lineage tracing results suggest that only 31% of OSNs that initially expressed Olfr151 went on to stabilize its expression in the wild-type. Curiously, the vast majority of OSNs that initially expressed Olfr1507 (MOR28) show stable expression in a similar experiment (Shykind et al., 2004; Dalton et al., 2013), suggesting that some ORs are more likely to be stabilized than other ORs.
The frequency with which any OR is chosen is different and the underlying cause for this difference remain unknown. Our results show that the number of OSNs expressing β2AR expressed from the Olfr151 locus is lower than the number expressing Olfr151. Our data call for future experiments to test whether protein sequences of ORs differentially influence initial OR gene choice.
RTP1,2DKO mice show diminished but not abolished responses to odors
Our data suggest that mice without RTP1 and RTP2 had diminished but not abolished responses to odorants. Even though we see a clear reduction in the number of mature OSNs and dramatically diminished responses to odors in RTP1,2DKO mice, functional OSNs are not completely eliminated. Previous studies suggest that anosmic mice often die in the first few days after birth (Brunet et al., 1996), but the RTP1,2DKO mice seem to have the sufficient olfactory ability for postnatal development. Our data show that OSNs expressing oORs are likely to mature and be functional in RTP1,2DKO, explaining the residual responses observed in the mutant. These OSNs become more abundant in the mutant OE as the mouse ages indicating that these ORs help maintain the olfactory ability of these mice. It will be interesting to assess olfactory-mediated behaviors of the RTP1,2DKO mice in which only a minor fraction of the ORs are functional, since this could address a fundamental question of why most mammals have so many ORs.
No GFP-positive glomeruli were observed in RTP1,2DKO;M71-IRES-tauGFP indicating OSNs expressing Olfr151 are unable to converge their axons without RTP1 and RTP2. Yet we observed small tdTomato-positive glomeruli in RTP1,2DKO; M71-IRES-cre; Rosa26-lox-stop-lox-tdTomato. These tdTomato-positive glomeruli could be formed by OSNs that initially chose Olfr151 and then switched and/or stabilized the expression of the same OR, presumably an oOR. Previous reports have shown that OR gene switching tend to take place within the same gene locus (Roppolo et al., 2007; Pacifico et al., 2012) leading us to test the hypothesis that the tdTomato-positive axons forming these glomeruli in RTP1,2DKO mice could be stabilizing the expression of Olfr143, an oOR near the Olfr151 gene locus. However, we found no tdTomato-positive OSNs that also expressed Olfr143. This suggests that at least a portion of tdTomato positive OSNs stabilize the expression of the same OR, which is probably an oOR other than Olfr143. Further investigation of the identities of tdTomato-positive OSNs could further our understanding of the OR gene switching mechanism.
Functional ORs expressed outside olfactory system
A number of reports have shown that ORs are expressed in various organs outside the olfactory system (Griffin et al., 2009; Feldmesser et al., 2006; Flegel et al., 2013; Kang and Koo, 2012), whereas expression of RTP1 and RTP2 appears to be confined to the peripheral olfactory tissues. One well-established example of a functional OR outside the olfactory system is Olfr78 (also known as MOR18-2) and its human ortholog OR51E2, which have been reported to function in the prostate, airway and kidney as well as the carotid body (Aisenberg et al., 2016; Chang et al., 2015; Wang et al., 2006; Pluznick et al., 2013). This receptor is an oOR and is trafficked to the cell surface in heterologous cells without the RTPs (Zhou et al., 2016; Pluznick et al., 2011; Neuhaus et al., 2009). It is also interesting to note that ORs expressed in the bladder and thyroid (Olfr544, Olfr558, and Olfr1386) are all oORs (Kang et al., 2015). Together, it is tempting to speculate that oORs expressed outside the OE are more likely to play chemosensory roles (Feingold et al., 1999).
In conclusion, our study suggests the importance of OR trafficking by the RTP family members in OSN function and the probability of OR gene choice. Our studies contribute to a broader understanding of how cells discern the presence of a GPCR on their cell surface post translationally and link protein trafficking to epigenetic modifications that give rise to changes in the cell’s expression profile. In the future, it will be interesting to ask whether the UPR pathway that links the functional trafficking of receptors to the cell surface with epigenetic modifications is utilized by other tissue types.
Materials and methods
Generation of RTP1 and RTP2 double knockout mice
Request a detailed protocolThe strategies for generating RTP1 and RTP2 knockout mice are illustrated in Figure 1. The coding regions of RTP1 and RTP2 were replaced by loxP and puromycin (Pac) cassette, respectively. Fragments used for the left and right arms were amplified by PCR using BAC clone derived from C57BL/6 mice as templates. For RTP1, ACN cassette and DT-A cassette were used for positive and negative selection, respectively. A targeting vector was electroporated into the ES cell line EF-1, which is 129/B6 hybrid (Bronson et al., 1996). Colonies were picked up in G418-containing medium. Genomic DNA was digested by Acc65 I and hybridized with a 500 nt external probe on Southern blots. To generate RTP1&2 double knockouts, the RTP2 targeting vector was electroporated into two ES cell lines in which RTP1 is replaced by loxP-Cre-neo-loxP cassette. Colonies were picked up in puromycin -containing medium. Genomic DNA was digested by Bgl II and hybridized with a 500 nt external probe on Southern blots. Targeted ES cell clones were injected into C57BL/6 blastocysts. Chimeric mice were bred with C57BL/6 mice. Mice with RTP1,2 mutant allele were maintained by backcrossing with C57BL/6. Primers for genotyping are: 5’-cggaattcatgtcaggctgcaacttc-3’ and 5’-gggccgatattgggttaggag-3’ for WT allele; 5’-agccagctcttaagtccttac-3’ and 5’-gctcgagatctagatatcgataccgt-3’ for mutated allele. Primers for genotyping RTP2 KO are: 5’-ccctgaagagtctcacccgctc-3’ and 5’-cacatataccccaacttctagg-3’ for WT allele; 5’-caaacagacgaaccctagcaattcccactg-3’ and 5’-cttcattctcagtattgttttgccaagttc-3’ for mutated allele.
Animal strains
Request a detailed protocolThe procedures of animal handling and tissue harvesting were approved by the Institutional Animal Care and Use Committee of Duke University. All experiments were carried out on both male and female mice. The following mouse strains were all obtained from The Jackson Laboratory: Olfr151tm14(Adrb2)Mom/MomJ (β2 Adrenergic Receptor-IRES-LacZ) (Feinstein et al., 2004b), Stock no. 006691 (RRID:IMSR_JAX:006691). B6;129P2-Olfr151tm26Mom/MomJ (M71-IRES-tauGFP) (Feinstein et al., 2004a), Stock no. 006676 (RRID:IMSR_JAX:006676). Olfr151tm28(cre)Mom/MomJ (M71-IRES-Cre) (Li et al., 2004), Stock no. 006677 (RRID:IMSR_JAX:006677). B6;129S6-Gt(ROSA)26Sortm9(CAG-tdTomato)Hze/J (ROSA26-loxP-stop-loxP-tdTomato) (Madisen et al., 2010), Stock no. 007905 (RRID:IMSR_JAX:007905).
DNA constructs
Request a detailed protocolOR ORFs were cloned into pCI (Promega) with a Rho tag at the N terminal. All plasmids were verified using Sanger’s sequencing (3100 Genetic Analyzer, Applied Biosystems).
Electroolfactrograms (EOG)
Request a detailed protocolMice aged 2–4 months were sacrificed by anesthesia followed by rapid decapitation. After removing the skin, the skull was hemisected along the midline with a razor blade, exposing the turbinates. Hemisections were stabilized with pins in a custom Sylgard chamber for recording. Electroolfactograms were measured using glass pipettes filled with ACSF (tip size15–20µm, resistance ~0.5 MΩ). For some recordings, electrode tips were filled with 0.5% agarose. Electrodes were placed on the anterior surface of turbinate II and referenced to an Ag/AgCl wire placed on the surrounding bone. Signals were filtered (0.1Hz - 100 Hz) and amplified (1000X) using a differential amplifier (DAM-80, WPI). High-pass filtering introduced a slight rebound above baseline for strong responses. Odors were delivered using a custom olfactometer at a final dilution of 0.01% (1:1000 vol/vol in mineral oil, and a further 1:10 via airflow by combining flow from odorant headspace with a moisturized deodorized airstream) and a total flow rate of 100 ml/min. Odorants were obtained from Sigma at the highest purity available and diluted in mineral oil. Epithelia were kept moisturized at all times. The blanks are an average of multiple interleaved trials interspersed within the series. This averaged blank is re-displayed for each different odorant for comparison.
Immunofluorescence
Request a detailed protocolMice were weaned at three weeks and their OE was dissected out and flash frozen. 20–25 micron thick sections were cut onto slides which were stored at −80°C and were subsequently thawed and fixed in 4% PFA for 20 min, washed for 1 min with 0.5% triton-x-100 in PBS and then rinsed twice in PBS. They were blocked for 1 hr (See table for blocking reagent and antibody concentration) in PBS with 0.1% triton-x-100 and then kept overnight at 4°C in primary antibody solution made in the same blocking reagent. The antibody solution was washed 3 times for 5 min in PBS and then subject to the secondary antibody for an hour. 0.001% bisbenzimide, used to visualize the nucleus, was added to the slides for 1 min followed by 3 × 5 min washes in PBS. The coverslip was mounted in 5–6 drops of Mowiol.
No. | Target | Source | Dilution | Company | Blocking | Catalog no. | RRID |
---|---|---|---|---|---|---|---|
1 | GFP | Rabbit | 1:400 | 5% milk | |||
2 | OMP | Goat | 1:400 | Wako | 5% milk | 019-22291 | AB_664696 |
3 | pS6 | Rabbit | 1:200 | ThermoFisher | 5% milk | 44-923G | AB_2533798 |
4 | Caspase 3 | Rabbit | 1:1000 | Cell Signaling | 5% milk | 9661 | AB_2341188 |
5 | LSD1 | Rabbit | 1:800 | abcam | 4% Donkey serum | ab17721 | AB_443964 |
6 | ATF5 | Goat | 1:1000 | Santa Cruz | 4% Donkey serum | sc-46934 | AB_2058761 |
7 | M71 | Guinea Pig | 1:3000 | Barnea et al. | 4% Donkey serum | ||
8 | LacZ | Mouse | 1:3000 | Promega | 4% Donkey serum | Z3781 | AB_430877 |
9 | Td Tomato | Rabbit | 1:10,000 | Rockland | 4% Donkey serum | 600-401-379 | AB_2209751 |
RNA in situ hybridization
Request a detailed protocolCandidate ORs which had less than 80% identity with all other ORs were chosen from the RNA-Seq data. Probes against their ORFs were prepared by the addition of the T3 start site to the 3’ end of the pCI plasmid primer followed by incorporation of digoxigenin (DIG) using T3 RNA polymerase (Promega) and alkaline degradation to get labeled probes of around 200 bp. Slides were prepared as described above and fixed in 4% PFA for 15 min and washed twice in PBS. They were acetylated in 1.2% triethanolamine and dropwise addition of acetic acid. They were washed in PBS for 5 min and prehybridized in buffer for an hour in large plates soaked in 5XSSC and 50% formamide at 58°C. 1 μl of the aforementioned labeled probe in 200 μl of the prehybridization buffer was pipetted on to the slides and covered with a layer of parafilm and kept overnight at 58°C. The slides were thoroughly washed in 5XSSC and then in 0.2XSSC twice for 30 min, 5 min in PBS and then blocked for an hour. The slides were stained with alkaline phosphatase conjugated antibody against DIG (Roche) for an hour and then kept in development buffer for 5 min before being subject to NBT-BCIP in development buffer overnight. Slides were then stained with bisbenzimide and mounted in Mowiol as described above.
Fluorescent in situ hybridization (FISH)
Request a detailed protocolSlides were prepared and treated using the in situ hybridization protocol until the antibody staining step. Horse radish peroxidase (HRP) conjugated antibody against DIG was used instead of the alkaline phosphatase conjugate. Antibody staining amplification was carried out using tyramide signal amplification (TSA) using cy3 as the flourophore (PerkinElmer) followed by antibody staining against nATF5 as described in the immunofluorescence section above.
Whole mount LacZ staining
Request a detailed protocolThree-week-old mice were sacrificed and the entire head was dissected and kept in 4%PFA for 30 min. The tissue was washed in buffer A (100 mM phosphate buffer [pH 7.4], 2 mM MgCl2, and 5 mM EGTA) once for 30 min and then for 5 min in buffer A and buffer B (100 mM phosphate buffer [pH7.4], 2 mM MgCl2, 0.01% sodium desoxycholate, and 0.02% Nonidet P40) and then kept in buffer C (buffer B, with 5 mM potassium-ferricyanide, 5 mM potassium-ferrocyanide, and 1 mg/ml of X-Gal) containing X gal overnight at 4°C. Whole mount tissue was then washed in PBS and imaged under a 5x objective (Mombaerts et al., 1996). The medial glomeruli were not observed using this method.
RNA-Seq
Request a detailed protocolThe whole olfactory mucosa was dissected out of 3-week-old mice (2 males and 1 female sex matched littermates) and homogenized in trizol. This solution was centrifuged at maximum speed for 10 min and the supernatant collected was treated with 0.2 ml chloroform for every 1 ml of trizol, hand shaken for 3 min and spun for 15 min at maximum speed. The aqueous phase was collected and added to isopropanol (500 μl per 1 ml of trizol) and shaken and kept for 5 min at room temperature and then spun for 10 min. The RNA pellet thus collected was washed once in 75% ethanol by adding 500 μl per 1 ml of trizol and centrifuged for 2 min and then again in 180 μl of 75% ethanol. The pellet was then briefly air-dried before being suspended in 50 μl of water and the concentration of RNA was determined using a spectrophotometer by taking 1 μl of the sample and diluting it with 9 μl of water. RNeasy cleanup kit was then used to process the sample. Library generation was carried out using Illumina TruSeq Stranded RNA-Seq kit and HiSeq Illumina sequencing was carried out at the Duke Sequencing and Genomic Technologies Core.
Reads were mapped, using kallisto (Bray et al., 2016), to the mouse transcripts which were downloaded from UCSC Genome Browser (https://genome.ucsc.edu) and whose OR genes were replaced with extended OR gene annotations. Reads assigned on each gene were counted by kallisto. The read count table thus generated was analyzed using EdgeR. DESeq was used to calculate the size factors of individual libraries, FDR (False Discovery Rate) was used to adjust multiple comparisons between the ORs as previously published (Jiang et al., 2015).
Cell culture
Request a detailed protocolHEK293T cells were grown in Minimal Essential Medium (MEM) containing 10% FBS (vol/vol) with penicillin-streptomycin and amphotericin B at 37°C, saturating humidity and 5% CO2. These cells were authenticated using polymorphic short tandem repeat (STR) at the Duke DNA Analysis Facility using GenePrint 10 (Promega) and shown to share profiles with the reference (ATCC). No mycoplasma infection was detected.
FACS
Request a detailed protocolHEK293T cells were grown to a 100% confluence before being trypsinized and resuspended and seeded onto 35 mm plates at 25% confluency. These plates were cultured overnight then transfected with rho tagged ORs in the plasmid pCI for 18 to 20 hr along with GFP expression vector to monitor the transfection efficiency. The cells were stripped with cell stripper and triturated before being kept in 5 mL round bottom polystyrene (PS) tubes (Falcon 2052) on ice. The cells were spun down at 4°C and resuspended in PBS containing 10 mM HEPES, 15 mM NaN3, and 2% FBS to wash the cell stripper. They were subjected to 30 min in primary antibody (mouse anti Rho [Laird and Molday, 1988]) and then washed, stained with phycoerythrin (PE)-conjugated donkey anti-mouse antibody (Jackson Immunologicals) in the dark. 7-Amino-actinomycin D (Calbiochem), a fluorescent, cell-impermeant DNA binding agent that selectively stains dead cells, was added to eliminate dead cells. The cells were analyzed using BD FACSCanto II FACS with gating allowing for 10,000 GFP positive, single, spherical, viable cells and the results were analyzed using Flowjo (Dey and Matsunami, 2011).
Image analysis
Request a detailed protocolAll images were captured on Zeiss Axioskop two fluorescent microscope using Q capture pro. The images were then analyzed using ImageJ. For OR bias, phosphorylated S6 and caspase three experiments, nuclear staining was quantified by selecting the OSN layer in the OE and using the maxima function in ImageJ to select and count all cells followed by quantification of OR or Caspase three positive OSNs by hand scoring using the cell counter function. Percent positive cells were calculated as (Positive cells/ total number of cells)*100. ATF5 was quantified by selecting the OSN layer and digitally straightening it using ImageJ followed by manually selecting the ATF5 positive cells and using the measure function to get the X and Y co-ordinates to measure the height from the basal end of the epithelium. All colocalization experiments were manually scored by selecting the OR positive cell’s nucleus and measuring the pixel intensity for the same selection. The pixel intensity of the neighboring area was subtracted to remove background and determine positive cells.
Marker gene analysis
Request a detailed protocolThe area occupied by OSNs in the OE was hand selected using nuclear staining in image J. The thickness of the sections was measured by straightening the OE followed by measuring the height of the straightened section in four places. The average height for every position was compared between the two genotypes. The OMP or ACIII or GAP43 positive area was selected using imageJ thresholding and these areas were measured and expressed as percentages in the figures. Images were taken at five roughly equivalent positions in the OE using the VNO and OB as landmarks (anterior VNO, middle VNO, posterior VNO, anterior OB and middle OB, Figure 2—figure supplement 1 shows an example of matching OE sections using the anterior OB). We further analyzed all the sections from RTP1,2DKO mice at 1 day, 21 day and 6 months and compared the difference in percent area occupied by the OMP positive layer, ACIII positive layer and GAP43 positive layer with that of their littermates. Each data point is obtained from one image from one matched section. The data sets were compared using a paired student t test or Mann Whitney U test as indicated in the figure legend.
Statistical analysis
Request a detailed protocolPercent positive cells were calculated by hand scoring positive cells and calculating percentage based on the total number of OSNs in the image counted using nuclear staining. Each individual section was counted as an individual data point. Height of positive cells was calculated by straightening OE sections and obtaining the Y co-ordinates of hand scored ATF5 positive cells. Multiple comparison data from our ANOVA analysis is included in supplementary file 3.
pS6 staining quantification
Request a detailed protocolEach OSN positive for OR FISH signal was selected in image J. These selections were used to measure the pixel intensity of pS6 staining. The average pixel intensity of the entire OE selection was used as background. The average intensity of each cell was normalized by the background followed by its subtraction.
Data availability
-
Olfactory receptor accessory proteins play crucial roles in receptor gene choicePublicly available at the NCBI Gene Expression Omnibus (accession no: GSE87694).
References
-
Defining an olfactory receptor function in airway smooth muscle cellsScientific Reports 6:38231.https://doi.org/10.1038/srep38231
-
Near-optimal probabilistic RNA-seq quantificationNature Biotechnology 34:525–527.https://doi.org/10.1038/nbt.3519
-
Olfactory horizontal basal cells demonstrate a conserved multipotent progenitor phenotypeJournal of Neuroscience 24:5670–5683.https://doi.org/10.1523/JNEUROSCI.0330-04.2004
-
Caspases 3 and 9 send a pro-apoptotic signal from synapse to cell body in olfactory receptor neuronsJournal of Neuroscience 21:7099–7109.
-
Truncation releases olfactory receptors from the endoplasmic reticulum of heterologous cellsJournal of Neurochemistry 72:2301–2311.https://doi.org/10.1046/j.1471-4159.1999.0722301.x
-
Mitral cells in the olfactory bulb are mainly excited through a multistep signaling pathJournal of Neuroscience 32:2964–2975.https://doi.org/10.1523/JNEUROSCI.5580-11.2012
-
Emerging Roles for the Unfolded Protein Response in the Developing Nervous SystemTrends in Neurosciences 39:394–404.https://doi.org/10.1016/j.tins.2016.04.002
-
External tufted cells: a major excitatory element that coordinates glomerular activityJournal of Neuroscience 24:6676–6685.https://doi.org/10.1523/JNEUROSCI.1367-04.2004
-
The olfactory transcriptomes of micePLoS Genetics 10:e1004593.https://doi.org/10.1371/journal.pgen.1004593
-
Molecular profiling of activated olfactory neurons identifies odorant receptors for odors in vivoNature Neuroscience 18:1446–1454.https://doi.org/10.1038/nn.4104
-
Evidence against the role of rhodopsin in rod outer segment binding to RPE cellsInvestigative Ophthalmology & Visual Science 29:419–428.
-
Unfolding the mystery of olfactory receptor gene expressionDevelopmental Cell 27:128–129.https://doi.org/10.1016/j.devcel.2013.10.012
-
Distribution and phosphorylation of the growth-associated protein GAP-43 in regenerating sympathetic neurons in cultureJournal of Neuroscience 8:2571–2581.
-
Activation of an olfactory receptor inhibits proliferation of prostate cancer cellsJournal of Biological Chemistry 284:16218–16225.https://doi.org/10.1074/jbc.M109.012096
-
The β2-adrenergic receptor as a surrogate odorant receptor in mouse olfactory sensory neuronsMolecular and Cellular Neuroscience 58:1–10.https://doi.org/10.1016/j.mcn.2013.10.010
-
Gene cluster lock after pheromone receptor gene choiceThe EMBO Journal 26:3423–3430.https://doi.org/10.1038/sj.emboj.7601782
-
Homodimerization of the beta2-adrenergic receptor as a prerequisite for cell surface targetingJournal of Biological Chemistry 279:33390–33397.https://doi.org/10.1074/jbc.M403363200
-
Mutually exclusive expression of odorant receptor transgenesNature Neuroscience 3:687–693.https://doi.org/10.1038/76641
-
Glomerular formation in the developing rat olfactory bulbThe Journal of Comparative Neurology 413:289–304.https://doi.org/10.1002/(SICI)1096-9861(19991018)413:2<289::AID-CNE9>3.0.CO;2-U
-
Single cell transcriptome analysis of mouse carotid body glomus cellsThe Journal of Physiology 594:4225–4251.https://doi.org/10.1113/JP271936
-
Synergism of accessory factors in functional expression of mammalian odorant receptorsJournal of Biological Chemistry 282:15284–15293.https://doi.org/10.1074/jbc.M700386200
Article and author information
Author details
Funding
National Institutes of Health (R01 DC014423)
- Hiroaki Matsunami
National Institutes of Health (R01 DC012095)
- Hiroaki Matsunami
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Gilad Barnea and Richard Axel for anti-M71 antibodies. Cheryl Bock and other members of Duke Cancer Institute Transgenic Mouse Facility for ES cell targeting and chimeric mouse production; Mengjue Jessica Ni for expert technical assistance. Simone Weyand, Ting Zhou, Claire de March, Xiaoyang Serene Hu, Tatjiana Abaffy, Kevin Zhu, Aashutosh Vihani, along with other members of the Matsunami lab and for critical input and discussion of the experiments and manuscript; William Wetsel for his help with statistical analysis; Doug Marchuk and Debby Silver for generously sharing equipment; Jianghai Ho and Mike Cook in Duke Flow Cytometry Core for help with FACS analysis and Duke Sequencing and Genomic Technologies Core for carrying out the RNA-Seq. This work was supported by grants from NIH.
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocols (A161-16-07) of the Duke Animal Care and Use program.
Copyright
© 2017, Sharma et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,333
- views
-
- 723
- downloads
-
- 31
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
Classical G-protein-coupled receptor (GPCR) signaling takes place in response to extracellular stimuli and involves receptors and heterotrimeric G proteins located at the plasma membrane. It has recently been established that GPCR signaling can also take place from intracellular membrane compartments, including endosomes that contain internalized receptors and ligands. While the mechanisms of GPCR endocytosis are well understood, it is not clear how well internalized receptors are supplied with G proteins. To address this gap, we use gene editing, confocal microscopy, and bioluminescence resonance energy transfer to study the distribution and trafficking of endogenous G proteins. We show here that constitutive endocytosis is sufficient to supply newly internalized endocytic vesicles with 20–30% of the G protein density found at the plasma membrane. We find that G proteins are present on early, late, and recycling endosomes, are abundant on lysosomes, but are virtually undetectable on the endoplasmic reticulum, mitochondria, and the medial-trans Golgi apparatus. Receptor activation does not change heterotrimer abundance on endosomes. Our findings provide a subcellular map of endogenous G protein distribution, suggest that G proteins may be partially excluded from nascent endocytic vesicles, and are likely to have implications for GPCR signaling from endosomes and other intracellular compartments.
-
- Cell Biology
Anionic lipid molecules, including phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2), are implicated in the regulation of epidermal growth factor receptor (EGFR). However, the role of the spatiotemporal dynamics of PI(4,5)P2 in the regulation of EGFR activity in living cells is not fully understood, as it is difficult to visualize the local lipid domains around EGFR. Here, we visualized both EGFR and PI(4,5)P2 nanodomains in the plasma membrane of HeLa cells using super-resolution single-molecule microscopy. The EGFR and PI(4,5)P2 nanodomains aggregated before stimulation with epidermal growth factor (EGF) through transient visits of EGFR to the PI(4,5)P2 nanodomains. The degree of coaggregation decreased after EGF stimulation and depended on phospholipase Cγ, the EGFR effector hydrolyzing PI(4,5)P2. Artificial reduction in the PI(4,5)P2 content of the plasma membrane reduced both the dimerization and autophosphorylation of EGFR after stimulation with EGF. Inhibition of PI(4,5)P2 hydrolysis after EGF stimulation decreased phosphorylation of EGFR-Thr654. Thus, EGFR kinase activity and the density of PI(4,5)P2 around EGFR molecules were found to be mutually regulated.