Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

RNA-Seq Analysis of Human Trigeminal and Dorsal Root Ganglia with a Focus on Chemoreceptors

Abstract

The chemosensory capacity of the somatosensory system relies on the appropriate expression of chemoreceptors, which detect chemical stimuli and transduce sensory information into cellular signals. Knowledge of the complete repertoire of the chemoreceptors expressed in human sensory ganglia is lacking. This study employed the next-generation sequencing technique (RNA-Seq) to conduct the first expression analysis of human trigeminal ganglia (TG) and dorsal root ganglia (DRG). We analyzed the data with a focus on G-protein coupled receptors (GPCRs) and ion channels, which are (potentially) involved in chemosensation by somatosensory neurons in the human TG and DRG. For years, transient receptor potential (TRP) channels have been considered the main group of receptors for chemosensation in the trigeminal system. Interestingly, we could show that sensory ganglia also express a panel of different olfactory receptors (ORs) with putative chemosensory function. To characterize OR expression in more detail, we performed microarray, semi-quantitative RT-PCR experiments, and immunohistochemical staining. Additionally, we analyzed the expression data to identify further known or putative classes of chemoreceptors in the human TG and DRG. Our results give an overview of the major classes of chemoreceptors expressed in the human TG and DRG and provide the basis for a broader understanding of the reception of chemical cues.

Introduction

Detection and processing of environmental chemosensory information is important for living. In addition to gustation and olfaction, the trigeminal system represents a third chemical sense and is provided by the Nervus trigeminus. The nasal cavity is innervated by the ophthalmic and maxillary branches of the trigeminal nerve. The trigeminal ganglia (TG) are located at the base of the scull and represent the cranial analogue of the spinal nerve-associated ganglia, the dorsal root ganglia (DRG). The cell bodies of the pseudounipolar TG and DRG neurons can terminate as free nerve endings in the facial skin and mucosae or the dermatomes of the trunk and extremities, respectively. The DRG efferents terminate in the dorsal horn of the spinal cord, and the TG efferents terminate in specific brainstem nuclei.

Somatosensory neurons in the TG and DRG act as sensors for various innocuous and noxious physical and chemical stimuli. Somatosensory neurons of the TG and DRG can serve as the detectors for a warning system in mammals. Sensations transmitted by the somatosensory system encompass temperatures ranging from freezing cold to painfully hot, stinging and burning elicited by plant compounds and mechanosensations such as touch, stroking, and itch [1,2]. We previously reported that the perception of astringency is also a trigeminal sensation [3].

On the molecular level, several classes of membrane receptors and ion channels are known to be critical for the chemosensory capacity of somatosensory neurons [2,412]. Central players involved in the detection of chemical cues by TG and DRG neurons include transient receptor potential (TRP) channels, two-pore potassium (K2P) channels, and acid-sensing ion channels (ACCN). These channels act as sensors of temperature, mechanical and chemical stimuli, and they are critical for nociception [13,11,5,2]. In addition to ion channels, G protein-coupled receptors (GPCRs) are essential for the detection of a large variety of chemicals [14]. These proteins are the largest superfamily of cell surface proteins and have seven transmembrane domains [15]. They can be activated by either exogenous ligands, such as odorants or tastants, or endogenous ligands, such as neurotransmitters or hormones; thus they play major roles in physiological and pathophysiological processes [16,17]. Known sensory ganglia-associated GPCRs are the mas-related G protein-coupled receptors (MRGPRs), which are, for example, involved in non-histaminic itch [18,19].

It is known that nearly all odorants, particularly at high concentrations, activate intranasal trigeminal nerve fibers and trigeminally innervated respiratory epithelia [20,21]. Odorant activation leads to chemically evoked activity patterns in the TG and higher order neurons in the brain [2224]. Most odorant molecules simultaneously stimulate both, the olfactory and the trigeminal system within the nasal cavity. Anosmic patients can detect and discriminate highly concentrated volatile odorants, supporting the idea that the trigeminal system essentially contributes to odor detection and discrimination [20,25,1,26]. Despite several studies which focused on odorant detection by the trigeminal system, the principles of molecular odor reception by TG neurons are only rudimentarily explored. Several studies showed that a variety of different odorants at high concentrations activate rodent TRP channels [2730]. For example, mammalian TRPV1 is weakly activated by vanillin, helional, heliotropyl acetone, citral, geraniol, thymol, and eugenol [27,3133,30], TRPV3 is activated by monoterpenoids [28], and TRPM8 is activated by linalool, geraniol, hydroxycitronellal and others [27]. Nevertheless, it is likely that additional unidentified trigeminal receptors exist that provide the perception of volatile cues. Our previous sequencing analysis of mouse somatosensory ganglia supports the idea that olfactory receptors (ORs) are expressed in murine sensory ganglia [34]; thus, we focused our chemoreceptor analysis of the TG and DRG on ORs.

To date, no systematic expression analysis of known or potential chemoreceptors expressed in human sensory ganglia is available, and expression analyses of human TG and DRG tissue samples are rare. The expression profiles of TG and DRG for different ion channels or GPCRs have mostly been validated and investigated in rodents [34,35,18,36]. Therefore, we generated and analyzed RNA-Seq data from several human TG and DRG samples. We established the first systematic gene expression profile of the human TG and DRG with a focus on the major classes of chemoreceptors. We subsequently performed microarray analysis, semi-quantitative RT-PCR (semi qRT-PCR), and immunohistochemical staining of human TG and DRG tissue to validate the expression profile of the chosen ORs.

Results and Discussion

General

In our study, RNA-Seq data from human TG and DRG samples were generated using the Illumina sequencing technique. We generated sequencing data for four individual human TG samples (TG 1–4) and a pool of 21 human DRG (Table 1). The data were analyzed with TopHat and Cufflinks software, and reads were mapped to the human reference genome (hg19). These tools are widely accepted and were often used in high resolution transcriptome studies. Cufflinks software normalized reads to account for the length of the gene and the depth of sequencing and ensures comparability of different datasets [37,38]. The quantitative expression values were calculated for each sample based on the number of fragments per kilobase of exon per million fragments mapped (FPKM) [37]. The expression values for all genes in every tissue sample are listed in S1 Table. On a rough scale, >0.1 FPKM corresponds to a very weak expression level, 1 FPKM to a weak expression level, 10 FPKM represents a moderate expression level and 100 FPKM indicates a high expression level. As a basis for comparison, we calculated the FPKM values for typical housekeeping genes. For example, the strongly expressed ß-actin gene yielded an expression value between ~1,000–3,000 FPKM whereas the weakly to moderately expressed TATA box binding protein (TBP) was detected at ~5–9 FPKM (S1 Fig). For an overview of the FPKM values for the expression of all genes, we plotted a histogram for the FPKM value distribution for the TG 3 sample (S2 Fig). Expression analysis revealed ~17,000 genes in the appropriate tissue (FPKM >0.1 out of ~23,000 genes). Gene expression generally reflects the relative abundance of a given gene on protein level, but may not always reflect it. RNA-Seq is an appropriate method to establish an overview for the expression of different receptor and has been used very often [3943,34,38,44]. For further investigations of receptors at the cellular resolution level, in situ hybridization or immunohistochemistry of whole tissue slices might be an additional approach to localize the expression of transcripts or proteins and would indicate whether there are specialized cell types or areas that highly express gene transcripts. We analyzed the transcriptomes of human TG and DRG to assess the expression of chemoreceptor genes in these ganglia. As RNA from the complete ganglia was used for the RNA-Seq experiments, no distinction could be made between different cell types e.g., the neurons and satellite cells that represent the most abundant cell types in sensory ganglia [45]. Thus, even if the total tissue FPKM levels were low, a single transcript’s level may be high in specific cell types. An example of this scenario was shown for the trace amine-associated receptor TAAR1 in brain [46] and for ORs in the olfactory epithelium [42]. The datasets were compared to existing data from a selected panel of human reference tissues (brain, colon liver, lung, skeletal muscle, and testis) [41] and datasets from murine sensory ganglia [34] and provides a basis for further investigations of molecular reception in human sensory ganglia. We confirmed the expression of chemoreceptors which are known to be expressed in sensory ganglia such as TRP or KCNK channels and our study allowed the comparison of different chemoreceptor classes for the first time in human sensory ganglia.

Olfactory Receptors (ORs)

With approximately 400 functional genes and 600 non-functional pseudogenes in humans, ORs form the largest superfamily of GPCRs [47,48]. Most OR genes consist of ~330 amino acids encoded by an intron-free reading frame of approximately 1,000 nucleotides [4951]. Initially, it was postulated that ORs, which detect volatile odorant molecules from the environment, are exclusively expressed in the olfactory epithelium, where they are located in the cilia of olfactory sensory neurons [49]. In 1992, one year after the discovery of ORs, Parmentier et al. found that mammalian OR genes are also expressed in a non-olfactory tissue (testis) by PCR [52]. Since then, a growing number of studies have shown OR expression in several non-olfactory human tissues [41,53,54] and described different physiological functions for these ectopically expressed ORs [5560].

We analyzed the expression of ORs in human TG and DRG samples in a comparative approach using RNA-Seq, microarray analysis, and semi qRT-PCR with a focus on TG 1 sample (Fig 1A). The summarized FPKM values (sFPKM) of the expressed OR genes in the TG (mean sFPKM value of samples TG 1–4) and DRG show that the overall OR expression levels in the TG and DRG are considerably higher (~ 24 sFPKM for TG and ~23 sFPKM for DRG) than in other tissues and comparable to OR expression level in testis (~32 sFPKM) (Fig 1A). Fig 1B shows the OR expression pattern detected by RNA-Seq sorted by the OR expression of the TG 1 sample. Out of the 58 detected ORs in TG 1 (FPKM >0.1), the 30 most highly expressed ORs in comparison to further sensory ganglia and different reference tissues is shown. Based on FPKM values, the RNA for most ORs is present in relatively low abundance compared with that of the housekeeping gene TBP. Out of the 387 annotated OR genes, we detected on average 48 different OR transcripts in a single TG sample and 26 in the DRG sample (FPKM >0.1, S3 Fig). A total of 10 OR transcripts showed FPKM values higher than 1 (OR6B3, OR7A5, OR6B2, OR2L13, OR2W3, OR52B6, OR2A1/42, OR7C1, OR4F21, and OR52H1), and 3 of these transcripts showed values higher than 3 (OR6B3, OR2L13, OR2W3) in one or more tissue samples. The five highest FPKM levels were found for OR2W3 in DRG 1 (6.8 FPKM), OR6B3 in TG 1 (5.2 FPKM), OR2L13 in DRG 1 (3.9 FPKM), OR7A5 in TG 1 (2.6 FPKM), and OR6B2 in TG 1 (2.5 FPKM). The TG and DRG exhibited coherent OR expression patterns and the highest OR transcript levels were consistent for both the TG and DRG. In total, 11 ORs were present in all TG samples, and 39 were detected in at least 2 different TG samples indicating a similar variation and pattern of ORs expressed as already described for human olfactory epithelia [61]. Interestingly, some trigeminally expressed ORs, e.g. OR2W3, were also broadly expressed in several other human tissues [41]. According to our data, some OR transcripts were exclusively present in the TG and DRG (OR6B3). In addition to sensory ganglion tissue, the transcripts of OR6B2 and OR2L13 were weakly detected in testis, colon, or brain. The idea of TG- and DRG-selective receptors is supported by recent RNA-Seq- and PCR-based analyses of the human olfactory epithelium [62,61]. These studies also confirmed that the abovementioned putative TG- and DRG-selective ORs are not or are only weakly expressed in human olfactory epithelium.

thumbnail
Fig 1. OR expression in the human TG and DRG.

A The bar diagram shows the summarized expression (sFPKM) of detected OR transcripts. The sFPKM values of OR gene transcripts in the TG and DRG are in the same range as the summarized expression values found in testis tissue. B The heatmap shows FPKM values for the 30 most abundant OR transcripts found in the TG 1 sample compared with other TG and DRG samples and reference tissues (brain, colon, liver, lung, s.muscle, and testis). Darker colors indicate higher FPKM values, and white indicates the absence of detectable transcripts. The ORs were sorted according to the FPKM values found in the TG 1 sample. As described in [41], the FPKM values of OR2A4 and OR2A7 as well as that of OR2A1 and OR2A42, which share high sequence similarity (99–100%), were analyzed together. C Validation of RNA-Seq results using semi qRT-PCR. The diagram shows the summarized semi qRT-PCR compared with 3 ng of genomic DNA. We investigated 22 different ORs that showed either high, moderate, or low expression in TG 1. D Microarray analysis of the ORs in TG 1. Intensity values for the TG 1 sample were normalized to the reference tissues (liver, heart, s. muscle, skin). Given are all ORs with ratios >1.5. ORs were also detected by RNA-Seq as indicated by red frames. t-test: *p<0.05; **p<0.00013 (after Bonferroni correction).

https://doi.org/10.1371/journal.pone.0128951.g001

We selected 22 ORs that were detectable in the TG 1 RNA-Seq data and confirmed all of the ORs using semi qRT-PCR. The expression levels of ORs determined by PCR correlated with those obtained by RNA-Seq (Fig 1C).

For a long time, microarray was the most commonly used method for transcriptome analysis. Despite the low levels of OR mRNA, microarrays were sufficiently sensitive for detecting ORs in our samples (Fig 1D). Using microarray, we identified 13 ORs with higher expression levels in TG 1 than that in reference tissues (normalized intensities >1.5). The microarray datasets validated the expression of the 5 most highly expressed ORs detected by RNA-Seq (Fig 1D). Furthermore, the microarray data suggested that, in addition to the abovementioned ORs, OR6Y1, OR8S1, OR5A2, OR4X1, OR52E8, and OR51B2 were expressed at higher levels in TG 1 than in reference tissues (Fig 1D). These receptors were not detected by NGS. However, only OR8S1 could be validated by RT-PCR (data not shown), indicating the relatively low reliability of the microarray assay for detecting rare transcripts.

In addition, we analyzed the mapped reads in the RNA-Seq data for the 5 most highly expressed OR genes using the Integrative Genomic Viewer (IGV) [63]. For the highly expressed OR genes OR6B3 and OR6B2, we observed aligned sequences in known exons, and we additionally identified intron-spanning reads that connect newly identified 5′ untranslated gene regions (5′UTRs) with open reading frame (ORF)-containing exons (Fig 2A). For OR2W3, we also detected previously described chimeric transcripts with the upstream Trim58 gene in the DRG but not in the TG samples [41]. For OR2L13, we detected previously annotated 5’UTRs in the TG and DRG via RNA-Seq. All spliced transcripts could be validated by RT-PCR (Fig 2B and S4 Fig).

thumbnail
Fig 2. Analysis of highly expressed OR genes in the human TG and DRG.

A Schematic representation of the newly identified 5’UTRs of OR6B3 and OR6B2. The gene is indicated by blue bars (exon) and thin lines (intron). The coding exon is indicated by ORF (open reading frame) and splice junctions as red arcs. We detected two unannotated 5’UTRs for OR6B3 and OR6B2. B 5’UTR-validation of OR-transcripts using RT-PCR with intron-spanning primers in the DRG sample. The newly identified 5’UTRs of the OR transcripts were confirmed by RT-PCR with a forward primer located in the identified exon and a reverse primer located in the ORF of the respective OR. OR6B3: Ex1 (forward primer in exon 1 of 5’UTR and reverse in OR6B3 ORF); Ex2 (forward primer in exon 2 of 5’UTR and reverse primer in OR6B3 ORF). OR6B2: Ex1 (forward primer in exon 1 of 5’UTR and reverse primer in OR6B2 ORF); Ex2 (forward primer in exon 2 of 5’UTR and reverse primer in OR6B2 ORF). The amplified PCR products were confirmed by Sanger sequencing.

https://doi.org/10.1371/journal.pone.0128951.g002

For the most highly expressed ORs identified in our study, ligands are only known for OR7A5 (myrac aldehyde and 4-hydroxy-2,5-dimethyl-3(2H)-furanone). This receptor was functionally characterized in human testis and spermatozoa [64,65]. OR7C1 showed reliable expression in all TG samples investigated (FPKM values up to 1.7) but not in the DRG sample. This receptor can be activated by androstadienone [66]. Interestingly, this testosterone metabolite is present in human male secretions such as saliva, sweat, and semen and has been implicated as a putative human pheromone [67,68]. However, androstadienone does not induce a specific response in the human vomeronasal duct in comparison with respiratory epithelium [69].

To localize cell types in sensory ganglia expressing OR genes, we performed initial immunohistochemical staining (Fig 3). Due to the limited availability of OR antibodies, we used only those targeted against the most highly expressed ORs (OR6B2, OR6B3, OR7A5). However, specificity tests for recombinantly expressed rho-tagged ORs in HANA3A cells showed that only the OR6B2 antibody specifically detects the corresponding OR protein (S5 and S6 Figs), whereas the other OR antibodies were found to be highly nonspecific and thus not used for further experiments. Due to the high sequence homology of OR6B2 and OR6B3 (94%), we cannot exclude that the specific OR6B2 antibody also detects the OR6B3 protein.

thumbnail
Fig 3. The OR6B2 protein is localized to satellite cells in human DRG.

Immunohistochemical staining of human DRG sections with the OR6B2 antibody indicated OR6B2 protein expression in the neuron-surrounding satellite cells. DAPI staining (blue) was used to determine the number and localization of cell nuclei. The secondary antibody alone did not produce any staining (S4 Fig). Scale bars: 10 μm.

https://doi.org/10.1371/journal.pone.0128951.g003

Subsequent immunohistochemical staining of human DRG slices indicated that OR6B2-positive cells form a thin envelope around unstained neurons (Fig 3). Stained cells seem to be satellite glial cells that tightly envelop neuronal cell bodies to form discrete anatomical units. As sensory ganglia are not protected from substances circulating in the blood, satellite glial cells may be important in the context of exposure to substances including large molecular weight compounds [70]. For example, it is known that satellite glial cells express ATP receptors. Furthermore, TRPV1 protein expression in satellite glial cells in the rat TG was suggested [71]. Sensory neurons or their surrounding cells monitor blood-borne chemicals, and it was proposed that they might play a chemosensory role on their own [72]. It is tempting to speculate that OR6B2 might be involved in this process but further experiments, such as in situ hybridizations, have to be done to give stronger evidences that ORs are expressed in satellite cells. Due to the limited availability of additional specific OR antibodies and TG slices, the question of which cell type, i.e. neurons or satellite cells within the TG and DRG, express other OR proteins remains to be answered.

This study is the first to show OR transcripts in the human TG and DRG. This subset of specific receptors might contribute to the chemosensory capacity of the TG. In addition, ORs might participate in general peripheral neuronal processes in surrounding satellite cells, which can come in contact with chemical substances via the blood vessels in the DRG and possibly in the TG but further experiments have to be done in the future

Further chemosensory GPCRs

In addition to ORs, further GPCRs in the TG and DRG bear the potential for detecting chemical cues. Trace amine-associated receptors (TAARs) are non-canonical ORs that are involved in the detection of social cues [7375]. Recently, it was demonstrated for the first time that human TAAR5 can be activated by amine trimethylamine [76]. TAAR expression is not detectable in the human TG and DRG, and these receptors do not appear to be involved in trigeminal chemosensation (S7 Fig).

The physiological function of human VNO-type chemoreceptors remains elusive. VN1R1 is present in human olfactory epithelium, brain, kidney, liver, and lung [77,78]. Weak expression of VNO-type chemoreceptors can also be detected in sensory ganglia (S8 Fig). Due to their widespread expression throughout the human body, we do not suggest a chemosensory contribution for VNO-type receptors in the human TG and DRG.

Taste receptors are expressed on the tongue where they serve as detectors for different gustatory stimuli. Murine TG and DRG transcriptome analysis exhibited only weak expression of taste receptors in these tissues (FPKM<0.5) [34]. We investigated the expression pattern of all annotated taste receptors in the human TG and DRG via RNA-Seq and detected only a very weak or no expression in human sensory ganglia (S9 Fig). Members of the TAS1R family of taste receptors function as molecular complexes [79]. The heterodimeric TAS1R2/TAS1R3 sweet taste receptor binds sweet stimuli, whereas the TAS1R1/TAS1R3 receptor recognizes amino acids. We detected low transcript levels for the Tas1R3 taste receptor in all sensory ganglia investigated, and for TAS1R1 in 3 of the 4 TG samples and the pool of DRG samples.

Twenty-five different TAS2R genes code for detectors of bitter compounds. Interestingly, it was shown that a variety of bitter taste substances activate rat TG neurons in calcium imaging experiments [80]. In contrast to that, human subjects did not perceive bitterness on the tongue with intact trigeminal innervation after lesion or anesthesia of the taste nerve [3]. However, we found low levels of Tas2R transcripts in sensory ganglia and various other tissues by RNA-Seq. Most detected transcripts of the TAS2R genes lay within the introns of the moderately expressed PRH1-PRR4 gene (S9 Fig). Therefore, the expression of these taste receptor transcripts remains unclear. We suggest that Tas2Rs do not have a specific function in chemosensation in human TG and DRG neurons.

Formyl Peptide Receptor-like proteins (FPRs) are GPCRs found in all mammals, and they are encoded by three genes in humans [81]. FPRs were attributed an olfactory function associated with the identification of pathogens (bacterial peptides) or the pathogenic state [82]. Fig 4 shows that this group of receptors is highly expressed throughout the human body as well as in murine TG and DRG (mTG and mDRG). In the TG and DRG, the expression of this receptor class has not been described. We could show that FPR1 is highly expressed in the human TG and DRG with FPKM values ≤28. For FPR2 transcripts, we found FPKM values of ≤3 and for FPR3, ≤4 in the sensory ganglia investigated. Only in lung tissue did we detect a higher expression of these three receptors and of FPR1 in the brain. It is tempting to speculate that FPRs may be involved in the trigeminally mediated reflexive stop in inspiration that serves to prevent the inhalation of potentially life-threatening pathogens.

thumbnail
Fig 4. Expression of FPRs in the human TG and DRG.

Expression of FPR transcripts could be detected in human (TG and DRG) and murine (mTG and mDRG) sensory ganglia and other reference tissues (brain, colon, liver, lung, s. muscle, and testis).

https://doi.org/10.1371/journal.pone.0128951.g004

MRGPR

The family of mas-related G protein-coupled receptors (MRGPRs) counts 8 members with intact coding sequences in humans [35]. In mice, this receptor class is mainly expressed in subpopulations of sensory neurons [35]. Human MRGPR expression has not been investigated in detail, but exclusive MGPRX expression has been shown in human DRG neurons [18]. Although MRGPRX1 and 2 have been detected in human mast cells [83,84], MRGPRs appear to be specifically expressed in sensory neurons. In general, MRGPRs can be activated by peptides [35,18,85] and chemicals [19]. Furthermore, MRGPR-expressing neurons are involved in detecting mechanical or thermal stimuli [8688].

Our data provide an overview of all MRGPRs expressed in the human TG and DRG in comparison with reference tissues (brain, colon, liver, lung, s.muscle, testis) (Fig 5). RNA-Seq data indicate that the most specific and predominantly expressed MRGPR in human sensory ganglia is MRGPRX1. It is known that this primate-specific MRGPRX1 is enriched in DRG neurons [18] which we can confirm in our datasets of human sensory ganglia (≤9 FPKM in the human TG and DRG) (Fig 5A). MRGPRX1 can be activated by chloroquine and the bovine adrenal medulla peptide BAM8-22, which both induce histamine-independent pruritus [18,19]. Furthermore, MRGPRX1 sensitizes and directly activates TRPV1 via distinct signaling pathways [89]. RNA-Seq data show that MRGPRD is also expressed in the human TG and DRG (≤1.5 FPKM), whereas this receptor is not detectable in human reference tissues (Fig 5A). An involvement of MRGPRD in ß-alanine-mediated pain transmission and its influence on the perception of thermal and mechanical stimuli was suggested [86,87]. Orphan MRGPRE is highly expressed in the human TG and DRG (≤88 FPKM) and nearly absent in reference tissues investigated. We detected high transcript levels for the orphan MRGPRX3 exclusively in the TG and DRG (≤8 FPKM) (Fig 5A). According to our data, MRGPRF is highly expressed in all tissues investigated and therefore the only MRGPR which is also expressed in non-sensory tissues (Fig 5A). In addition, we investigated the distribution of mapped reads in the IGV. We found that human MRGPR transcripts also contain long 3’UTRs as described in mice (S10 Fig).

thumbnail
Fig 5. Expression of MRGPRs in the human TG and DRG.

A MRGPR transcripts were detected in the TG and DRG, whereas they were mostly absent in the reference tissues investigated (brain, colon, liver, lung, s. muscle, and testis). Only MRGPRF transcripts could be detected in all reference tissues. NO = no ortholog B Investigation of the sFPKM values of all MRGPRs with the exception of MRGPRF in all tissues. The expression of MRGPRs is restricted to the sensory ganglia investigated.

https://doi.org/10.1371/journal.pone.0128951.g005

We show for the first time that in humans most MRGPRs are specifically expressed in the TG and DRG and thus are likely to be involved in somatosensation, e.g. itch and pain (Fig 5B).

Ion Channels

Several ion channels, including TRP channels, two-pore potassium channels, and calcium-activated chloride channels are directly and indirectly involved in the perception of chemical and physical stimuli by TG and DRG neurons [5,13,9092]. We examined the expression patterns of different ion channel classes in the human TG and DRG that play a role in chemosensation.

TRP channels

TRP channels comprise a highly conserved superfamily of cation-permeable ion channels [93]. The TRP superfamily is categorized into 7 subfamilies, TRPC, TRPV, TRPM, TRPML, TRPP, TRPA, and TRPN, of which the latter has not been described in mammals and fish [94]. These channels serve a multitude of physiological functions, including the maintenance of intracellular Ca2+ homeostasis and the detection of environmental stimuli. Several TRP channels act as thermosensors, covering a wide range of temperatures, including TRPV1 (> 42°C), TRPV2 (> 52°C), TRPV3 (> 33°C), TRPV4 (25–34°C), TRPA1 (< 17°C), and TRPM8 (< 28°C) [95]. These thermo TRPs include multimodal sensors that also detect either chemical stimuli: TRPV1 (e.g. capsaicin and piperine), TRPV 3 (e.g., thymol and carvacrol), TRPA1 (e.g., cinnamal and allyl-isocyanate), TRPM8 (e.g., menthol and eucalyptol), or mechanical stimulation: TRPV4 [5,96101,2730]. The structures, localization, and physiological functions of TRP channels have been intensely studied during the last fifteen years mostly using rodent model organisms. Despite their multiple important somatosensory functions, limited information for TRP channel gene expression in different human tissues, particularly in human somatosensory ganglia is available. TRP channel expression was shown in several human non-chemosensory tissues such as prostate [102], the gastrointestinal system [103105], eye [106], and brain [107]. The first evidence for functional TRPV1 receptors in human sensory nerves and ganglia was obtained by radioimmunography using [3H]-labeled resiniferatoxin, a TRPV1 ligand [108]. More recent studies show immunoreactivity for TRPA1 and TRPV1 in human lingual nerve neuromas [109,110], for TRPV1 in the human tongue [111] and nerve fibers of the nasal mucosa [112], for TRPA1 in the human DRG [113], for TRPV1, TRPV2, and TRPM8 in human skin nerve fibers [114], and for TRPV1, TRPV3, TRPV4, and TRPM8 in the human DRG and skin [115]. In a functional approach, however, human keratinocytes were resistant to vanilloid-mediated stimulation of TRPV1 and showed a much lower TRPV1 expression level than human sensory ganglia [116] suggesting that capsaicin and other TRPV1 agonists act on peripheral nerves but not on keratinocytes in humans. Data on the TRP gene expression or protein immunolocalization in human trigeminal ganglia is nearly absent. In their study, Hou and coworkers showed that VR1 (TRPV1) expression and VR1 immunoreactivity colocalizes with CGRP in human trigeminal tissue [117]. These available fragmentary data indicate not only the presence but also the functional importance of TRP channels in human sensory tissues.

Our RNA-Seq analysis of the human TG and DRG reveals the tissue selective expression of all TRP channel genes in comparison with non-chemosensory tissues (Fig 6). Here, we found moderate TRPV1 and TRPV2 expression (~10–30 FPKM), which corresponds to data obtained for the mTG (Fig 6A). TRPA1 was also expressed at moderate levels in human TG samples (~3–30 FPKM), which were lower than that in the mTG (Fig 6A). Most strikingly, TRPM8 expression was weak to moderate in our human TG and DRG samples (~0.5–2.5 FPKM), whereas it showed strong expression in the mTG and mDRG (Fig 6A). According to this observation, TRPM8 might play a less prominent role in human than mouse somatosensation. In contrast, TRPV3 expression was moderate in the human TG and DRG but weak in the corresponding mouse ganglia (Fig 6A). In mouse, the warmth sensitive TRPV3 channel is highly expressed in keratinocytes [99]. Upon TRPV3-mediated heat activation, keratinocytes release ATP onto sensory neurons, that in turn, are insensitive to direct activation by heat [118]. According to our expression profiling and immunohistochemistry data [115], TRPV3 might actually be a warmth sensor in human sensory fibers. The thermosensitive TRPV4 channel is a swell-activated osmolarity sensor [97,98,119]. TRPV4 immunoreactivity was shown in human nerve fibers [120,115] and human skin [120,121]. Our RNA-Seq analysis revealed weak TRPV4 expression in all four human TG samples and the DRG with FPKM values comparable to that of TRPV3 (Fig 6A). Similar levels of TRPV4 expression were detected in the mTG and mDRG and human brain, liver, lung, and testis. Furthermore, we found moderate to strong expression of TRPC1, TRPM2, TRPM3, TRPM7, and TRPML1 in the human TG and DRG (Fig 6A). The high expression levels of TRPC1, TRPM7, and TRPML1 were not specific to the TG and DRG but were also found in several non-somatosensory tissues (Fig 6A). Among the non-thermo TRPs, TRPC3 showed weak expression in the human TG and DRG, whereas it was expressed at moderate levels in the corresponding mouse tissues (Fig 6A). We found weak expression of the TRPP channel family member PKD2L1 in human TG and DRG samples and human testis (Fig 6A). This gene is expressed at comparable levels in the mDRG and mTG [34]. In mouse taste buds, PKD2L1 heterodimerizes with PKD1L3 giving rise to acid-sensing channels [122], and its role in somatosensory nerve fibers is unknown.

thumbnail
Fig 6. Expression of TRPs in the human TG and DRG.

A TRP transcripts were detected in the human TG and DRG. B Investigation of the sFPKM values of the TRP transcripts across all tissues. The expression of TRP is higher in human sensory ganglia than in the human reference tissues (brain, colon, liver, lung, s. muscle, and testis) but lower than that in the mTG and mDRG.

https://doi.org/10.1371/journal.pone.0128951.g006

Ca2+-activated chloride channels

ANO.

ANO channels are Ca2+-activated Cl- channels integral to olfactory signaling in rodents [123,124] and somatosensation [125]. We recently showed that Ca2+-activated Cl- currents contribute to signal amplification in murine TG neurons upon the detection of chemical cues [92,3].

To date, no expression data or immunohistochemical analysis of ANO type Ca2+-activated Cl- channels in the human TG and DRG is available. Here, we analyzed the expression of ANO1-10 in human somatosensory ganglia (Fig 7). In general, ANO transcript levels were comparable between the TG and DRG samples. ANO3 through ANO7 are putatively intracellularly localized channels [126]. Due to their intracellular localization, a role in chemosensory signaling is unlikely. FPKM values for these transcripts ranged from <1 FPKM (ANO7), >1 FPKM (ANO3 and ANO4), and >10 FPKM (ANO5 and ANO6) (Fig 7). In RNA-Seq analysis of mTG tissue, Ano3, Ano4, and Ano6 showed the highest FPKM values [92]. ANO1, ANO2, ANO8, and ANO10 give rise to transmembrane Ca2+-activated Cl- currents [127], indicating a possible role in fast signaling events in somatosensory neurons. ANO1 is involved in the bradykinin-mediated depolarization of DRG neurons [128] and acts as a heat sensor in these cells [125]. In the mTG, Ano1 is expressed at low levels [92]. In the human TG and DRG, ANO1 and ANO2 have FPKM values >1, corresponding to weak to moderate expression levels (Fig 7). ANO8 and ANO10 show moderate to high expression in the human TG and DRG. Interestingly, the high FPKM values for ANO8 are specific to the sensory ganglia, whereas the high FPKM values for ANO10 were also detected in the colon, lung, and testis (Fig 7). FPKM values for ANO8 and ANO10 are comparable to that of the mTG and mDRG [92,88], indicating that both Ca2+-activated Cl- channels have the same physiological functions in these sensory ganglia.

thumbnail
Fig 7. Expression of ANOs in the human TG and DRG.

Transcription of ANO channel members was detected in the TG and DRG.

https://doi.org/10.1371/journal.pone.0128951.g007

TTYH.

Tweety channels (TTYH) form another group of Ca2+-activated Cl- channels. We found expression of all three TTYH channels in the human TG and DRG (≤117 FPKM; Fig 8). The expression intensities in the human TG and DRG are relatively high and comparable to that of the mTG and mDRG as well as that of the human brain and testis. To date, the functions of the TTYH channels in human sensory neurons remain unknown.

thumbnail
Fig 8. Expression of TTYHs in the human TG and DRG.

The transcription of all TTYH channel members was detected in the TG and DRG. The expression of the TTYH genes is higher in the human sensory ganglia than the human reference tissues with the exception of the brain and testis.

https://doi.org/10.1371/journal.pone.0128951.g008

ACCN.

Amiloride-sensitive cation channels (ACCN), also known as acid-sensing ion channels (ASICs), are voltage-insensitive, proton-gated Na+ channels belonging to a superfamily including the epithelial Na+ channel. ACCN subunits possess two putative TM domains and are assembled as hetero- or homotrimers [129,130]. ACCN proteins 1–3 were previously detected in a sub-population of small-diameter sensory neurons in the human DRG [131]. Our datasets confirm that these three receptors are highly expressed in the human TG and DRG with FPKM values ≤27 (Fig 9). ACCN3 is known to be involved in multimodal sensory perception in the sensory ganglia, including nociception, mechanosensation, and chemosensation [132]. ACCN4 and 5 are not or are only weakly expressed in the human TG and DRG (Fig 9).

thumbnail
Fig 9. Expression of ACCN in the human TG and DRG.

ACCNs 1–3 are highly expressed in sensory ganglia.

https://doi.org/10.1371/journal.pone.0128951.g009

Piezo receptors

FAM38A and B, also known as Piezo1 and 2, are mechanically activated cation channels and may be involved in sensations such as touch and pain [133]. Their expression has not been previously described in the human TG and DRG. FAM38A and B are highly expressed in the TG and DRG with FPKM values ≤17 (Fig 10).

thumbnail
Fig 10. Expression of piezo receptors in the human TG and DRG.

FAM38A and B are highly expressed in both sensory ganglia.

https://doi.org/10.1371/journal.pone.0128951.g010

Purinergic receptors

P2X receptors are involved in a wide range of pathophysiological pain mechanisms, such as inflammatory, neuropathic, acute, and migraine-induced pain [134]. These receptors can define chemosensory properties in rat trigeminal neurons and may contribute to odorant perception via the trigeminal system [10]. Furthermore, we previously showed that keratinocytes act as chemosensors linking the environment and the trigeminal system via ATP signaling [135].

We identified the expression of the ionotropic purinergic receptor class P2X in the human TG and DRG (Fig 11). The expression of P2RX2-7 has been well described in the sensory ganglia of rodents but not in humans. Our RNA-Seq results revealed that P2RX3 is highly expressed in sensory ganglia, whereas it is almost absent in reference tissues investigated (Fig 11). The expression of P2RX3 has been described for human fetal DRG, where it is generally localized in small-sized neurons (<20 μm diameter) [136]. Moreover, the high P2RX3 expression is consistent with its expression in mTG and mDRG and for all of the additionally expressed P2RXs [34] (Fig 11). P2RX1 and P2RX4-7 are broadly and moderately expressed in almost all human tissues investigated, but they are highly expressed in sensory ganglia (Fig 11).

thumbnail
Fig 11. Expression of P2XRs in the human TG and DRG.

The ionotropic purinergic receptors are highly expressed in sensory ganglia. Five of the seven P2XRs are expressed higher in the TG and DRG than in any other tissue.

https://doi.org/10.1371/journal.pone.0128951.g011

Potassium channels

In addition to TRP channels, potassium channels are involved in somatosensation and nociception in sensory nerves [137,138]. Potassium channels are subclassified as voltage-gated (KCNA-KCND, KCNF-KCNH, KCNQ and KCNS), calcium-activated (KCNM-KCNN), inward rectifying (KCNJ), and background/leak, 2 pore (KCNK) channels [139]. Expression analysis of all potassium channels can be found in S11 Fig.

KCNK

The 15 KCNK channel members are a class of thermo-, mechano-, and chemosensitive membrane proteins that are involved in nociception and somatosensation [137]. KCNKs 3, 9, and 18 function as chemoreceptors for hydroxyl α-sanshool in murine TG neurons [90]. Additionally, these three human KCNKs are inhibited by the pungent substances piperine, capsaicin, 6-gingerol, and polygodial [91]. Furthermore, human KCNKs 2, 10, and 18 are blocked by high concentrations of (-)-nicotine [140]. Expression in human DRG was already described for KCNK1-7, 10 and, 18 [141,142].

Our RNA-Seq data allowed for the complete expression analysis of human DRG and TG and revealed the expression of virtually all human annotated KCNK channels in all sensory ganglia investigated (Fig 12). As an exception, KCNK16 was detected only in one TG sample (TG 2). We detected the highest mean transcript levels (mean FPKM) in the TG and DRG for KCNK12 (mean FPKM: 62.3), KCNK3 (mean FPKM: 7.1), and KCNK1 (mean FPKM: 6.5). For KCNK12, which was most highly and specifically expressed, the function is still unknown [143]. The three most highly expressed KCNK genes in murine sensory ganglia differed in the human expression results (Kcnk1, Kcnk18, and Kcnk2). Our data indicate the exclusive expression for KCNK18 in the human TG and DRG (mean FPKM: 2.8) (Fig 12), which was also detected with RT-PCR experiments by Lafrenière and colleagues [142]. Murine RNA-Seq data also indicate selective KCNK18 expression in the mTG and mDRG [34]. A mutation within the KCNK18 open reading frame leads to a non-functional protein and has been found to be associated with migraines [142].

thumbnail
Fig 12. Expression of KCNK channels in the human TG and DRG.

All KCNK genes are expressed in the TG and DRG. The KCNK genes are sorted by the mean of their expression values across all human sensory ganglia.

https://doi.org/10.1371/journal.pone.0128951.g012

Other GPCRs and ion channels

Furthermore, we analyzed the transcript levels of all annotated non-olfactory GPCRs. A list of non-olfactory GPCR genes (375 different genes) was generated based on the IUPHAR database. In total, we detected 273 and 260 GPCRs expressed in the human TG (mean number of TG 1–4) and DRG (FPKM > 0.1). We generated a ranking of the 30 most selectively expressed GPCRs in the human TG and DRG (Fig 13). Among the 30 most selectively expressed GPCR transcripts in the human TG and DRG, we detected GPCRs known to be involved in chemoreception (e.g., MRGPRX1 [19]), nociception (e.g., MRGPRD [86]), migraines (e.g., HTR1D [144]), and inflammation (e.g., F2RL2 [145]). The three most selectively expressed GPCRs in the human TG are MRGPRX1 (mean FPKM: 1.8), MRGPRE (mean FPKM: 17) and MRGPRX3 (mean FPKM: 3.4). In the DRG, we detected GPR139 (FPKM: 1.1), MRGPRX1 (FPKM: 8.7) and MRGPRE (FPKM: 66.9) as the three most selectively expressed GPCRs. The human TG and DRG share 18 of the 30 most highly expressed GPCRs.

thumbnail
Fig 13. Expression analysis of the 30 most selectively expressed GPCRs in the human TG and DRG.

A The heatmap shows the 30 most selectively expressed GPCR transcripts in the TG in comparison with the DRG and reference tissues (brain, colon, liver, lung, s. muscle, and testis). Genes are ranked according to their selective expression in the TG (Quotient of the mean FPKM value (TG1-4) and the mean FPKM value for all reference tissues; only the transcripts detected in all four TG samples with FPKM >0.1 are taken into account). MRGPRs are the most selective GPCRs that have been detected. B Shown are the 30 most selectively expressed GPCRs in the DRG. Genes are ranked according to their selective expression in DRG. The ranking was calculated by the quotient of the FPKM value of DRG (FPKM >0.1) and the mean FPKM value for all reference tissues. GPR139 is the most specific GPCR followed by three MRGPRs, which have been selectively detected in the human DRG.

https://doi.org/10.1371/journal.pone.0128951.g013

Next, we investigated the expression of the most common ion channels in the human TG and DRG (351 different genes) based on [146]. In total, we found 304 ion channels in the TG and 299 in the DRG to be expressed (0.1> FPKM). Fig 14 shows the 30 most selectively expressed ion channels in the human TG and DRG. The expression pattern of the selectively expressed ion channel genes is similar between the TG and DRG, and 21 ion channels are the most selectively expressed ion channels in both tissues. SCN10A, KCNK18, and SCN11A are the three most selective ion channel genes in both, the TG and DRG. More than half of the 30 most selectively expressed ion channels in the human TG and DRG play a role in diverse sensory functions. Among these channels are the previously described KCNK, TRP, P2RX, and ACCN channels and the FAM38 receptors. SCN10A and 11A, which belong to the family of voltage-gated sodium channel α-subunits SCNA were found to be highly selectively expressed in the human TG and DRG. These channels are known to be involved in orofacial or trigeminal neuropathic pain and toothache [147150]. In the TG and DRG, we detected the same four members of the nicotinic acetylcholine receptors (CHRNA6, 7, 9, and B3), which possibly act as chemosensors for (-)-nicotine [140]. Ionotropic kainate glutamate receptor 3 (GRIK3) is selectively expressed in the human TG and DRG. In insects, this receptor is associated with chemoreception [151]. In humans, polymorphisms within this gene are associated with several diseases such as schizophrenia [152].

thumbnail
Fig 14. Expression analysis of the 30 most selectively expressed ion channels in the human TG and DRG.

A The heatmap shows the 30 most selectively expressed ion channels in the TG compared with the DRG and reference tissues (brain, colon, liver, lung, s.muscle, and testis). Genes are ranked according to their selective expression in the TG (Quotient of the mean FPKM value (TG1-4) and the mean FPKM value of all reference tissues, only those genes expressed in all four TG with FPKM >0.1 were taken into account), and of which SCN10A is the most selectively expressed ion channel in the human TG. B Shown are the 30 most selectively expressed ion channels in the DRG. The genes are ranked according to their selective expression in the DRG. The ranking was calculated by the quotient of the FPKM value in the DRG (FPKM >0.1) and the mean FPKM value for all reference tissues. SCN10A is the most selective ion channel detected in the DRG.

https://doi.org/10.1371/journal.pone.0128951.g014

Genes associated with olfactory signaling

Next, we determined the expression pattern of the signaling components of a potential olfactory signal transduction cascade. In olfactory neurons, ORs couple to a cAMP-mediated second messenger cascade that involves Gαolf, adenylyl cyclase III, and the subunits of the CNG channel CNGA2, CNGA4 and CNGB1 [153]. Two of the essential components of olfactory signal transduction, the Gα subunit Gαolf and adenylyl cyclase III (ADCY3), are highly expressed in the human sensory ganglia, but no pronounced expression of the corresponding olfactory CNG channel was detected (Fig 15). Lately, the olfactory CNG channel CNGA2 was identified in mTG tissue by RNA-Seq [34] indicating the presence of a cAMP target channel in mTG neurons. In human sensory ganglia, ORs could couple to Golf and increase cAMP by ADCY3 activation; however, cAMP does not trigger Ca2+ influx similar to olfactory sensory neurons by opening olfactory CNG channels. In other tissues, ectopically expressed ORs use different signaling mechanisms. In keratinocytes OR activation leads to the opening of the rod CNG channel [59], and TRPV6 is involved in the OR signal transduction cascade in prostate cancer cells [57].

thumbnail
Fig 15. Expression of olfactory signal transduction components in the human TG and DRG.

Expression pattern of signal transduction components including adenylyl cyclase III (ADCY3), the calcium-activated chloride channel (ANO2), the classical olfactory and rod CNG channel subunits (CNGA1, CNGA2, CNG3, CNGA4, CNGB1 and CNGB3), and Gαolf (GNAL). We also investigated the expression of accessory proteins including receptor-transporting proteins (RTP1 and RTP2) and receptor-enhancing proteins 1 (REEP1) as well as the expression of olfactory marker protein (OMP), a specific marker for olfactory sensory neurons.

https://doi.org/10.1371/journal.pone.0128951.g015

Conclusions

Somatosensation by the TG and DRG is an important characteristic of human sensory physiology that is relevant to nutrition, the enjoyment of food, peripheral mechanisms, such as itch and pain, the regulation of temperature, physiology of the skin, environmental toxicology, eye pathologies, and respiratory diseases. Our results provide a framework for insight into the molecular basis of these different sensations and processes. We detected and described a broad panel of different chemoreceptors, and the receptor expression patterns align with the known functions of the TG and DRG. The expression of some receptors, such as ORs, was described here for the first time. In addition to testis tissue, the TG could be defined as a tissue with the highest number of expressed OR genes.

Materials and Methods

Ethics Statement

Human trigeminal ganglia samples were obtained from the archives of the institution of pathology of the Ruhr-University Bochum following appropriate institutional approval by the ethical committee of the medical faculty of the Ruhr-University Bochum (file number 5145–14). The ethical committee waived the need for written consent. The samples were received and analyzed anonymously. The committee did not raise any objections against this study.

Human tissues and RNA isolation

A total RNA pool from 21 different human DRG was purchased from a commercial source (Clontech, CA, USA). Human trigeminal ganglia samples were obtained at autopsy from male adult subjects at the institution of pathology of the Ruhr-University Bochum. The isolation of RNA was performed using the RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The RNA samples were subjected to DNaseI treatment using the TURBO DNA-free Kit (Life Technologies, Carlsbad, CA, USA) according to the standard protocol.

RNA-Seq by next generation sequencing

At the Cologne Center for Genomics NGS unit, libraries for sequencing were constructed from total or mRNA. Libraries from total RNA were subjected to DSN normalization. RNA-Seq was performed on the Illumina GAIIx (read length: 36 bp single and paired-end) or HiSeq 2000 sequencing platforms (101 bp paired-end). The Body Map 2.0 data used in this study as a reference were obtained as previously described [41].

Alignment of RNA-Seq reads using TopHat

We analyzed the RNA-Seq data as previously described [41]. RNA-Seq reads were aligned to the hg19 reference genome by TopHat v1.2.0. Aligned data were visualized with the Integrative Genomic Viewer [63]. The command line used for TopHat was as follows:

tophat—output-dir <name output>—GTF <hg19refseq.gtf> <indexes> tissue.fq

Alignment assembly and gene expression using Cufflinks

Cufflinks v1.0.3 software was used to calculate the abundance of transcripts based on the RefSeq gene model as described in [41]. The relative abundance of transcripts was reported in FPKM (fragments per kilobase of exon per million fragments mapped) units [154]. The Cufflinks parameters are listed below.

cufflinks—output-dir <name output2>—GTF <hg19refseq.gtf>—multi-read-correct—compatible-hits-norm—minfrags-per-transfrag 1—frags-bias-correct <hg19.fa> sorted.bam

cDNA synthesis and semi quantitative RT-PCR

cDNA synthesis was performed according to the previously described protocol of [41]. An equivalent of 50 ng of total RNA was used for each RT-PCR experiment. To validate the expression of different ORs, we designed primers that detect ~100–300 bp fragments of OR ORFs. PCR was performed using GoTaq qPCR Master Mix (Promega, Madison, WI, USA) with the Mastercycler RealPlex2 (Eppendorf, Hamburg, Germany) (20 μl total volume, 40 cycles of 95°C, 59°C, and 72°C at 45 s each). All experiments were conducted in triplicate. For a semi-quantitative analysis we normalized probes to 3 ng of genomic DNA. A list of primers used can be found in S12 Fig.

Microarray analysis

Gene expression analysis of trigeminal ganglia (n = 3) and reference tissues (liver, heart, muscle, skin) was performed using a customized Agilent Whole Human Genome Microarray 8x15k from Miltenyi Biotech GmbH (Bergisch Gladbach, Germany). The microarray contains triplicate probes 50 nt in length for ~850 OR genes or pseudogenes. Total RNA was isolated as described above. The quality of the total RNA was controlled for using the Agilent 2100 Bioanalyzer System. Linear amplification of the RNA and hybridization reactions as performed according to Agilent’s standard protocols. Pre-analysis of the data, including median normalization of the intensity values, was performed by the Miltenyi Biotec GmbH using Agilent Feature Extraction Software. For the expression analysis of ORs in TG 1, normalized mean intensities of three oligonucleotides per array were normalized to different reference tissues.

Cell culture and transfection

Hana3A cells were used for recombinant OR expression to test specificity of the OR antibodies. The Hana3A cells were kindly provided by H. Matsunami (Duke University Medical Center, Durham, NC). The Hana3A cells were maintained under standard conditions as previously described [76] and grown on cover slips in 24-well plates. Cells were transfected with Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol with 300 ng of OR plasmid and 60 ng of mRTP1S plasmid (constructed using standard PCR methods).

For the antibody specificity studies, the OR6B2 coding sequence was PCR amplified from human genomic DNA using specific primer pairs (for: 5’-GCATATGAATTCATGAGTGGGGAGAATGTCACC-3’ and rev: 5’-GCATATGCGGCCGCTCAGTGTGAAGTTTGACCCAAGC-3’) that amplify the complete open reading frame and contain the EcoRI and NotI restriction sites for further subcloning into the pCI plasmid (Promega, Madison, USA), which contain the coding sequence for the N-terminal rhodopsin tag (rho-tag, first 20 amino acids of rhodopsin). All plasmid constructs and PCR products were verified by Sanger sequencing.

Immunocytochemical and—histochemical staining

Transfected Hana3A cells grown on cover slips were fixed with 4% paraformaldehyde. Commercially available human DRG paraffin tissue sections with a thickness of 5 μm (Zyagen, San Diego, CA, USA) were deparaffinized by dipping slides in xylene (2 changes, 5 min each), rehydration in a serial ethanol solutions, and finally washing in PBS-/-. For the deparaffinized slides, antigen retrieval with citrate buffer was performed (20 min, 90°C). Afterward, fixed cells and deparaffinized tissue slides were permeabilized with PBS-/- + 0.1% Triton X-100, washed with PBS-/-, and incubated with blocking reagent (PBS-/- + 0.1% Triton X-100, 5% NGS, and 1% fish gelatine) for 1 h. Cells and tissue slides were incubated overnight (4°C) with primary antibody and then incubated with the indicated Alexa-conjugated secondary antibodies (Invitrogen, Carlsbad, USA) and DAPI for 45 min at room temperature in the dark. The following antibodies were used: anti-OR6B2 (polyclonal, 1:30, C-terminal, NBP1-71360, Novus Biologicals, Cambridge, UK) and anti-Rhodopsin 4D2 (monoclonal, 1:250, Abcam, Cambridge, UK). The samples were mounted in ProLong Gold (Invitrogen, Carlsbad, USA). Images were obtained using a confocal fluorescent microscope (LSM 500 Meta; Carl Zeiss, Oberkochen, Germany).

Supporting Information

S1 Fig. Expression of housekeeping genes in the human TG, DRG, and reference tissues.

Shown are highly expressed genes (ß-actin (ACTB) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH)), moderately expressed genes (ribosomal protein L29 (RPL29) and ribosomal protein L13A (RPL13A)) and weakly expressed genes (β-glucuronidase (GUSB), transferrin receptor (TFRC), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and TATA box binding protein (TBP)).

https://doi.org/10.1371/journal.pone.0128951.s001

(TIF)

S2 Fig. Distribution of the FPKM values in the human TG.

To gain an estimate of the FPKM values for expressed genes, a histogram of FPKM distribution for the TG 3 sample was calculated. Values <0.3 can be regarded as indicating weak expression, 0.3–3 as weakly expressed and 3–30 as moderately expressed. Values of 30–100 indicate high expression, and values >100 indicate extremely high expression. Of the ~23,000 analyzed genes, expression at >0.1 FPKM was detected for ~17,000 genes; mRNA for ~500 of these genes were extremely highly expressed with FPKM values >100.

https://doi.org/10.1371/journal.pone.0128951.s002

(TIF)

S3 Fig. Number of expressed OR transcripts in the human TG and DRG.

Each bar represents the number of OR genes that were expressed in each tissue with an FPKM value >0.1. For the TG, the average number of ORs is shown.

https://doi.org/10.1371/journal.pone.0128951.s003

(TIF)

S4 Fig. Validation of detected OR splicing events for OR2W3 and OR2L13 by RT-PCR.

Detected splicing events could be validated by RT-PCR; and example is shown for DRG. OR2W3: Ex3 (forward primer in exon 3 of TRIM58 and reverse primer in OR2W3 ORF); Ex5 (forward primer in exon 5 of TRIM58 and reverse primer in OR2W3 ORF). OR2L13: Ex1 (forward primer in known exon 1 of 5’UTR and reverse primer in OR2L13 ORF); Ex2 (forward primer in known exon 2 of 5’UTR and reverse primer in OR2L13 ORF). The amplified PCR products were confirmed by Sanger sequencing.

https://doi.org/10.1371/journal.pone.0128951.s004

(TIF)

S5 Fig. The OR6B2 antibody specifically detects recombinantly expressed OR6B2 in Hana3A cells.

Immunostaining of Hana3A cells transiently transfected with OR6B2. The cells were stained with a specific OR6B2 antibody (OR6B2, green) and a rhodopsin-antibody (rho, red). DAPI staining (blue) was used to confirm the number and localization of cell nuclei. Scale bars: 20 μm.

https://doi.org/10.1371/journal.pone.0128951.s005

(TIF)

S6 Fig. Secondary antibody control in human DRG slides.

Control staining was performed without primary antibody, showing non specific staining by the secondary antibody (Alexa Fluor 488 Goat Anti-Rabbit, Control). Scale bar: 20 μm.

https://doi.org/10.1371/journal.pone.0128951.s006

(TIF)

S7 Fig. Expression of TAARs in the human TG and DRG.

No TAAR transcripts could be detected in the sensory ganglia investigated.

https://doi.org/10.1371/journal.pone.0128951.s007

(TIF)

S8 Fig. Expression of VNRs in the human TG and DRG.

Shown are the FPKM values for VNR transcripts in sensory ganglia compared to different reference tissues (brain, colon, liver, lung, s. muscle, and testis).

https://doi.org/10.1371/journal.pone.0128951.s008

(TIF)

S9 Fig. Expression of taste receptors in the human TG and DRG.

The FPKM values for TAS1Rs and TAS2Rs are shown. The sequencing read distribution was checked in the IGV, and in some cases, unclear expression in the human TG and DRG was revealed. Most of the detected TAS2R lay within introns in the moderately expressed PRH1-PRR4 gene. Intronic reads within the PRH1-PRR4 gene could stem from unprocessed transcripts of PRH1-PRR4 and not TAS2R transcripts as described in [155].

https://doi.org/10.1371/journal.pone.0128951.s009

(TIF)

S10 Fig. The MRGPR transcript in the human TG.

A large 3’UTR could be detected in transcripts from human MRGPRs (the example shown is MRGPRE), which is found in mice.

https://doi.org/10.1371/journal.pone.0128951.s010

(TIF)

S11 Fig. Expression of potassium channels in the human TG and DRG.

Genes are sorted by the mean of their expression values across all human sensory ganglia.

https://doi.org/10.1371/journal.pone.0128951.s011

(TIF)

S12 Fig. Primer sequences used for PCR and splice transcript validation.

https://doi.org/10.1371/journal.pone.0128951.s012

(TIF)

S1 Table. Expression of all genes in the human TG and DRG.

Shown are the FPKM values for all expressed genes in sensory ganglia compared to different reference tissues (brain, colon, liver, lung, s. muscle, and testis).

https://doi.org/10.1371/journal.pone.0128951.s013

(XLSX)

Acknowledgments

We thank T. Lichtleitner, S. Pyschny, A. Stöck, and A. Sommer for excellent technical support.

Author Contributions

Conceived and designed the experiments: HH GG CF. Performed the experiments: CF JA CB. Analyzed the data: CF GG. Contributed reagents/materials/analysis tools: HH AT. Wrote the paper: CF NS.

