|
|
||||||||
1 Johnson and Johnson Pharmaceutical Research and Development, Department of Internal Medicine, Beerse, Belgium
2 Holburn Group of Companies, Bowmanville, Ontario, Canada
3 Johnson and Johnson Pharmaceutical Research and Development, Department of Functional Genomics, Beerse, Belgium
| ABSTRACT |
|---|
|
|
|---|
microarray; vagal and spinal afferents; nociception; spectral map analysis
| INTRODUCTION |
|---|
|
|
|---|
Vagal and spinal afferents supplying the GI tract also differ in the pattern of their terminal innervation, which in part determines the stimulus-response properties of the peripheral endings (4). Vagal afferents terminate close to the mucosal epithelium, where they are exposed to chemicals (e.g., nutrients) absorbed from the lumen or mediators released from enteroendocrine cells or immune cells in the lamina propria. Vagal afferents also form intramuscular arrays and intraganglionic laminar endings that are thought to detect mechanical activity. Spinal afferents also innervate the mucosa, submucosa, and myenteric plexus. Additionally, projections of DRG neurons terminate in the serosa and mesenteric attachments, often in association with blood vessels. These endings are mechanosensitive, but the basis of this mechanosensitivity at the molecular level is unknown. Both vagal and spinal afferents respond to distension and contraction, but while vagal afferent endings respond to levels of distension that occur during the normal course of digestion, many spinal afferents have thresholds for activation that when applied in humans give rise to discomfort or pain (20).
These observations are the basis for the common view that vagal and spinal afferents have different functional roles: spinal afferents play a major role in nociception, while vagal afferents mediate physiological responses and behavioral regulation, particularly in relation to food intake, satiety, anorexia, and emesis. However, vagal and spinal afferents share a number of features: both have a large proportion of unmyelinated axons that can be ablated by capsaicin, and both express the capsaicin receptor (TRPV1), considered a hallmark of nociceptive neurons (44). Both NG and DRG neurons have been shown to become sensitized after inflammation, demonstrating plasticity in the mechanisms that regulate neuronal excitability, which has implications for pain processing (14).
Studies employing PCR, in situ hybridization, and immunohistochemistry have yielded some insight into the extent to which regulatory molecules form common pathways for the processing of GI sensory information, and those that are unique to vagal or spinal afferents. However, the advent of expression profiling and the availability of microarray chips that cover the mouse genome facilitate the exciting opportunity to examine the whole range of receptors, channels, transporters, and kinases that determine and shape sensory signal transduction. Here we have compared the expression profiles of NG and DRG neurons projecting to the mouse peritoneal cavity.
| MATERIAL AND METHODS |
|---|
|
|
|---|
Assessment of Numbers of Neurons in Sensory Ganglia
In a second of series of mice (n = 4), CTB488 was administered ip (250 µg in 100 µl), and the animals were euthanized 4 days later. One of each pair of DRG (from T10 to T13) or NG was frozen before sectioning on the cryostat at 10 µm. The paired ganglia from the contralateral side were squashed on slides beneath coverslips, as described above. Photomicrographs of at least 10 sections per ganglion and the squash preparations were prepared, and the numbers of fluorescent cells were counted from the resultant photomicrographs. Subsequently, the sections were stained with methylene blue, and the total numbers of neurons (ganglion cells containing a recognizable nucleus) were counted. From these measurements, the percentage of fluorescent neurons in the sections (no. of fluorescent neurons x 100/total no. of neurons) was determined.
Laser-Capture Microdissection
All ganglia for microarray studies were isolated 34 days after one ip injection of CTB488. NG and DRG (T10T13) were harvested, placed in tissue-freezing medium (TFM; Triangle Biomedical Sciences, Durham, NC), frozen, and stored at 80°C.
Cryostat sections (12 µm) were attached to RNase-free PEN membrane-covered glass slides (PALM Microlaser Technologies, Bernried, Germany), fixed with 100% ethanol and air dried before laser-capture microdissection (LCM). Microdissection was performed on a PALM microbeam-equipped microscope (Axiovert 135; Zeiss, Gottingen, Germany). Fluorescent cells were excised after Nissl staining [0.5% cresyl violet acetate (Sigma-Aldrich, St. Louis, MO)/0.1 M sodium acetate (Fluka, Buchs, Switzerland)]. For each animal, a total of 200300 cells were captured: 100150 from DRG and 100150 from NG. The microdissected cells from the DRG and the NG, respectively, were catapulted in two separate tubes, each one filled with 75 µl of RNeasy lysis buffer (RLT; Qiagen, Hilden, Germany) containing 0.14 M ß-mercaptoethanol and 200 ng polyinosinic acid (Sigma).
RNA Isolation
Laser-captured samples were incubated at 42°C for 20 min and then chilled on ice. An equal volume of 70% ethanol was added to each sample and then transferred to RNeasy MinElute Spin Columns (Qiagen). RNA was cleaned up according to the manufacturer's instructions, eluted in 14 µl of RNase-free water, and adjusted to 4 µl by vacuum drying.
RNA Amplification
As "spike-in" controls, the GeneChip Poly-A RNA control kit (Affymetrix, Santa Clara, CA) was used. Serial dilutions were made of the prokaryotic Poly-A control using the following dilution steps: 1:20, 1:50, 1:50, 1:20, and 1:10. First-strand cDNA was prepared as described by the Affymetrix two-cycle cDNA synthesis protocol, except for the use of Superscript III (Invitrogen, Carlsbad, CA) and incubation at 50°C for 30 min. Second-strand master mix consisted of 1 µl of 10x Bst polymerase buffer (Epicentre, Madison, WI), 1 µl of 10 mM dNTP (Invitrogen), 0.5 µl (1 U) thermostable RNaseH (Invitrogen), 1 µl (5 U) Bst DNA polymerase (Epicentre), and water to 10 µl. This master mix was added to the first-strand cDNA reaction and incubated at 65°C for 10 min before heat inactivation at 80°C for 3 min. Subsequently 2 µl of exonuclease mix were added containing ExoI and ExoVII and incubated at 37°C for 10 min, followed by heat inactivation at 80°C for 3 min. Double-stranded cDNA was transcribed at 42°C for 3 h using the AmpliScribe T7 High Yield Transcription Kit (Epicentre) in a total volume of 100 µl. The resulting amplified RNA was incubated with DNase I (4 U/µl) at 37°C for 15 min. Amplified RNA was purified after addition of 100 ng of polyinosinic acid using RNeasy MinElute Cleanup Kit (Qiagen). RNA was eluted in 14 µl of RNase-free water and adjusted to 4 µl by vacuum drying. The second round of amplification was performed as described above, except that 50 ng of random hexamer primers were used to prime the reverse transcription reaction and the second-strand cDNA reaction was primed with 0.25 ng of T7 oligonucleotide.
RNA Labeling and Microarray Hybridization
The third round of amplification, including biotin labeling, was performed on 500 ng of second round-amplified RNA. First-strand cDNA synthesis was performed as described above. Second-strand cDNA synthesis was performed using 1 µl of T7 oligonucleotide dT24 (Affymetrix, 100 pmol/µl) annealed for 5 min at 70°C, and the reaction was then incubated at 42°C for 10 min. A master mix was prepared, consisting of 10x second-strand buffer, dNTPs (200 mM final), Escherichia coli RNaseH (2 U), and 10 U of E. coli DNA polymerase (Invitrogen), and added to the first-strand reaction to obtain a 50-µl reaction volume. After incubation at 37°C for 10 min, denaturation was done at 80°C for 3 min. Cleanup was performed with Qiagen PCR purification kit. Synthesis of biotin-labeled RNA was performed using the BioArray High Yield RNA Transcript Labeling Kit (Enzo Life Sciences, Farmingdale, NY). Cleanup was done using the RNeasy Mini Kit (Qiagen). Labeled RNA was hybridized to either mouse genome MG-U74Av2 or MG-430_2.0 GeneChip arrays (Affymetrix). Hybridization of microarrays was performed using 12.5 µg of biotin-labeled RNA at 45°C for 16 h under continuous rotation. Arrays were stained in Affymetrix Fluidics stations using streptavidin-phycoerythrin (SAPE) followed by staining with anti-streptavidin antibody and a second SAPE staining. Subsequently, arrays were scanned with a Agilent Laserscanner (Affymetrix), and data were analyzed with the Microarray Suite (MAS) software version 5.0 (Affymetrix). No scaling or normalization was performed at this stage. All data are available at the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) with accession number GSE2917.
Data Analysis and Selection of Genes
Normalization.
Genes that were called absent according to Affymetrix MAS 5.0 software (P value > 0.06) were removed from further analysis. Raw intensities were log2 transformed to get data normally distributed. Subsequently, all data were quantile normalized per type of ganglion (1). Following the groupwise quantile normalization, a second quantile normalization was carried out across both DRG and NG samples. This alignment sets the range of intensities of one array to the range measured across all arrays, compensating for array-to-array variations in hybridization, washing and staining, ultimately allowing a reasonable comparison between arrays.
Spectral map analysis.
This unsupervised multivariate projection method, applied to the normalized data set, helps to reduce the complexity of highly dimensional data (n genes vs. p samples) (46) and provides an unbiased means to visually identify the predominant clusters of genes or subjects in the data. The aim of the technique is to retrieve the most predominant differences in the data set, disregarding genes that do not contribute to the difference.
Significance analysis.
Individual genes with different expression levels between groups (NG vs. DRG) were identified using significance analysis of microarray data (SAM) (43), and the false discovery rate was determined using q-values (40). Whereas a P value is commonly used to evaluate a single significance test, the q-value allows for multiple tests performed simultaneously, as in microarray experiments. We applied a 1% threshold (q = 0.01) for our analysis.
Fold-difference filtering.
For comparison purposes, an arbitrary fold-difference (FD) filter was applied excluding all genes that exhibited a difference in expression below 50% (1.5 FD filter).
Effect of CTB488 Injection and Amplification on Gene Expression
Effect of CTB488.
The effect of CTB488 labeling on gene expression profiles was assessed by comparing ganglia isolated from three vehicle-treated animals with three combined intradermal- and ip-injected mice (labeling most neurons). Although a clear difference in expression profile was observed between NG and DRG (see below), no significant effect of the dye injection was noted (data not shown).
Effect of amplification.
Efficiency and sensitivity of amplification were assessed with spike-in controls. In agreement with previous reports, spike-in controls revealed a detection limit of 1 copy in 1,000,000 and a direct correlation between signal intensity and copy number (data not shown). Arrays that did not show this correlation were excluded from data analysis.
Quantitative RT-PCR
Microarray data were confirmed using real-time PCR analysis. First-strand cDNA synthesis was performed on 50 ng of second round-amplified RNA using random hexamer primers and SuperscriptII RT (Invitrogen). Quantitative PCR (qPCR) was performed on an ABI Prism 7900 cycler (Applied Biosystems, Foster City, CA) using a Taqman PCR kit (Applied Biosystems). Serial dilutions of cDNA were used to generate standard curves of threshold cycles vs. the logarithms of concentration for ATPSase and the genes of interest (see Table 1 for sequences of primers; Eurogentec, Seraing, Belgium).
|
| RESULTS |
|---|
|
|
|---|
|
Three separate microarray experiments were performed on RNA isolated from 1) whole ganglia on MG-U74Av2 GeneChips, 2) laser-captured neurons on MG-U74Av2 GeneChips, and 3) laser-captured neurons on MG-430v2.0 GeneChips.
1) Whole ganglia using MG-U74Av2 arrays.
RNA extracted from the entire DRG and NG was hybridized to MG-U74Av2 GeneChips. Graphical exploration of microarray data using spectral map analysis (SPM) revealed overt differences in gene expression between whole DRG and NG (Fig. 2). The first component of the principal component analysis explained 60% of the total variance present in the data set of 7,744 genes that were reliably detected, and this separated DRG from NG samples. To identify differentially expressed individual genes, SAM (q-value < 0.01) and FD filtering (>1.5 FD) were applied. Positioning of the genes at the extremities of the SPM biplot was also taken into account (data not shown). On the basis of these three criteria, 628 genes were identified to be differentially expressed in the entire DRG vs. NG, including 17 G protein-coupled receptors (GPCRs) and 18 ion channels (Table 2).
|
|
Of the 4,956 genes that showed no differences in expression between whole NG and DRG (i.e., FD < 1.5 or q-value > 0.01), 467 were significantly different when measured in laser-captured cells. Striking examples of these differences include galanin, TrkC, and the
2c-adrenergic receptor (Fig. 3). Conversely, for some genes, such as the tetrodotoxin (TTX)-resistant, voltage-gated, Nav1.8 sodium channel, differences determined in whole ganglia were not found when comparing laser-captured NG and DRG neurons. In general, of 5,031 transcripts that showed no differences in expression in laser-captured cells, 579 were significantly different when measured in whole NG and DRG.
|
As a proof of concept, the expression levels of several genes known to be differentially expressed in spinal vs. vagal visceral afferents were examined in the array data of the third experiment and confirmed by qPCR. Of the several 5-HT receptors present on the arrays, only mRNA for the 5-HT3A receptor was readily detected in both NG and DRG neurons (Fig. 4), with a ninefold higher abundance in NG. qPCR revealed a greater difference in expression (21-fold difference). Messenger RNA for the 5-HT3B receptor was clearly expressed in NG but not in DRG neurons, in contrast with mRNA levels for the 5-HT1D receptor, which were markedly higher in DRG neurons than in NG cells. Cholecystokinin (CCK)A receptor mRNA was also detected more abundantly in vagal compared with spinal afferents (Fig. 4). In total, 18 genes were evaluated for FD in expression level using qPCR and compared with array experiments (Table 3). In general, FD was higher with qPCR compared with array data. Very high differences (>8-fold) seem to level off in array experiments. However, a good correlation was present between FD measured by microarray and qPCR (Fig. 5).
|
|
|
-subunits were reliably detected (Fig. 6). Messenger RNA for Nav1.8 (Scn10a) was equally present in vagal and spinal afferents, whereas Nav1.1 (Scn1a), Nav1.6 (Scn8a), Nav1.7 (Scn9a), and Nav1.9 (Scn11a) were enriched in DRG. The expression level of the sodium level-sensing channel Nax (encoded by Scn7a) was 6.5-fold higher in visceral DRG vs. NG neurons. In contrast, high levels of the cardiac channel Nav1.5 mRNA (Scn5a) were detected in vagal but not spinal afferents.
|
| DISCUSSION |
|---|
|
|
|---|
Spinal sensory neurons supplying somatic and visceral receptive fields have a number of features in common but differ in terms of morphology, conduction parameters, and expression of ion channels and receptors (29, 36). Vagal neurons are also heterogeneous in their terminal innervation to cardiac, pulmonary, and GI structures and express a wide range of ion channels and receptors, some of which are specific for different populations (47). Our fiber-tracing experiments provide quantitative information on the extent to which abdominal viscera-projecting neurons contribute to the total pool of sensory neurons in DRG and NG. In T10T13 DRG,
3% of neurons were labeled after ip CTB injection. Previous studies in the rat estimate the proportion of visceral afferents to be <7% of the total (10). As demonstrated before, ip administration of tracer allows enrichment of sensory afferents projecting within the ip cavity (2, 32). Clearly, the extent to which endings in the peritoneal lining are labeled cannot be determined (42). All im CTB-labeled neurons were also labeled by ip injection, and extrapolation of the number of neurons that supply a 5-cm segment of jejunum to the whole bowel gives a close approximation of the total population labeled by ip injection. The same holds true for NG. The importance of this information is twofold. First, it gives confidence that a specific subset of sensory neurons that supply the abdominal viscera can be labeled, harvested, and enriched for without any prior surgical intervention. Such surgical procedures are used to elicit paralytic ileus, which recent studies have shown to be dependent on an inflammatory response within the bowel wall (15, 16). Because inflammation leads to long-term changes in sensory neuron excitability, we felt it was necessary to develop a method of labeling that was relatively noninvasive to establish a baseline of gene expression in nodose and DRG neurons. Second, it gives an estimate of the extent to which expression from these signals is likely to be diluted when whole ganglia expression is assessed. In this respect, the comparison of data obtained from whole DRG and NG with data obtained from laser-captured neurons using the same microarrays is important. Even in whole DRG vs. NG, a large number of genes are differentially expressed (8% of the reliably detected genes). This may be not surprising, given that DRG and NG neurons have different embryonic origins (neural crest and epibranchial placodes, respectively; Refs. 13, 30), innervate different target organs (largely somatic structures vs. visceral organs), and generate different sensory signals (conscious perception vs. reflex and behavior control; Ref. 9). Notable genes that appear in the analysis of whole ganglia reflect these differences. These include the 5-HT3A receptor, the Nav1.8 sodium channel, and the CCKA receptor. Some of these differences in gene expression observed in whole ganglia are preserved when laser-captured neurons are compared (e.g., 5-HT3A and CCKA). Other genes drop out (e.g., Nav1.8), illustrating the extent to which somatic neurons dominate expression in whole ganglia. A third category of genes is not differentially expressed when whole ganglia are compared, but is clearly different in laser-captured material. This reflects the extent to which dilution of the signal can obscure important differences in the specific population of neurons that innervate the abdominal viscera.
Important validation of our study comes from expression profiles of genes such as those encoding 5-HT3A, Nav1.8 and CCKA, for which the presence in visceral afferents has previously been documented (5, 7, 37). In addition, we have shown a good correlation between expression measured by qPCR and microarray technology for 18 genes.
In the further discussion below, we have highlighted genes that play a role in sensory signaling. One example is the capsaicin receptor Trpv1, considered important in nociception (8). Its expression in both DRG and NG confirms earlier immunohistochemical data (25, 35, 44) and may indicate a role for both vagal and spinal afferents in nociception. Trpa1 was also readily detectable in both DRG and NG neurons. This channel is gated by mechanical stimuli as well as by cold temperatures, mustard oil, and cannabinoids (12, 27). As such, Trpa1 could confer mechanosensitivity on the peripheral endings of both spinal and vagal afferents.
Other gene families of particular importance for sensory neurons are those that determine membrane excitability. Sodium and potassium channels are particularly important, especially given recent experimental data showing hypersensitivity resulting from upregulation of sodium and downregulation of potassium channels after inflammation or injury (5, 26, 39). For the purpose of this discussion, we have focused only on voltage-gated sodium channels. The
-subunits determine much of the biophysical and pharmacological properties, and nine
-subunits have been identified for expressed mammalian sodium ion channels (34). In unselected DRG, Nav1.7 and, to a lesser extent, Nav1.1, Nav1.2, and Nav1.6 carry much of the TTX-sensitive current (6). In GI-projecting DRG and NG neurons, the TTX-sensitive current has biophysical properties resembling Nav1.7 (encoded by Scn9a). In this respect, Scn9a is expressed in both DRG and NG. A recent study using the Cre-loxP system to generate nociceptor-specific knockouts for Nav1.7 found increased thresholds to mechanical and thermal pain in the Nav1.7-deficient mice (33). TTX-resistant channels Nav1.8 and Nav1.9 have received a great deal of attention recently because of their preferential localization to small, unmyelinated, nociceptor-like neurons (18). Scn10a (Nav1.8) mRNA expression and immunocytochemical staining for Nav1.8 have been demonstrated in DRG neurons innervating the colon, ileum and stomach, and Nav1.8-like currents are present in most DRGs innervating these organs (41). Only one study has systematically looked for persistent Nav1.9-like currents in GI afferents, showing that only a small proportion (
13%) of colonic DRGs had a persistent, low threshold current (5). Thus it seems likely that Nav1.8 carries the major TTX-resistant current in GI afferents; however, evidence demonstrating either mRNA or immunohistochemically detectable Nav1.8 or Nav1.9 in gut-specific neurons is lacking. In our study we found that Nav1.9 (Scn11a) mRNA is abundant in DRG and, to a lesser extent, also expressed in NG, while Nav1.8 mRNA expression level is comparable in both neuronal populations. The functional significance of Nav1.5 expression in NG neurons needs to be explored, as it has been described as the cardiac sodium channel (17). Of note, mRNA for one of the auxiliary ß-subunits (Scn3b) that alters the kinetics of Nav1.5 is also enriched in NG neurons (19). Detection of Nav1.5 mRNA expression in NG but not in DRG visceral neurons may indicate that this channel contributes to the sodium ion current in visceral vagal afferents (38).
Glutamate is a major transmitter at first-order synapses in the brainstem and spinal cord. There are at least six NMDA receptors, nine non-NMDA receptors, and eight metabotropic glutamate receptors (mGluRs) (11, 45). Both pre- and postsynaptic mGluRs have been implicated in shaping autonomic signal transmission in the nucleus tractus solitarius (21, 24). Similarly, NMDA and non-NMDA receptors influence the release of both glutamate and substance P from spinal nociceptive neurons (3, 31). Therefore, an examination of the expression profile of glutamate receptors in DRG and NG neurons is important for understanding sensory transmission. The abundant expression of the kainate receptor GluR5 (Grik1) in DRG compared with NG neurons strongly supports the view that this is the predominant presynaptic regulator of visceral afferent input to the spinal cord (28). Gria1, encoding the AMPA receptor GluR1, is also more abundant in DRG compared with NG. Differential incorporation of specific AMPA subunits has been shown to modulate pain-related behavior (23). Our expression data suggest that, in visceral afferents, heteromers of GluR1/GluR2 or GluR1/GluR4 might predominate at the spinal level, whereas at the vagal level, GluR3-containing heteromers may be more common. Specific mGluR subtypes present at the synapses also modulate glutaminergic transmission. In general, group I mGluRs (subtypes-1 and -5) increase neuronal excitability, while group II (subtypes-2 and -3) and group III (subtypes-4, -6, -7, and -8) receptors inhibit glutamate release to reduce synaptic transmission (11). Our current findings suggest that the presynaptic mechanisms controlling transmitter release may be different for DRG and NG neurons, with mGluR5 and mGluR8 being dominant in the spinal transmission and none of the mGluRs being enriched in vagal afferents.
In conclusion, the present study provides a genome-wide insight into the molecular signatures that underlie both functional differences and similarities between NG and DRG sensory neurons projecting to the peritoneal cavity. These divergences may underlie the contrasting roles of NG and DRG and contribute to differences in the control of excitability between these visceral sensory neurons. Moreover, these findings will offer novel insight into the mode of action of pharmacological agents modulating visceral sensation and pave the way for additional, more directed studies investigating the role of novel proteins abundantly expressed in visceral afferents.
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: P. J. Peeters, Johnson & Johnson Pharmaceutical Research and Development, a division of Janssen Pharmaceutica NV, Turnhoutseweg 30, B-2340 Beerse, Belgium (e-mail: ppeeter3{at}prdbe.jnj.com).
1 The Supplemental Material for this article (Supplemental Fig. S1) is available online at http://physiolgenomics.physiology.org/cgi/content/full/00169.2005/DC1. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. Nassenstein, K. Kwong, T. Taylor-Clark, M. Kollarik, D. M. MacGlashan, A. Braun, and B. J. Undem Expression and function of the ion channel TRPA1 in vagal afferent nerves innervating mouse lungs J. Physiol., March 15, 2008; 586(6): 1595 - 1604. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Hillsley, J.-H. Lin, A. Stanisz, D. Grundy, J. Aerssens, P. J. Peeters, D. Moechars, B. Coulie, and R. H. Stead Dissecting the role of sodium currents in visceral sensory neurons in a model of chronic hyperexcitability using Nav1.8 and Nav1.9 null mice J. Physiol., October 1, 2006; 576(1): 257 - 267. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |