Different expression domains for two closely related amphibian TAARs generate a bimodal distribution similar to neuronal responses to amine odors

TL;DR

This study shows how two similar amphibian olfactory receptors have different expression patterns, leading to a bimodal response to amine odors.

Contribution

The paper reveals a bimodal distribution in amine odor responses linked to distinct expression domains of two TAARs in amphibians.

Findings

TAAR4b and TAAR4a show distinct spatial expression patterns in two dimensions.

The combined expression of TAAR4a and TAAR4b creates a bimodal distribution similar to neuronal amine responses.

Quantitative analysis suggests other receptors may also contribute to amine odor perception in amphibians.

Abstract

Olfactory perception is mediated by a multitude of olfactory receptors, whose expression in the sensory surface, the olfactory epithelium, is spatially regulated. A common theme is the segregation of different olfactory receptors in different expression domains, which in turn leads to corresponding segregation in the neuronal responses to different odor groups. The amphibian olfactory receptor gene family of trace amine associated receptors, in short TAARs, is exceedingly small and allows a comprehensive analysis of spatial expression patterns, as well as a comparison with neuronal responses to the expected ligands for this receptor family, amines. Here we report that TAAR4b exhibits a spatial expression pattern characteristically different in two dimensions from that of TAAR4a, its close homolog. Together, these two genes result in a bimodal distribution resembling that of amine responses as visualized by calcium imaging. A stringent quantitative analysis suggests the involvement of additional olfactory receptors in amphibian responses to amine odors.

Figures

Full text

Amphibians possess both a main and an accessory olfactory system, similar to mammals. However, the amphibian main olfactory epithelium (MOE) still contains a large number of microvillous sensory neurons, in addition to ciliated sensory neurons1. These two sensory neuron populations are to some extent segregated within the MOE234, resulting in corresponding inhomogeneities of the olfactory receptors, since the olfactory receptor gene families expressed in microvillous neurons are different from those found in ciliated neurons25. Neuronal responses to different odorants may be expected to reflect such inhomogeneities. Indeed, responses to amino acids are preferentially localized in the basolateral region of the MOE (microvillous neurons, lateral stream), whereas responses to alcohols, aldehydes and ketones are enriched in the apical and medial region of the MOE (ciliated neurons, medial stream), see6.

Lateral and medial stream segregate more sharply in the olfactory bulb. Amine-responsive glomeruli can be found in both the lateral and the medial olfactory bulb6, suggesting that amine responses may be bimodal, and carried by both odor-processing streams. Due to the sparse occurrence of amine-responding cells their spatial distribution in the amphibian olfactory epithelium has not been investigated yet.

Candidates for amphibian amine receptors would be trace amine-associated receptors (TAARs), since mammalian and fish TAARs have been shown to respond sensitively and specifically to amines789. Taar genes show large evolutionary dynamics resulting in repertoires of over one hundred genes in some ray-finned vertebrate species, e.g. zebrafish, whereas much smaller repertoires are found in the lobe-finned lineage (tetrapods, see10). In amphibians such as Xenopus tropicalis only a minute TAAR family of three genes was reported10, one of which (TAAR1) is not even expressed in the olfactory epithelium of Xenopus6 and other species710. We wished to investigate how such a small gene repertoire becomes integrated into spatial expression patterns dominated by much larger gene families such as odorant receptors (ORs26) and vomeronasal type-2 receptors (V2Rs2). We report here that the two Xenopus taar genes expressed in olfactory tissue, TAAR4a and TAAR4b, show two-dimensional expression patterns clearly distinct from each other and from those of OR and V2R receptors. TAAR4b occupies a basal expression zone, which is homogenous in the lateral-to-medial dimension, in contrast to the lateral-enriched and preferentially apical expression pattern of TAAR4a. Furthermore we show that amine-responsive cells follow a bimodal spatial pattern that could be partially explained by the spatial expression pattern of taar4a and taar4b genes. Our results are consistent with taar genes possessing their own unique expression zones, and possibly mediating amine responses, but also point to the existence of additional amine receptors in amphibian olfaction.

Results

Phylogenetic analysis identifies two additional Xenopus tropicalis taar genes