References

  1. 1. Laska M, Distel H, Hudson R. Trigeminal perception of odorant quality in congenitally anosmic subjects. Chem. Senses 1997; 22(4):447–56. pmid:9279467
  2. 2. Viana F. Chemosensory properties of the trigeminal system. ACS Chem Neurosci 2011; 2(1):38–50. pmid:22778855
  3. 3. Schöbel N, Radtke D, Kyereme J, Wollmann N, Cichy A, Obst K, et al. Astringency Is a Trigeminal Sensation That Involves the Activation of G Protein-Coupled Signaling by Phenolic Compounds. Chem. Senses 2014.
  4. 4. Silver WL, Farley LG, Finger TE. The effects of neonatal capsaicin administration on trigeminal nerve chemoreceptors in the rat nasal cavity. Brain Res. 1991; 561(2):212–6. pmid:1724948
  5. 5. Caterina MJ, Schumacher MA, Tominaga M, Rosen TA, Levine JD, Julius D. The capsaicin receptor: a heat-activated ion channel in the pain pathway. Nature 1997; 389(6653):816–24. pmid:9349813
  6. 6. Bonaventure P, Voorn P, Luyten WH, Leysen JE. 5HT1B and 5HT1D receptor mRNA differential co-localization with peptide mRNA in the guinea pig trigeminal ganglion. Neuroreport 1998; 9(4):641–5. pmid:9559931
  7. 7. Liu L, Chang GQ, Jiao YQ, Simon SA. Neuronal nicotinic acetylcholine receptors in rat trigeminal ganglia. Brain Res. 1998; 809(2):238–45. pmid:9853116
  8. 8. Zhu JJ, Uhlrich DJ. Cellular mechanisms underlying two muscarinic receptor-mediated depolarizing responses in relay cells of the rat lateral geniculate nucleus. Neuroscience 1998; 87(4):767–81. pmid:9759965
  9. 9. Alimohammadi H, Silver WL. Evidence for nicotinic acetylcholine receptors on nasal trigeminal nerve endings of the rat. Chem. Senses 2000; 25(1):61–6. pmid:10667995
  10. 10. Spehr J, Spehr M, Hatt H, Wetzel CH. Subunit-specific P2X-receptor expression defines chemosensory properties of trigeminal neurons. Eur. J. Neurosci. 2004; 19(9):2497–510. pmid:15128403
  11. 11. McKemy DD, Neuhausser WM, Julius D. Identification of a cold receptor reveals a general role for TRP channels in thermosensation. Nature 2002; 416(6876):52–8. pmid:11882888
  12. 12. Gu Y, Huang LY. Modulation of glycine affinity for NMDA receptors by extracellular Ca2+ in trigeminal neurons. J. Neurosci. 1994; 14(7):4561–70. pmid:8027794
  13. 13. Bautista DM, Siemens J, Glazer JM, Tsuruda PR, Basbaum AI, Stucky CL, et al. The menthol receptor TRPM8 is the principal detector of environmental cold. Nature 2007; 448(7150):204–8. pmid:17538622
  14. 14. Chen G, Way J, Armour S, Watson C, Queen K, Jayawickreme CK, et al. Use of constitutive G protein-coupled receptor activity for drug discovery. Mol. Pharmacol. 2000; 57(1):125–34. pmid:10617687
  15. 15. Zhang Y, Devries ME, Skolnick J. Structure modeling of all identified G protein-coupled receptors in the human genome. PLoS Comput. Biol. 2006; 2(2):e13. pmid:16485037
  16. 16. Pierce KL, Premont RT, Lefkowitz RJ. Seven-transmembrane receptors. Nat. Rev. Mol. Cell Biol. 2002; 3(9):639–50. pmid:12209124
  17. 17. George SR, O'Dowd BF, Lee SP. G-protein-coupled receptor oligomerization and its potential for drug discovery. Nat Rev Drug Discov 2002; 1(10):808–20. pmid:12360258
  18. 18. Lembo PMC, Grazzini E, Groblewski T, O'Donnell D, Roy M, Zhang J, et al. Proenkephalin A gene products activate a new family of sensory neuron—specific GPCRs. Nat. Neurosci. 2002; 5(3):201–9. pmid:11850634
  19. 19. Liu Q, Tang Z, Surdenikova L, Kim S, Patel KN, Kim A, et al. Sensory neuron-specific GPCR Mrgprs are itch receptors mediating chloroquine-induced pruritus. Cell 2009; 139(7):1353–65. pmid:20004959
  20. 20. Doty RL, Brugger WE, Jurs PC, Orndorff MA, Snyder PJ, Lowry LD. Intranasal trigeminal stimulation from odorous volatiles: psychometric responses from anosmic and normal humans. Physiol. Behav. 1978; 20(2):175–85. pmid:662939
  21. 21. Kobal G, Hummel C. Cerebral chemosensory evoked potentials elicited by chemical stimulation of the human olfactory and respiratory nasal mucosa. Electroencephalogr Clin Neurophysiol 1988; 71(4):241–50. pmid:2454788
  22. 22. Hummel T, Doty RL, Yousem DM. Functional MRI of intranasal chemosensory trigeminal activation. Chem. Senses 2005; 30 Suppl 1:i205–6. pmid:15738116
  23. 23. Rothermel M, Ng BSW, Grabska-Barwińska A, Hatt H, Jancke D. Nasal chemosensory-stimulation evoked activity patterns in the rat trigeminal ganglion visualized by in vivo voltage-sensitive dye imaging. PLoS ONE 2011; 6(10):e26158. pmid:22039441
  24. 24. Lübbert M, Kyereme J, Rothermel M, Wetzel CH, Hoffmann K, Hatt H. In vivo monitoring of chemically evoked activity patterns in the rat trigeminal ganglion. Front Syst Neurosci 2013; 7:64. pmid:24115922
  25. 25. Mason JR, Silver WL. Trigeminally mediated odor aversions in starlings. Brain Res. 1983; 269(1):196–9. pmid:6871699
  26. 26. Thuerauf N, Kaegler M, Dietz R, Barocka A, Kobal G. Dose-dependent stereoselective activation of the trigeminal sensory system by nicotine in man. Psychopharmacology (Berl.) 1999; 142(3):236–43. pmid:10208315
  27. 27. Behrendt H, Germann T, Gillen C, Hatt H, Jostock R. Characterization of the mouse cold-menthol receptor TRPM8 and vanilloid receptor type-1 VR1 using a fluorometric imaging plate reader (FLIPR) assay. Br. J. Pharmacol. 2004; 141(4):737–45. pmid:14757700
  28. 28. Vogt-Eisele AK, Weber K, Sherkheli MA, Vielhaber G, Panten J, Gisselmann G, et al. Monoterpenoid agonists of TRPV3. Br. J. Pharmacol. 2007; 151(4):530–40. pmid:17420775
  29. 29. Sherkheli MA, Vogt-Eisele AK, Bura D, Beltrán Márques LR, Gisselmann G, Hatt H. Characterization of selective TRPM8 ligands and their structure activity response (S.A.R) relationship. J Pharm Pharm Sci 2010; 13(2):242–53. pmid:20816009
  30. 30. Lübbert M, Kyereme J, Schöbel N, Beltrán L, Wetzel CH, Hatt H. Transient receptor potential channels encode volatile chemicals sensed by rat trigeminal ganglion neurons. PLoS ONE 2013; 8(10):e77998. pmid:24205061
  31. 31. Xu H, Delling M, Jun JC, Clapham DE. Oregano, thyme and clove-derived flavors and skin sensitizers activate specific TRP channels. Nat. Neurosci. 2006; 9(5):628–35. pmid:16617338
  32. 32. Stotz SC, Vriens J, Martyn D, Clardy J, Clapham DE. Citral sensing by Transient [corrected] receptor potential channels in dorsal root ganglion neurons. PLoS ONE 2008; 3(5):e2082. pmid:18461159
  33. 33. Ohkawara S, Tanaka-Kagawa T, Furukawa Y, Nishimura T, Jinno H. Activation of the human transient receptor potential vanilloid subtype 1 by essential oils. Biol. Pharm. Bull. 2010; 33(8):1434–7. pmid:20686244
  34. 34. Manteniotis S, Lehmann R, Flegel C, Vogel F, Hofreuter A, Schreiner BSP, et al. Comprehensive RNA-Seq Expression Analysis of Sensory Ganglia with a Focus on Ion Channels and GPCRs in Trigeminal Ganglia. PLoS ONE 2013; 8(11):e79523. pmid:24260241
  35. 35. Dong X, Han S, Zylka MJ, Simon MI, Anderson DJ. A diverse family of GPCRs expressed in specific subsets of nociceptive sensory neurons. Cell 2001; 106(5):619–32. pmid:11551509
  36. 36. Vandewauw I, Owsianik G, Voets T. Systematic and quantitative mRNA expression analysis of TRP channel genes at the single trigeminal and dorsal root ganglion level in mouse. BMC Neurosci 2013; 14:21. pmid:23410158
  37. 37. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008; 5(7):621–8. pmid:18516045
  38. 38. Ramsköld D, Wang ET, Burge CB, Sandberg R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput. Biol. 2009; 5(12):e1000598. pmid:20011106
  39. 39. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009; 10(1):57–63. pmid:19015660
  40. 40. Garber M, Grabherr MG, Guttman M, Trapnell C. Computational methods for transcriptome annotation and quantification using RNA-seq. Nat. Methods 2011; 8(6):469–77. pmid:21623353
  41. 41. Flegel C, Manteniotis S, Osthold S, Hatt H, Gisselmann G. Expression profile of ectopic olfactory receptors determined by deep sequencing. PLoS ONE 2013; 8(2):e55368. pmid:23405139
  42. 42. Shiao M, Chang AY, Liao B, Ching Y, Lu MJ, Chen SM, et al. Transcriptomes of mouse olfactory epithelium reveal sexual differences in odorant detection. Genome Biol Evol 2012; 4(5):703–12. pmid:22511034
  43. 43. Kanageswaran N, Demond M, Nagel M, Schreiner BSP, Baumgart S, Scholz P, et al. Deep sequencing of the murine olfactory receptor neuron transcriptome. PLoS ONE 2015; 10(1):e0113170. pmid:25590618
  44. 44. Ibarra-Soria X, Levitin MO, Saraiva LR, Logan DW. The olfactory transcriptomes of mice. PLoS Genet. 2014; 10(9):e1004593. pmid:25187969
  45. 45. Hanani M. Satellite glial cells in sensory ganglia: from form to function. Brain Res. Brain Res. Rev. 2005; 48(3):457–76. pmid:15914252
  46. 46. Borowsky B, Adham N, Jones KA, Raddatz R, Artymyshyn R, Ogozalek KL, et al. Trace amines: identification of a family of mammalian G protein-coupled receptors. Proc. Natl. Acad. Sci. U.S.A. 2001; 98(16):8966–71. pmid:11459929
  47. 47. Glusman G, Yanai I, Rubin I, Lancet D. The complete human olfactory subgenome. Genome Res. 2001; 11(5):685–702. pmid:11337468
  48. 48. Firestein S. How the olfactory system makes sense of scents. Nature 2001; 413(6852):211–8. pmid:11557990
  49. 49. Buck L, Axel R. A novel multigene family may encode odorant receptors: a molecular basis for odor recognition. Cell 1991; 65(1):175–87. pmid:1840504
  50. 50. Ben-Arie N, Lancet D, Taylor C, Khen M, Walker N, Ledbetter DH, et al. Olfactory receptor gene cluster on human chromosome 17: possible duplication of an ancestral receptor repertoire. Hum. Mol. Genet. 1994; 3(2):229–35. pmid:8004088
  51. 51. Asai H, Kasai H, Matsuda Y, Yamazaki N, Nagawa F, Sakano H, et al. Genomic structure and transcription of a murine odorant receptor gene: differential initiation of transcription in the olfactory and testicular cells. Biochem. Biophys. Res. Commun. 1996; 221(2):240–7. pmid:8619840
  52. 52. Parmentier M, Libert F, Schurmans S, Schiffmann S, Lefort A, Eggerickx D, et al. Expression of members of the putative olfactory receptor gene family in mammalian germ cells. Nature 1992; 355(6359):453–5. pmid:1370859
  53. 53. Feldmesser E, Olender T, Khen M, Yanai I, Ophir R, Lancet D. Widespread ectopic expression of olfactory receptor genes. BMC Genomics 2006; 7:121. pmid:16716209
  54. 54. Zhang X, La Cruz O de, Pinto JM, Nicolae D, Firestein S, Gilad Y. Characterizing the expression of the human olfactory receptor gene family using a novel DNA microarray. Genome Biol. 2007; 8(5):R86. pmid:17509148
  55. 55. Spehr M, Gisselmann G, Poplawski A, Riffell JA, Wetzel CH, Zimmer RK, et al. Identification of a testicular odorant receptor mediating human sperm chemotaxis. Science 2003; 299(5615):2054–8. pmid:12663925
  56. 56. Braun T, Voland P, Kunz L, Prinz C, Gratzl M. Enterochromaffin cells of the human gut: sensors for spices and odorants. Gastroenterology 2007; 132(5):1890–901. pmid:17484882
  57. 57. Neuhaus EM, Zhang W, Gelis L, Deng Y, Noldus J, Hatt H. Activation of an olfactory receptor inhibits proliferation of prostate cancer cells. J. Biol. Chem. 2009; 284(24):16218–25. pmid:19389702
  58. 58. Gu X, Karp PH, Brody SL, Pierce RA, Welsh MJ, Holtzman MJ, et al. Volatile-Sensing Functions for Pulmonary Neuroendocrine Cells. Am. J. Respir. Cell Mol. Biol. 2013.
  59. 59. Busse D, Kudella P, Grüning N, Gisselmann G, Ständer S, Luger T, et al. A Synthetic Sandalwood Odorant Induces Wound Healing Processes in Human Keratinocytes via the Olfactory Receptor OR2AT4. J. Invest. Dermatol. 2014.
  60. 60. Maßberg D, Simon A, Häussinger D, Keitel V, Gisselmann G, Conrad H, et al. Monoterpene (-)-citronellal affects hepatocarcinoma cell signaling via an olfactory receptor. Arch. Biochem. Biophys. 2015; 566:100–9. pmid:25513961
  61. 61. Verbeurgt C, Wilkin F, Tarabichi M, Gregoire F, Dumont JE, Chatelain P. Profiling of olfactory receptor gene expression in whole human olfactory mucosa. PLoS ONE 2014; 9(5):e96333. pmid:24800820
  62. 62. Keydar I, Ben-Asher E, Feldmesser E, Nativ N, Oshimoto A, Restrepo D, et al. General olfactory sensitivity database (GOSdb): candidate genes and their genomic variations. Hum. Mutat. 2013; 34(1):32–41. pmid:22936402
  63. 63. Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief. Bioinformatics 2013; 14(2):178–92. pmid:22517427
  64. 64. Veitinger T, Riffell JR, Veitinger S, Nascimento JM, Triller A, Chandsawangbhuwana C, et al. Chemosensory Ca2+ dynamics correlate with diverse behavioral phenotypes in human sperm. J. Biol. Chem. 2011; 286(19):17311–25. pmid:21454470
  65. 65. Hartmann C, Triller A, Spehr M, Dittrich R, Hatt H, Buettner A. Sperm-Activating Odorous Substances in Human Follicular Fluid and Vaginal Secretion: Identification by Gas Chromatography–Olfactometry and Ca2+ Imaging. ChemPlusChem 2013; 78(7):695–702.
  66. 66. Mainland JD, Keller A, Li YR, Zhou T, Trimmer C, Snyder LL, et al. The missense of smell: functional variability in the human odorant receptor repertoire. Nat. Neurosci. 2014; 17(1):114–20. pmid:24316890
  67. 67. Grosser BI, Monti-Bloch L, Jennings-White C, Berliner DL. Behavioral and electrophysiological effects of androstadienone, a human pheromone. Psychoneuroendocrinology 2000; 25(3):289–99. pmid:10737699
  68. 68. Savic I, Berglund H, Gulyas B, Roland P. Smelling of odorous sex hormone-like compounds causes sex-differentiated hypothalamic activations in humans. Neuron 2001; 31(4):661–8. pmid:11545724
  69. 69. Hummel T, Schultz S, Witt M, Hatt H. Electrical responses to chemosensory stimulation recorded from the vomeronasal duct and the respiratory epithelium in humans. Int J Psychophysiol 2011; 81(2):116–20. pmid:21619899
  70. 70. Abram SE, Yi J, Fuchs A, Hogan QH. Permeability of injured and intact peripheral nerves and dorsal root ganglia. Anesthesiology 2006; 105(1):146–53. pmid:16810006
  71. 71. Ichikawa H, Sugimoto T. Vanilloid receptor 1-like receptor-immunoreactive primary sensory neurons in the rat trigeminal nervous system. Neuroscience 2000; 101(3):719–25. pmid:11113320
  72. 72. Devor M. Unexplained preculiarities of the dorsal root ganglion. Pain Supplement 6 1999:27–35.
  73. 73. Liberles SD, Buck LB. A second class of chemosensory receptors in the olfactory epithelium. Nature 2006; 442(7103):645–50. pmid:16878137
  74. 74. Liberles SD. Trace amine-associated receptors are olfactory receptors in vertebrates. Ann. N. Y. Acad. Sci. 2009; 1170:168–72. pmid:19686131
  75. 75. Ferrero DM, Lemon JK, Fluegge D, Pashkovski SL, Korzan WJ, Datta SR, et al. Detection and avoidance of a carnivore odor by prey. Proc. Natl. Acad. Sci. U.S.A. 2011; 108(27):11235–40. pmid:21690383
  76. 76. Wallrabenstein I, Kuklan J, Weber L, Zborala S, Werner M, Altmüller J, et al. Human trace amine-associated receptor TAAR5 can be activated by trimethylamine. PLoS ONE 2013; 8(2):e54950. pmid:23393561
  77. 77. Rodriguez I, Greer CA, Mok MY, Mombaerts P. A putative pheromone receptor gene expressed in human olfactory mucosa. Nat. Genet. 2000; 26(1):18–9. pmid:10973240
  78. 78. Shirokova E, Raguse JD, Meyerhof W, Krautwurst D. The human vomeronasal type-1 receptor family—detection of volatiles and cAMP signaling in HeLa/Olf cells. FASEB J. 2008; 22(5):1416–25. pmid:18096815
  79. 79. Scott K. Taste recognition: food for thought. Neuron 2005; 48(3):455–64. pmid:16269362
  80. 80. Liu L, Simon SA. Responses of cultured rat trigeminal ganglion neurons to bitter tastants. Chem. Senses 1998; 23(2):125–30. pmid:9589160
  81. 81. Migeotte I, Communi D, Parmentier M. Formyl peptide receptors: a promiscuous subfamily of G protein-coupled receptors controlling immune responses. Cytokine Growth Factor Rev. 2006; 17(6):501–19. pmid:17084101
  82. 82. Rivière S, Challet L, Fluegge D, Spehr M, Rodriguez I. Formyl peptide receptor-like proteins are a novel family of vomeronasal chemosensors. Nature 2009; 459(7246):574–7. pmid:19387439
  83. 83. Tatemoto K, Nozaki Y, Tsuda R, Konno S, Tomura K, Furuno M, et al. Immunoglobulin E-independent activation of mast cell is mediated by Mrg receptors. Biochem. Biophys. Res. Commun. 2006; 349(4):1322–8. pmid:16979137
  84. 84. Subramanian H, Kashem SW, Collington SJ, Qu H, Lambris JD, Ali H. PMX-53 as a dual CD88 antagonist and an agonist for Mas-related gene 2 (MrgX2) in human mast cells. Mol. Pharmacol. 2011; 79(6):1005–13. pmid:21441599
  85. 85. Han S, Dong X, Hwang J, Zylka MJ, Anderson DJ, Simon MI. Orphan G protein-coupled receptors MrgA1 and MrgC11 are distinctively activated by RF-amide-related peptides through the Galpha q/11 pathway. Proc. Natl. Acad. Sci. U.S.A. 2002; 99(23):14740–5. pmid:12397184
  86. 86. Shinohara T, Harada M, Ogi K, Maruyama M, Fujii R, Tanaka H, et al. Identification of a G protein-coupled receptor specifically responsive to beta-alanine. J. Biol. Chem. 2004; 279(22):23559–64. pmid:15037633
  87. 87. Rau KK, McIlwrath SL, Wang H, Lawson JJ, Jankowski MP, Zylka MJ, et al. Mrgprd enhances excitability in specific populations of cutaneous murine polymodal nociceptors. J. Neurosci. 2009; 29(26):8612–9. pmid:19571152
  88. 88. Vrontou S, Wong AM, Rau KK, Koerber HR, Anderson DJ. Genetic identification of C fibres that detect massage-like stroking of hairy skin in vivo. Nature 2013; 493(7434):669–73. pmid:23364746
  89. 89. Solinski HJ, Zierler S, Gudermann T, Breit A. Human sensory neuron-specific Mas-related G protein-coupled receptors-X1 sensitize and directly activate transient receptor potential cation channel V1 via distinct signaling pathways. J. Biol. Chem. 2012; 287(49):40956–71. pmid:23074220
  90. 90. Bautista DM, Sigal YM, Milstein AD, Garrison JL, Zorn JA, Tsuruda PR, et al. Pungent agents from Szechuan peppers excite sensory neurons by inhibiting two-pore potassium channels. Nat. Neurosci. 2008; 11(7):772–9. pmid:18568022
  91. 91. Beltrán LR, Dawid C, Beltrán M, Gisselmann G, Degenhardt K, Mathie K, et al. The pungent substances piperine, capsaicin, 6-gingerol and polygodial inhibit the human two-pore domain potassium channels TASK-1, TASK-3 and TRESK. Front Pharmacol 2013; 4:141. pmid:24302912
  92. 92. Schöbel N, Radtke D, Lübbert M, Gisselmann G, Lehmann R, Cichy A, et al. Trigeminal ganglion neurons of mice show intracellular chloride accumulation and chloride-dependent amplification of capsaicin-induced responses. PLoS ONE 2012; 7(11):e48005. pmid:23144843
  93. 93. Clapham DE, Runnels LW, Strübing C. The TRP ion channel family. Nat. Rev. Neurosci. 2001; 2(6):387–96. pmid:11389472
  94. 94. Flockerzi V. An introduction on TRP channels. Handb Exp Pharmacol 2007; (179):1–19. pmid:17217048
  95. 95. Patapoutian A, Peier AM, Story GM, Viswanath V. ThermoTRP channels and beyond: mechanisms of temperature sensation. Nat. Rev. Neurosci. 2003; 4(7):529–39. pmid:12838328
  96. 96. Caterina MJ, Rosen TA, Tominaga M, Brake AJ, Julius D. A capsaicin-receptor homologue with a high threshold for noxious heat. Nature 1999; 398(6726):436–41. pmid:10201375
  97. 97. Liedtke W, Choe Y, Martí-Renom MA, Bell AM, Denis CS, Sali A, et al. Vanilloid receptor-related osmotically activated channel (VR-OAC), a candidate vertebrate osmoreceptor. Cell 2000; 103(3):525–35. pmid:11081638
  98. 98. Strotmann R, Harteneck C, Nunnenmacher K, Schultz G, Plant TD. OTRPC4, a nonselective cation channel that confers sensitivity to extracellular osmolarity. Nat. Cell Biol. 2000; 2(10):695–702. pmid:11025659
  99. 99. Peier AM, Moqrich A, Hergarden AC, Reeve AJ, Andersson DA, Story GM, et al. A TRP channel that senses cold stimuli and menthol. Cell 2002; 108(5):705–15. pmid:11893340
  100. 100. Smith GD, Gunthorpe MJ, Kelsell RE, Hayes PD, Reilly P, Facer P, et al. TRPV3 is a temperature-sensitive vanilloid receptor-like protein. Nature 2002; 418(6894):186–90. pmid:12077606
  101. 101. Xu H, Ramsey IS, Kotecha SA, Moran MM, Chong JA, Lawson D, et al. TRPV3 is a calcium-permeable temperature-sensitive cation channel. Nature 2002; 418(6894):181–6. pmid:12077604
  102. 102. van Haute C, Ridder D de, Nilius B. TRP channels in human prostate. ScientificWorldJournal 2010; 10:1597–611. pmid:20730379
  103. 103. Faussone-Pellegrini M, Taddei A, Bizzoco E, Lazzeri M, Vannucchi MG, Bechi P. Distribution of the vanilloid (capsaicin) receptor type 1 in the human stomach. Histochem. Cell Biol. 2005; 124(1):61–8. pmid:16041554
  104. 104. Dömötör A, Peidl Z, Vincze A, Hunyady B, Szolcsányi J, Kereskay L, et al. Immunohistochemical distribution of vanilloid receptor, calcitonin-gene related peptide and substance P in gastrointestinal mucosa of patients with different gastrointestinal disorders. Inflammopharmacology 2005; 13(1–3):161–77. pmid:16354389
  105. 105. Guarino MPL, Cheng L, Ma J, Harnett K, Biancani P, Altomare A, et al. Increased TRPV1 gene expression in esophageal mucosa of patients with non-erosive and erosive reflux disease. Neurogastroenterol. Motil. 2010; 22(7):746–51, e219. pmid:20456759
  106. 106. Martínez-García MC, Martínez T, Pañeda C, Gallego P, Jimenez AI, Merayo J. Differential expression and localization of transient receptor potential vanilloid 1 in rabbit and human eyes. Histol. Histopathol. 2013; 28(11):1507–16. pmid:23709255
  107. 107. Mezey E, Tóth ZE, Cortright DN, Arzubi MK, Krause JE, Elde R, et al. Distribution of mRNA for vanilloid receptor subtype 1 (VR1), and VR1-like immunoreactivity, in the central nervous system of the rat and human. Proc. Natl. Acad. Sci. U.S.A. 2000; 97(7):3655–60. pmid:10725386
  108. 108. Szallasi A. Autoradiographic visualization and pharmacological characterization of vanilloid (capsaicin) receptors in several species, including man. Acta Physiol Scand Suppl 1995; 629:1–68. pmid:8801775
  109. 109. Biggs JE, Yates JM, Loescher AR, Clayton NM, Boissonade FM, Robinson PP. Vanilloid receptor 1 (TRPV1) expression in lingual nerve neuromas from patients with or without symptoms of burning pain. Brain Res. 2007; 1127(1):59–65. pmid:17109831
  110. 110. Morgan CR, Bird EV, Robinson PP, Boissonade FM. TRPA1 expression in human lingual nerve neuromas in patients with and without symptoms of dysaesthesia. Neurosci. Lett. 2009; 465(2):189–93. pmid:19715741
  111. 111. Yilmaz Z, Renton T, Yiangou Y, Zakrzewska J, Chessell IP, Bountra C, et al. Burning mouth syndrome as a trigeminal small fibre neuropathy: Increased heat and capsaicin receptor TRPV1 in nerve fibres correlates with pain score. J Clin Neurosci 2007; 14(9):864–71. pmid:17582772
  112. 112. Seki N, Shirasaki H, Kikuchi M, Sakamoto T, Watanabe N, Himi T. Expression and localization of TRPV1 in human nasal mucosa. Rhinology 2006; 44(2):128–34. pmid:16792172
  113. 113. Anand U, Otto WR, Facer P, Zebda N, Selmer I, Gunthorpe MJ, et al. TRPA1 receptor localisation in the human peripheral nervous system and functional studies in cultured human and rat sensory neurons. Neurosci. Lett. 2008; 438(2):221–7. pmid:18456404
  114. 114. Axelsson HE, Minde JK, Sonesson A, Toolanen G, Högestätt ED, Zygmunt PM. Transient receptor potential vanilloid 1, vanilloid 2 and melastatin 8 immunoreactive nerve fibers in human skin from individuals with and without Norrbottnian congenital insensitivity to pain. Neuroscience 2009; 162(4):1322–32. pmid:19482060
  115. 115. Facer P, Casula MA, Smith GD, Benham CD, Chessell IP, Bountra C, et al. Differential expression of the capsaicin receptor TRPV1 and related novel receptors TRPV3, TRPV4 and TRPM8 in normal human tissues and changes in traumatic and diabetic neuropathy. BMC Neurol 2007; 7:11. pmid:17521436
  116. 116. Pecze L, Szabó K, Széll M, Jósvay K, Kaszás K, Kúsz E, et al. Human keratinocytes are vanilloid resistant. PLoS ONE 2008; 3(10):e3419. pmid:18852901
  117. 117. Hou M, Uddman R, Tajti J, Kanje M, Edvinsson L. Capsaicin receptor immunoreactivity in the human trigeminal ganglion. Neurosci. Lett. 2002; 330(3):223–6. pmid:12270633
  118. 118. Mandadi S, Sokabe T, Shibasaki K, Katanosaka K, Mizuno A, Moqrich A, et al. TRPV3 in keratinocytes transmits temperature information to sensory neurons via ATP. Pflugers Arch. 2009; 458(6):1093–102. pmid:19669158
  119. 119. Suzuki M, Watanabe Y, Oyama Y, Mizuno A, Kusano E, Hirao A, et al. Localization of mechanosensitive channel TRPV4 in mouse skin. Neurosci. Lett. 2003; 353(3):189–92. pmid:14665413
  120. 120. Gopinath P, Wan E, Holdcroft A, Facer P, Davis JB, Smith GD, et al. Increased capsaicin receptor TRPV1 in skin nerve fibres and related vanilloid receptors TRPV3 and TRPV4 in keratinocytes in human breast pain. BMC Womens Health 2005; 5(1):2. pmid:15755319
  121. 121. Fusi C, Materazzi S, Minocci D, Maio V, Oranges T, Massi D, et al. Transient Receptor Potential Vanilloid 4 (TRPV4) Is Downregulated in Keratinocytes in Human Non-Melanoma Skin Cancer. J. Invest. Dermatol. 2014; 134(9):2408–17. pmid:24643128
  122. 122. Ishimaru Y, Inada H, Kubota M, Zhuang H, Tominaga M, Matsunami H. Transient receptor potential family members PKD1L3 and PKD2L1 form a candidate sour taste receptor. Proc. Natl. Acad. Sci. U.S.A. 2006; 103(33):12569–74. pmid:16891422
  123. 123. Rasche S, Toetter B, Adler J, Tschapek A, Doerner JF, Kurtenbach S, et al. Tmem16b is specifically expressed in the cilia of olfactory sensory neurons. Chem. Senses 2010; 35(3):239–45. pmid:20100788
  124. 124. Pifferi S, Pascarella G, Boccaccio A, Mazzatenta A, Gustincich S, Menini A, et al. Bestrophin-2 is a candidate calcium-activated chloride channel involved in olfactory transduction. Proc. Natl. Acad. Sci. U.S.A. 2006; 103(34):12929–34. pmid:16912113
  125. 125. Cho H, Yang YD, Lee J, Lee B, Kim T, Jang Y, et al. The calcium-activated chloride channel anoctamin 1 acts as a heat sensor in nociceptive neurons. Nat. Neurosci. 2012; 15(7):1015–21. pmid:22634729
  126. 126. Duran C, Qu Z, Osunkoya AO, Cui Y, Hartzell HC. ANOs 3–7 in the anoctamin/Tmem16 Cl- channel family are intracellular proteins. Am. J. Physiol., Cell Physiol. 2012; 302(3):C482–93. pmid:22075693
  127. 127. Kunzelmann K, Schreiber R, Kmit A, Jantarajit W, Martins JR, Faria D, et al. Expression and function of epithelial anoctamins. Exp. Physiol. 2012; 97(2):184–92. pmid:21908539
  128. 128. Liu B, Linley JE, Du X, Zhang X, Ooi L, Zhang H, et al. The acute nociceptive signals induced by bradykinin in rat sensory neurons are mediated by inhibition of M-type K+ channels and activation of Ca2+-activated Cl- channels. J. Clin. Invest. 2010; 120(4):1240–52. pmid:20335661
  129. 129. Jasti J, Furukawa H, Gonzales EB, Gouaux E. Structure of acid-sensing ion channel 1 at 1.9 A resolution and low pH. Nature 2007; 449(7160):316–23. pmid:17882215
  130. 130. Gonzales EB, Kawate T, Gouaux E. Pore architecture and ion sites in acid-sensing ion channels and P2X receptors. Nature 2009; 460(7255):599–604. pmid:19641589
  131. 131. Yiangou Y, Facer P, Smith JA, Sangameswaran L, Eglen R, Birch R, et al. Increased acid-sensing ion channel ASIC-3 in inflamed human intestine. Eur J Gastroenterol Hepatol 2001; 13(8):891–6. pmid:11507351
  132. 132. Li W, Xu T. ASIC3 channels in multimodal sensory perception. ACS Chem Neurosci 2011; 2(1):26–37. pmid:22778854
  133. 133. Coste B, Mathur J, Schmidt M, Earley TJ, Ranade S, Petrus MJ, et al. Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science 2010; 330(6000):55–60. pmid:20813920
  134. 134. Burnstock G. Purinergic receptors and pain. Curr. Pharm. Des. 2009; 15(15):1717–35. pmid:19442186
  135. 135. Sondersorg AC, Busse D, Kyereme J, Rothermel M, Neufang G, Gisselmann G, et al. Chemosensory information processing between keratinocytes and trigeminal neurons. J. Biol. Chem. 2014.
  136. 136. Pan A, Wu H, Li M, Lu D, He X, Yi X, et al. Prenatal expression of purinergic receptor P2X3 in human dorsal root ganglion. Purinergic Signal. 2012; 8(2):245–54. pmid:22052556
  137. 137. Gerhold KA, Bautista DM. Molecular and cellular mechanisms of trigeminal chemosensation. Ann. N. Y. Acad. Sci. 2009; 1170:184–9. pmid:19686135
  138. 138. Takeda M, Tsuboi Y, Kitagawa J, Nakagawa K, Iwata K, Matsumoto S. Potassium channels as a potential therapeutic target for trigeminal neuropathic and inflammatory pain. Mol Pain 2011; 7:5. pmid:21219657
  139. 139. Gutman GA, Chandy KG, Grissmer S, Lazdunski M, McKinnon D, Pardo LA, et al. International Union of Pharmacology. LIII. Nomenclature and molecular relationships of voltage-gated potassium channels. Pharmacol. Rev. 2005; 57(4):473–508. pmid:16382104
  140. 140. Schreiner BSP, Lehmann R, Thiel U, Ziemba PM, Beltrán L, Sherkheli MA, et al. Direct action and modulating effect of (+)- and (-)-nicotine on ion channels expressed in trigeminal sensory neurons. Eur. J. Pharmacol. 2014.
  141. 141. Medhurst AD, Rennie G, Chapman CG, Meadows H, Duckworth MD, Kelsell RE, et al. Distribution analysis of human two pore domain potassium channels in tissues of the central nervous system and periphery. Brain Res. Mol. Brain Res. 2001; 86(1–2):101–14. pmid:11165386
  142. 142. Lafrenière RG, Cader MZ, Poulin J, Andres-Enguix I, Simoneau M, Gupta N, et al. A dominant-negative mutation in the TRESK potassium channel is linked to familial migraine with aura. Nat. Med. 2010; 16(10):1157–60. pmid:20871611
  143. 143. Enyedi P, Czirják G. Molecular background of leak K+ currents: two-pore domain potassium channels. Physiol. Rev. 2010; 90(2):559–605. pmid:20393194
  144. 144. Johnson KW, Phebus LA, Cohen ML. Serotonin in migraine: theories, animal models and emerging therapies. Prog Drug Res 1998; 51:219–44. pmid:9949863
  145. 145. Coughlin SR, Camerer E. PARticipation in inflammation. J. Clin. Invest. 2003; 111(1):25–7. pmid:12511583
  146. 146. Gabashvili I, Sokolowski B, Morton C., Giersch A. Ion Channel Gene Expression in the Inner Ear. J Assoc Res Otolaryngol 2007; 8(3):305–28. pmid:17541769
  147. 147. Catterall WA. From ionic currents to molecular mechanisms: the structure and function of voltage-gated sodium channels. Neuron 2000; 26(1):13–25. pmid:10798388
  148. 148. Dib-Hajj SD, Cummins TR, Black JA, Waxman SG. Sodium channels in normal and pathological pain. Annu. Rev. Neurosci. 2010; 33:325–47. pmid:20367448
  149. 149. Strickland IT, Martindale JC, Woodhams PL, Reeve AJ, Chessell IP, McQueen DS. Changes in the expression of NaV1.7, NaV1.8 and NaV1.9 in a distinct population of dorsal root ganglia innervating the rat knee joint in a model of chronic inflammatory joint pain. Eur J Pain 2008; 12(5):564–72. pmid:17950013
  150. 150. Padilla F, Couble M, Coste B, Maingret F, Clerc N, Crest M, et al. Expression and localization of the Nav1.9 sodium channel in enteric neurons and in trigeminal sensory endings: implication for intestinal reflex function and orofacial pain. Mol. Cell. Neurosci. 2007; 35(1):138–52. pmid:17363266
  151. 151. 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. pmid:21220098
  152. 152. Kilic G, Ismail Kucukali C, Orhan N, Ozkok E, Zengin A, Aydin M, et al. Are GRIK3 (T928G) gene variants in schizophrenia patients different from those in their first-degree relatives? Psychiatry Res 2010; 175(1–2):43–6. pmid:20045197
  153. 153. Bönigk W, Bradley J, Müller F, Sesti F, Boekhoff I, Ronnett GV, et al. The native rat olfactory cyclic nucleotide-gated channel is composed of three distinct subunits. J. Neurosci. 1999; 19(13):5332–47. pmid:10377344
  154. 154. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010; 28(5):511–5. pmid:20436464
  155. 155. Ameur A, Zaghlool A, Halvardson J, Wetterbom A, Gyllensten U, Cavelier L, et al. Total RNA sequencing reveals nascent transcription and widespread co-transcriptional splicing in the human brain. Nat. Struct. Mol. Biol. 2011; 18(12):1435–40. pmid:22056773