The taar gene family of Xenopus tropicalis has been predicted to consist of just 3 genes, TAAR1, TAAR4a and TAAR4b10. To search for potential new additions to the databanks we have performed a thorough bioinformatic analysis using the Ensembl genome browser (http://www.ensembl.org/Xenopus_tropicalis/Info/Index), and the NCBI genomic databank (http://www.ncbi.nlm.nih.gov/assembly/GCF_000004195.2). Taar genes from several species (teleost fish as well as mammals) were used as queries. We found two additional members of the TAAR family that show the expected sequence motifs, including the TAAR fingerprint motif11 and that group with the established taar genes in the phylogenetic analysis (Fig. 1). Both genes belong to the TAAR4 subfamily, and were therefore named TAAR4c and TAAR4d. We have, however, not been able to clone orthologs of TAAR4c and TAAR4d from Xenopus laevis, unlike the previously described taar genes TAAR1, TAAR4a, and TAAR4b, for which Xenopus laevis orthologs were obtained by PCR-based homology cloning without difficulties. While we cannot exclude technical difficulties, it is also possible that gene losses have occurred in the Xenopus laevis TAAR4 subfamily compared to Xenopus tropicalis. A gene gain restricted to the Xenopus tropicalis lineage appears less likely, due to the considerable divergence within the TAAR4 subfamily (60–75% identity for pairwise comparisons).

Two closely related taar genes exhibit very different expression zones

In previous experiments246 we had identified two axes, which together describe the spatial distribution of neurons expressing particular olfactory receptor genes. These are laminar height (apical-to-basal dimension) and perpendicular to this axis the position in the medial-to-lateral dimension. We found for several olfactory receptor genes that their distribution is broad, but not random, with preferred positions in both dimensions26, as observed previously for teleosts12.

Here we analysed the spatial patterns of TAAR4a and TAAR4b-expressing olfactory neurons by quantitative in situ hybridization in serial sections of the MOE in larval Xenopus laevis. Both genes show sparse expression as expected for olfactory receptor genes (Fig. 2). For TAAR4a-expressing neurons we find an apical and lateral-enriched distribution (Figs 3 and 4), consistent with previous observations26. However, the closely related TAAR4b shows a very different distribution. No enrichment is seen in the medial-to-lateral dimension, and in the apical-to-basal dimension the distribution of TAAR4b is centered more basally (Figs 3 and 4, Table 1). The differences between TAAR4a and 4b distributions are significant in both dimensions, as analysed by chi square test (p < 10^−7^, medial-to-lateral dimension) and Kolmogorov-Smirnov test (p < 0.001, laminar height) Table 1. Moreover, TAAR4b is expressed somewhat more frequently than TAAR4a, with 1.21 ± 0.08 cells/section compared to 0.27 ± 0.04 cells/section for TAAR4a (mean +/−SEM, n = 139 and 238 respectively; sections devoid of cells are included for calculation of the average). These frequencies are in the range observed for other single genes2.

Amine responses are not lateralized but show a bimodal height distribution

Rodent TAAR receptors appear to be specialized for detection of various amines78, and recently a zebrafish TAAR receptor has also been shown to respond sensitively to some amines9. It is therefore conceivable that detection of amines constitutes an evolutionary conserved feature within the TAAR family. Indeed, amine responses have been shown for Xenopus laevis tadpoles613, but their underlying receptors are not known. If amine responses were mediated by TAARs, we would expect corresponding spatial patterns.

Here we performed a spatial analysis of amine responses in the MOE of Xenopus laevis tadpoles. Odor-induced neuronal activity was measured as intracellular calcium increase in acute slices of the MOE (Fig. 2), using a confocal microscope and fluo-4 as calcium indicator. A mix of thirteen different amines, among them the ligand for mouse TAAR4814, phenylethylamine, was used as odor stimulus. Amine-responding cells were sparse, 3.4 +/− 0.6 cells/section (mean +/−SEM, n = 23), confirming and extending previous observations6. In comparison, amino acid responses were over twofold more frequent, 8.1 +/− 1.3 cells/section (mean +/−SEM, n = 22) Thus, amine responses can be considered to result from a minor neuronal population within the epithelium, i.e. from expression of a few different receptors.

In the medial-to-lateral dimension amine responses showed no pronounced enrichment or depletion, in contrast to amino acid responses which show lateral enrichment (Figs 2, 3, 4). However, in the second dimension analysed, laminar height (basal-to-apical), we observed an unexpected bimodal distribution of amine-responding cells. Two frequency peaks at 0.4 and 0.7 relative height are separated by an intervening region with very few amine responses (Figs 3 and 4, Table 1). The two peaks in the histogram do not depend on the exact size and position of bins used, and are also visible (as plateau) in the empirical cumulative distribution function (ECDF) that does not involve binning (Fig. 3). Such a bimodal distribution suggests a heterogenous composition of the population of amine-responding cells. We would like to emphasize, that no difference in height distribution was observed between the three different compartments in the medial-to-lateral dimension (Fig. 3, Table 1), in contrast to the distribution observed for amino acid responses, where the lateral segment shows a basal peak and the medial segment an apical peak (Fig. 3), see also2. Thus, preferred positions of amine-responding cells appear to be specified independently for each dimension (apical-to-basal and medial-to-lateral).

Partial overlap of TAAR expression zones with amine-responses

Remarkably, the apical peak of the height distribution of amine-responsive cells overlaps largely with the height distribution of TAAR4a-expressing cells, whereas the position of the basal peak is similar to that of TAAR4b-expressing cells (Fig. 3). In both cases the odor response peak is slightly shifted towards more basal positions compared to the peak seen for TAAR-expressing receptor neurons. Using identical assay methods2 we observed previously that the response to an agonist for ciliated neurons was likewise shifted towards slightly more basal positions, compared to the peak seen for a marker of ciliated neurons. In this case it seems clear that the slight difference in peak position ought to be due to technical reasons. Thus we conclude that the apparent small shift between the peak of TAAR-expressing neurons and of amine responses is most likely due to the difference in experimental conditions for both methods.

Both TAAR4b expression and amine-responsive cells show a homogenous distribution in the second dimension, medial-to-lateral. In other words, TAAR4b and the basal partition of amine-responsive cells overlap in both dimensions, consistent with the hypothesis that TAAR4b could mediate (some) amine responses in this region. In contrast, TAAR4a expression is centered in the lateral region, and consequently TAAR4a expression covers only a subregion of the apical partition of amine-responsive cells. Thus TAAR4a may carry lateral-and-apical amine responses, but for medial-and-apical amine responses additional receptors have to be postulated (Fig. 4).

In comparison to two large response domains described previously26, the height peak of the TAAR4a domain and the apical peak of the amine responses are identical to those observed for OMP and forskolin, respectively, both designating the ciliated neuron domain, whereas the height peak for the TAAR4b domain lies in between the ciliated neuron domain and the microvillous neuron domain as defined by TRPC2 expression246.

Discussion

Amines are one of the main odor groups described for amphibians, and have been shown to mediate aversive reactions in mammals1415 and fishes9. In these species olfactory receptors of the TAAR family have been deorphanized as amine receptors79. The amphibian TAAR family has been reported to consist of only three genes (one of which is not expressed in olfactory tissue6), allowing a comprehensive study of expression patterns and a thorough quantitative comparison with the spatial pattern of amine-responsive cells.

We report here that unexpectedly the spatial expression patterns of the remaining two receptors, the closely related TAAR4a and TAAR4b, are very different, both in terms of the height distribution and with respect to medial/lateral enrichment. TAAR4a is expressed apical and lateral, whereas TAAR4b is expressed basal and homogenous in the medial-to-lateral dimension. Thus, two distinctly different expression zones have resulted within the TAAR4 subfamily. Both zones are situated within the neuronal region of the Xenopus MOE, which is sandwiched between a thin layer of supporting cells contacting the lumen, and a basal layer of proliferating cells and immature neurons16.

The TAAR4b pattern is intercalated with respect to height between the spatial expression patterns of MOE-specific v2r genes (basal centered) and of or genes (apical centered), but shares a homogenous distribution in the medal-to-lateral dimension with OMP, a marker for ciliated neurons (note that another marker, G_olf_ is laterally depleted) and TRPC2, a marker for microvillous neurons246, see also (Fig. 4). Since TAAR4b-expressing cells are a minor cell population, they could conceivably be a constituent part of either cell type. However, both mammalian and zebrafish TAARs are known to be expressed in ciliated neurons1718 (and unpublished observation, respectively), suggesting that this feature might be conserved in amphibians.

TAAR4a-expressing cells share their height distribution with ciliated neurons, suggesting an expression in ciliated neurons, but their distinct lateral enrichment in the medial-to-lateral dimension (Fig. 4) is not reflected in the distributions for ciliated marker genes (see above). However, ciliated neurons are present all along the medial-to-lateral dimension (including the lateral segment), and TAAR4a-expressing cells are an even smaller population than TAAR4b-positive cells, so their distribution would not influence the distribution of all ciliated neurons to any measurable degree. Thus, these data are consistent with TAAR4a being expressed in ciliated neurons.

Expression zones for both receptors together cover a large proportion of the amine-response region. Furthermore, amine-responding cells are rare, about 3 cells/explant, consistent with few receptors carrying the amine response. Conceivably, TAAR4a and 4b could mediate a considerable portion of the total amine response in Xenopus laevis. Together these receptors amount to about 1.5 cells/section. However, a direct comparison with the frequency of amine-responding cells is not possible due to the different assay paradigms. In any case, additional amine receptors have to be postulated to explain the apical and medial segment of amine responses, since this region is spared in the TAAR expression patterns analysed here. These receptors might be the putative orthologs of Xenopus tropicalis TAAR4c and 4d, which we were unable to clone, or may be found among both classes of ORs, which show an apical position and varying degrees of lateral depletion6 – such patterns would complement the laterally enriched TAAR4a. Unfortunately, a direct test of this hypothesis is beyond the scope of the present investigation, since the OR family comprises several hundred genes in Xenopus.

Taken together, we have shown that two closely related TAAR receptors exhibit expression zones distinctly different in two dimensions. Together these receptors cover a major proportion of the bimodal amine-responsive region, consistent with at least some of the amine responses being mediated by the Xenopus laevis TAAR4a and TAAR4b.

Materials and Methods

Sequence data mining and phylogenetic analysis

Using Xenopus tropicalis TAAR4a and TAAR4b, as well as representative TAAR sequences from two teleost fishes (zebrafish, medaka) and three mammals (opossum, cow, and mouse) as queries, we searched with tblastn (https://blast.ncbi.nlm.nih.gov/Blast.cgi) for taar genes in the Xenopus tropicalis genome version 4.1 (Joint Genome Institute, JGI, http://genome.jgi-psf.org/Xentr4/Xentr4.info.html). Homology regions above 200 amino acid length were considered further. Candidate sequences had to be located inside the phylogenetic tree with branch support over 80%, had to map to a unique, non-overlapping genomic position, and had to contain the TAAR fingerprint motif11 characteristic for that gene family. Two new sequences fulfilled all criteria, and were named according to phylogenetic relationship as TAAR4c and TAAR4d.

For the phylogenetic analysis published mouse and zebrafish TAAR sequences10 were included as ingroup, and several aminergic receptors as outgroup. The complete set of sequences was aligned using MAFFT, version 6, and E-INS-i strategy with default parameters19. The tree was constructed using a Maximum likelihood algorithm (PhyML-aLRT) with SPR setting for tree optimization and chi square-based aLRT for branch support20, as implemented on the phylemon2 server21. Figtree (http://tree.bio.ed.ac.uk/software/figtree/) was used for visualization of the phylogenetic tree. Branch support above 80% was considered significant.

Animal Handling

All procedures for animal handling were approved by the governmental animal care and use office (Niedersächsisches Landesamt für Verbraucherschutz und Lebensmittelsicherheit, Oldenburg, Germany, Protocol No. T24.07) and were in accordance with the German Animal Welfare Act as well as with the guidelines of the Göttingen University Committee for Ethics in Animal Experimentation. Larval Xenopus laevis, stages 50–54, staged after22, were cooled to produce complete immobility and killed by transection of the brain at its transition to the spinal cord.

Cloning and In Situ Hybridization

Nonambiguous primers for TAAR4b were designed based on published sequence information of homologous sequences in Xenopus tropicalis, taking care to avoid regions of possible cross-reactivity with TAAR4a6. An annealing temperature of 56 °C was used with genomic DNA as template, and a fragment with the predicted length of 516 bp was obtained, cloned into pGEM-T (Promega, Mannheim, Germany) and later confirmed by sequencing.

In situ hybridization was performed using digoxigenin-labeled (DIG; Roche Molecular Biochemicals, Mannheim, Germany) RNA probes prepared from the cloned DNA by using the same forward primers and reverse primers with a T3 promoter site attached to their 5′ end. For in situ hybridization, tissue blocks containing MOE and vomeronasal organ were cut horizontally, fixed in 4% (wt/vol) formaldehyde solution for 2 h at room temperature, equilibrated in 30% saccharose, and embedded in Jung tissue-freezing medium (Leica, Bensheim, Germany). Cryostat sections of 10–12 μm (Leica CM1900) were dried at 55 °C and postfixed in 4% (wt/vol) paraformaldehyde for 10–15 min at room temperature. Hybridizations were performed overnight at 60 °C in 50% (vol/vol) formamide using standard protocols. Anti-DIG primary antibody coupled to alkaline phosphatase and NBT-BCIP (4-nitro blue tetrazolium chloride, and 5-bromo-4-chloro-3-indolyl-phosphate, both from Roche Molecular Biochemicals) were used for signal detection.

Calcium Imaging

Odor responses were measured as changes of intracellular calcium concentrations of individual olfactory neurons in vibratome slices of the olfactory organs using Fluo-4/AM as calcium indicator dye, essentially as described previously223. A mix of thirteen amines (2-Phenylethylamine, tyramine, butylamine, cyclohexylamine, hexylamine, 3-methylbutylamine, N,N-dimethylethylamine, 2-methylbutylamine, 1-formylpiperidine, 2-methylpiperidine, N-ethylcyclohexylamine, 1-ethylpiperidine, piperidine) was used to analyse amine responses, (cf.6). For comparison an L-amino acid mix (19 amino acids, cf.6) was used. Odorant stimuli were prepared at 100 μM concentration per compound and applied by gravity feed into a small carrier flow, i.e. final concentration was somewhat lower than 100 μM. This concentration was chosen in accordance with previous publications using the same assay system236 and has also been used for other aquatic species24. The minimum interstimulus interval was 2 minutes to avoid adaptation. Reproducibility of neuronal responses was verified by repeating the application of each stimulus at least twice. Fluorescence images of the whole MOE were acquired at 1 Hz (excitation at 488 nm; emission above 505 nm) using a laser-scanning confocal microscope (LSM 510/Axiovert 100 M; Zeiss, Jena, Germany) and analyzed using custom programs written in MATLAB (MathWorks, Natick, USA). Active olfactory neurons were identified as regions of high cross-correlation between the fluorescence signals of neighboring pixels, as described in25. The diameter of such regions was typically 6–10 μm, consistent with these signals emanating from the somata of individual olfactory neurons26. Optical section thickness was chosen to ensure that observed signals originated from single cells. Ten images before the onset of stimulus application were taken as control (F_0_). In accordance with previous publications23131623 a response was considered significant if three criteria were fulfilled: it had to occur within 10 s after stimulus application, the first two fluorescence values after stimulus arrival at the mucosa, F_1_ and F_2_, had to be larger than the maximum of the F_0_ values, and F_2_ had to be larger than F_1_.

Analysis of Spatial Distribution

The position of cells was evaluated in two dimensions perpendicular to each other, medial-to-lateral and basal-to-apical. Position in the first dimension was determined according to6, and in the second dimension according to2. The relative height of the cell was defined as distance of the cell soma center from the basal border of the epithelium divided by total thickness of the epithelial layer at the position of the cell (h_rel_ = h_cell_/h_layer_). Cell positions were measured using ImageJ (http://rsbweb.nih.gov/ij/) and/or manually on printouts. Distributions are visualized as histograms with 10 bins (x value given corresponds to the bin center), or unbinned as empirical cumulative distribution function (ECDF). Median, skewness, and half-width of the spatial distributions were calculated from unbinned values using Open Office (version 3.2; www.openoffice.org/). Half-width of a height distribution was defined as difference between the height values for the upper quartile and the lower quartile. The peak value was taken from the graphical representation of the histograms. To estimate whether two height distributions were significantly different, we performed Kolmogorov–Smirnov tests on the unbinned distributions, see27. This test is particularly suitable for continuous distributions and makes no assumptions about the nature of the distributions investigated, which is essential because the observed distributions are not Gaussian. Due to the sensitive nature of the test for large distributions (n > 100), we selected P < 0.01 as cutoff criterion for significant difference.

Additional Information

How to cite this article: Syed, A. S. et al. Different expression domains for two closely related amphibian TAARs generate a bimodal distribution similar to neuronal responses to amine odors. Sci. Rep. 5, 13935; doi: 10.1038/srep13935 (2015).