AJP - Lung Journal of Neurophysiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Lung Cell Mol Physiol 292: L704-L715, 2007. First published November 10, 2006; doi:10.1152/ajplung.00383.2006 Free Article
1040-0605/07 $8.00
This Article
Free upon publication Free Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Tables
Right arrow All Versions of this Article:
292/3/L704    most recent
00383.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Balbir, A.
Right arrow Articles by Shirahata, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Balbir, A.
Right arrow Articles by Shirahata, M.

A search for genes that may confer divergent morphology and function in the carotid body between two strains of mice

Alexander Balbir,1 Hannah Lee,2 Mariko Okumura,1 Shyam Biswal,2 Robert S. Fitzgerald,1 and Machiko Shirahata1

Department of Environmental Health Sciences, Divisions of 1Physiology and 2Toxicology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland

Submitted 27 September 2006 ; accepted in final form 9 November 2006


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The carotid body (CB) is the primary hypoxic chemosensory organ. Its hypoxic response appears to be genetically controlled. We have hypothesized that: 1) genes related to CB function are expressed less in the A/J mice (low responder to hypoxia) compared with DBA/2J mice (high responder to hypoxia); and 2) gene expression levels of morphogenic and trophic factors of the CB are significantly lower in the A/J mice than DBA/2J mice. This study utilizes microarray analysis to test these hypotheses. Three sets of CBs were harvested from both strains. RNA was isolated and used for global gene expression profiling (Affymetrix Mouse 430 v2.0 array). Statistically significant gene expression was determined as a minimum six counts of nine pairwise comparisons, a minimum 1.5-fold change, and P ≤ 0.05. Our results demonstrated that 793 genes were expressed less and that 568 genes were expressed more in the A/J strain vs. the DBA/2J strain. Analysis of individual genes indicates that genes encoding ion channels are differentially expressed between the two strains. Genes related to neurotransmitter metabolism, synaptic vesicles, and the development of neural crest-derived cells are expressed less in the A/J CB vs. the DBA/2J CB. Through pathway analysis, we have constructed a model that shows gene interactions and offers a roadmap to investigate CB development and hypoxic chemosensing/chemotransduction processes. Particularly, Gdnf, Bmp2, Kcnmb2, Tph1, Hif1a, and Arnt2 may contribute to the functional differences in the CB between the two strains. Bmp2, Phox2b, Dlx2, and Msx2 may be important for the morphological differences.

hypoxia; microarray; development; glomus cell


ABNORMALITIES IN THE VENTILATORY response to hypoxia and hypercapnia during sleep are hallmark features in disorders such as congenital central hypoventilation syndrome, sudden infant death syndrome, and obstructive sleep apnea. The role of the carotid body (CB) in the autonomic reflex response to hypoxia is very well defined (35). Dysfunction, as well as the size of the CB, may contribute to the aforementioned disorders (19, 23, 48). Hypoxia increases ventilation, but this response is variable among individuals (28, 106, 109). Genetics may play a critical role in explaining this variability. Indeed, longitudinal studies in the same individuals and studies in twins do demonstrate the role of genetics in the hypoxic ventilatory response (20, 53). Studies utilizing inbred strains of mice have also demonstrated the effect of genetics on the response to hypoxia (102). In particular, two strains of mice, DBA/2J and A/J, were identified as high and low responders to hypoxia, respectively. To further elucidate phenotypic differences between DBA/2J and A/J mice, Rubin et al. (93) examined these two strains in a model of sleep-induced hypoxia. The DBA/2J strain demonstrated a higher sensitivity to hypoxia during sleep, as evidenced by shorter time to arousal, compared with that of the A/J strain. Many factors may contribute to this phenotype of depressed hypoxic sensitivity in the A/J strain. However, a prominent contributor to this phenotype appears to be the CB. We (114) have shown that the volume of the A/J CB is significantly smaller than that of the DBA/2J CB. Furthermore, the A/J CB contains significantly fewer chemosensory cells (glomus cells) compared with its DBA/2J counterpart (114). In addition to these morphological differences, the glomus cell sensitivity to acetylcholine (ACh) is lower in the A/J mice than that in the DBA/2J mice (114). Therefore, both morphological and functional abnormalities in the CB may correlate to decreased hypoxic sensitivity in the A/J strain.

Exposure to environmental stressors plays a critical role in CB structure and function. Studies have indicated that exposing animals to chronic hypoxia results in complex morphological, genetic, and functional changes (13, 59, 85). Some of the changes include enlargement of the CB (46), increased or decreased CB response to acute hypoxia (59), and modified neurotransmitter responses (13, 43, 44). From a developmental standpoint, exposure to hyperoxia at early neonatal periods induces hypoplasia of the CB and decreases systemic response to hypoxia during the adult stage (38). Also, postnatal exposure to hypoxia decreases glomus cell sensitivity to hypoxia (98).

However, no environmental stressors were imposed on the DBA/2J and A/J mice preceding the studies described above. Therefore, we can presume that differences in the CB structure, as well as hypoxic sensitivity of the CB, between these two strains are a direct result of genetics. However, the genes responsible for the phenotypes described above are relatively unknown. We hypothesize that gene expression in the CBs of these two strains of mice are divergent, thus contributing to the disparity in the phenotypes. More specifically, 1) genes related to CB function are expressed less in the A/J mice compared with the DBA/2J mice, and 2) gene expression levels of morphogenic and trophic factors of the CB are significantly lower in the A/J mice compared with the DBA/2J mice. To test these hypotheses, we utilized microarray analysis to examine transcriptional differences between CBs of both strains of mice. This is the first attempt in the field to demonstrate differences in basal gene expression in the CBs between the hypoxic sensitive mice (DBA/2J) and the hypoxic insensitive mice (A/J). These genes may regulate components related to hypoxic sensitivity in the CBs of these mice.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Animals. Male DBA/2J and A/J strains were initially purchased from Jackson Laboratories (Bar Harbor, ME), later bred in the animal facility of the Johns Hopkins Bloomberg School of Public Health (Baltimore, MD), and housed with proper temperature (~22°C) and regulated light cycle (12:12-h light-dark cycle). Food (Agway Pro-Lab RMH 1000) and water were available ad libitum before tissue harvesting. All animal experiments were conducted in accordance with the Guiding Principles of the American Physiological Society for the Care and Use of Animals. Protocols for animal experimentation were approved by the Johns Hopkins University Animal Care and Use Committee.

Tissue preparation. Due to the low amount of RNA per CB, it was necessary to harvest large quantities of CBs to perform microarray analysis. For both strains, three sets of CBs were established (Table 1) to perform adequate statistical analyses (see Microarray data analysis; 3 x 3 comparison). Mice were anesthetized with urethane by peritoneal injection (1.25 g/kg). Bilateral carotid sinuses were exposed via a midline incision of the neck, thus permitting visualization of CB region. The heart was excised to prevent excessive bleeding into the CB region. Immediately following heart excision, carotid bifurcations were harvested and deposited into chilled phosphate-buffered saline. CBs were cleaned of extraneous tissue, deposited in ice-cold TRIzol reagent (Invitrogen), and frozen at –80°C until use. CBs were pooled into appropriate sets for subsequent RNA isolation.


View this table:
[in this window]
[in a new window]

 
Table 1. Mice used for microarray analysis

 
RNA isolation for microarray analysis. RNA was isolated utilizing the TRIzol method (Invitrogen). Briefly, each set of CBs was homogenized in 0.8 ml of TRIzol reagent using a pestle homogenizer. Insoluble tissue was removed via centrifugation of the sample tube at 12,000 rpm for 10 min at 4°C. Clear homogenate solution was then transferred to a fresh microcentrifuge tube. The homogenate solution was allowed to sit at room temperature (~23°C) for 5 min to further dissociate nucleoprotein units. Chloroform (0.16 ml) was added to the homogenate, agitated, and allowed to sit at room temperature for 2–3 min. The sample was again centrifuged at 12,000 rpm for 15 min at 4°C. The RNA-containing, aqueous phase was placed in a fresh microcentrifuge tube. Two microliters of glycogen (5 µg/µl) was added to facilitate RNA precipitation, followed by an addition of isopropanol (0.4 ml). The sample was incubated at room temperature for 10 min and centrifuged at 12,000 rpm for 10 min at 4°C. The RNA pellet was identified, and the overlaying supernatant was discarded. The RNA pellet was washed with 0.8 ml of 75% ethanol, vortexed for mixing, and centrifuged at 7,600 rpm for 10 min at 4°C. The ethanol was discarded, and the RNA pellet was air-dried for 5–10 min. RNA was eluted in 10 µl of TE buffer and incubated for 10 min at 55°C.

Amplification, labeling, fragmentation, and hybridization for microarray. Transcriptional profiling via microarray has been previously described (89, 103). Briefly, total RNA for each set was analyzed for purity and concentration (Agilent 2100 Bioanalyzer) (Table 1). Microarray experiments were performed with the Affymetrix Mouse Genome 430 v2.0 GeneChip Array. Hybridization method was adapted from previous work (18, 89, 103). Briefly, the Two-Cycle Target Labeling protocol for GeneChip Expression Analysis (Affymetrix) was used. Beginning with 100 ng of total RNA for the first cycle of amplification, the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen) was used to convert RNA into first cycle, first and second strand cDNA. After purification with the GeneChip Sample Cleanup Module (Affymetrix), first cycle cDNA was converted to cRNA for amplification using the MEGAscript T7 Kit (Ambion). The first cycle cRNA was cleaned with the RNeasy Mini Kit (Qiagen) and converted to second cycle, first and second strand cDNA in a similar manner to first cycle cDNA except that a T7-oligo(dT) promoter primer (Affymetrix) was incorporated. The second cycle cDNA was cleaned in a similar manner to the first cycle cDNA and then converted to biotin-labeled cRNA with the BioArray High Yield RNA Transcript Labeling Kit (ENZO). Labeled cRNA was cleaned with the GeneChip Sample Cleanup Module (Affymetrix). Fragmentation was achieved by addition of 5x Fragmentation Buffer (Affymetrix). External control was Control Oligo B2 (3 nM), Eukaryotic Hybridization control (Affymetrix). Hybridization protocol consisted of 200 µl of hybridization cocktail containing 10 µg of fragmented cRNA placed on a microarray chip. Hybridization occurred for 18 h in the GeneChip Hybridization Oven 640 (Affymetrix). After hybridization, the chip was washed and stained using two streptavidin stains and one antibody stain on the GeneChip Fluidics Station 400 (Affymetrix) for 1.5 h. Arrays were scanned with a GCS300 Scanner (Affymetrix) for images and intensities.

Microarray data analysis. Data extraction was processed by Affy GeneChip Operating Software v1.3 (Affymetrix). Global scaling was equal to 500. Normalization factor was equal to 1. Analysis of differentially expressed transcripts was accomplished using a pairwise comparative analysis (Affy Data Mining Tool v3.1; Affymetrix). Analysis compares the differences in values of perfect match to mismatch of each probe pair in the A/J strain array to its matching probe pair on the DBA/2J strain array. P values were determined by the Wilcoxon signed rank test and denoted as increase, decrease, or no change. Signal log ratio, obtained from Data Mining Tool, estimates the magnitude and direction of change in a transcript when two arrays are compared (A/J vs. DBA/2J). Signal log ratio output was converted to fold change using the following equation recommended by Affymetrix:

Formula 1(1)

In our current study, we examined nine pairwise comparisons (A/J, n = 3 sets vs. DBA/2J, n = 3 sets). Significant, differentially expressed genes were selected based on the following inclusion criteria: 1) minimum six counts out of nine confirming pairwise comparisons; 2) >1.5-fold change average; 3) Mann-Whitney pairwise comparison test: P ≤ 0.05 (test performed to rank the results by concordance as a calculation of significance of each identified change in gene expression) (103). The expressed sequence tags of the selected genes were deciphered for their annotation in the NetAffx Analysis Center of the Affymetrix site (http://www.affymetrix.com). Also, gene ontology (GO), and subsequent GO mining information, was obtained at the NetAffx Analysis Center. The microarray data was submitted to the GEO database (accession no. GSE5313).

Real-time RT-PCR analysis. RNA was isolated from DBA/2J and A/J CBs utilizing the same TRIzol method as described above. However, for this analysis, a pair of CBs was individually used for RNA isolation, and RNA was treated with DNase I (Invitrogen). Briefly, 3 µl of 10x DNase I Reaction Buffer and 3 µl of DNase I Amp Grade was added to the RNA solution for a total volume of 26 µl. The solution was incubated at room temperature for 15 min followed by enzyme inactivation for 10 min at 65°C. cDNA synthesis was accomplished with the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen). RNA-primer mixture was composed of 2 µl of primer (50 ng/µl random hexamers) and 2 µl of 10 mM deoxynucleotide triphosphate (dNTP) mix plus 16 µl of RNA sample for a total solution volume of 20 µl. The RNA-primer mixture was incubated for 5 min at 65°C followed by a 1-min ice bath treatment. The cDNA synthesis mix solution was composed of 4 µl of 10x RT buffer, 8 µl of 25 mM MgCl2, 4 µl 0.1 M DTT, 2 µl of RNaseOUT (40 U/µl), and 2 µl of SuperScript III RT (200 U/µl) for a total volume of 20 µl. The cDNA synthesis mix solution was added to 20 µl of the RNA-primer mixture for a total volume of 40 µl and incubated for 10 min at 25°C followed by 50 min at 50°C. The reaction was terminated at 85°C for 5 min. The reaction mixture was chilled in an ice bath for 1 min, followed by the addition of 2 µl of RNase H, and incubated for 20 min at 37°C. The final volume of cDNA available for real-time PCR assay was 42 µl.

TaqMan gene expression assays (Applied Biosystems), consisting of primer/probe solutions, were used to detect transcriptional levels of endogenous control and experimental genes. The endogenous control was mouse beta-actin (Mouse ACTB). Experimental genes included glial cell line-derived neurotrophic factor (Gdnf; Mm00599849_m1), bone morphogenetic protein-2 (Bmp2; Mm01340178_m1), potassium large conductance calcium-activated channel, subfamily M, beta-member-2 (Kcnmb2; Mm00511481_m1), tyrosine hydroxylase (Th; Mm00447546_m1), choline acetyltransferase (Chat; Mn01221876_m1), tryptophan hydroxylase-1 (Tph1; Mm00493794_m1), hypoxia-inducible factor-1, {alpha}-subunit (Hif1a; Mm00468869_m1), aryl hydrocarbon receptor nuclear translocator-2 (Arnt2; Mm00476009_m1), paired-like homeobox-2b (Phox2b; Mm00435872_m1), distal-less homeobox-2 (Dlx2; Mm00438427_m1), and homeobox msh-like-2 (Msx2; Mm00442992_m1).

PCR amplification of cDNA was performed with a total volume of 25 µl, comprising 12.5 µl of TaqMan Universal PCR Master Mix (Applied Biosystems), 1.25 µl of TaqMan primer/probe assay, 10.25 µl of diethyl pyrocarbonate (DEPC)-treated water (Quality Biological), and 1 µl of template cDNA. Real-time PCR amplification was performed in an ABI PRISM 7000 Sequence Detection System (Applied Biosystems) according to the following protocol: 2 min at 50°C, 10 min at 95°C, and 50 cycles of 95°C for 15 s followed by 60°C for 1 min. Individual cDNA samples were run in triplicates for each gene of interest, including endogenous control. After completion of real-time amplification, threshold designation was determined on the amplification plot using ABI PRISM 7000 Sequence Detection System Software (Applied Biosystems). Reaction efficiencies were similar for all real-time assays due in part to adherence to manufacturers' suggestions (8). Due to the small quantity of RNA per pair of CBs (~6 ng/µl), we utilized a "semiquantitative" approach for analysis. The threshold for all genes was established at the ~1/3 point of the linear phase of replication. Mean cycle threshold (Ct) was calculated for each gene transcript. Then the expression value of each gene from a pair of CBs was expressed as {Delta}Ct, where {Delta}Ct = (mean Ctgene) – (mean Ctbeta-actin). A cycle difference of 1 corresponds to a twofold change in gene expression. The data from several animals were expressed as mean {Delta}Ct ± SE, and the Mann-Whitney test was performed between the A/J and the DBA/2J. In addition, we used the Relative Expression Software Tool v1.9.12 (REST 2005, Corbett Life Science; http://rest-2005.gene-quantification.info) to compare Ct values for genes of interest between two samples (A/J vs. DBA/2J). This semiquantitative analysis for real-time RT-PCR is based on a mathematical randomization model that has been previously described (83).


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Animals. The age of mice (29.5 ± 0.18 days, mean ± SE, n = 54 in DBA/2J; 29.3 ± 0.06 days, mean ± SE, n = 60 in A/J) and body weight (13.7 ± 0.28 g, mean ± SE, in DBA/2J; 13.6 ± 0.23 g, mean ± SE, in A/J) were not significantly different (unpaired t-test, P > 0.05 for both variables). The ages and body weights among the three experimental sets were not statistically different (one-way ANOVA, P > 0.05; Table 1).

Microarray profile. Microarray quality and hybridization information for each experimental set can be found in Supplement 1. Normalized and summarized microarray data for the comparative CB analysis can be found in Supplement 2. Data supplements are available online at the AJP-Lung Cellular and Molecular Physiology web site.

Gene expression. A total of 10,261 probe sets were expressed less and 7,285 probe sets were expressed more in the A/J CB vs. the DBA/2J CB. Based on inclusion criteria for significance, 738 genes were identified as expressed less in the A/J strain vs. the DBA/2J strain, and 554 genes were expressed more in the A/J strain vs. the DBA/2J strain. Gene expression profiles, increased and decreased, were subsequently organized into their respective GO branch: Cellular Component, Molecular Function, and Biological Process. Within each GO branch, genes were further sorted according to functional annotation terms (Fig. 1). A large number of genes involved in signal transducer activity, catalytic activity, cell, and development are expressed less in the A/J CB. Genes associated with binding, catalytic activity, envelope, and cellular process are expressed more in the A/J CB. Furthermore, we examined individual genes and their relevance to the CB. Several genes from which products are confirmed or suggested to be involved in CB function were identified. Those include genes involved in ion channels, the metabolism of neurotransmitters, neurotransmitter receptors, transporters, and putative chemosensors. Some of these genes were expressed less in the CB of the A/J mice than that of the DBA/2J mice (Table 2). Additionally, we examined other genes that may also be involved in ion movement, the metabolism of neurotransmitters, neurotransmitter receptors, synaptic vesicles, second messengers, and putative chemosensors. Table 3 lists these genes that were differentially expressed in the CBs of the two strains. These genes may play a role in CB function. The final genes identified were those which are considered to be involved in morphogenesis, differentiation, and degeneration in neural crest-derived tissues. Many genes were differentially expressed between the two strains of mice (Table 4). A brief description of each gene in Tables 24 is in Supplement 3. To understand the interactions among differentially expressed genes, pathway analysis (Ingenuity Pathway Analysis, Ingenuity Systems) was performed. Based on this analysis and our independent literature searches, we have constructed a simplified model that explains possible pathways contributing to the morphological and functional phenotypic differences in the CB between the A/J and DBA/2J mice (Fig. 2).


Figure 1
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 1. Gene Ontology (GO) comparison of the differentially expressed genes in the carotid body (CB). A: genes expressed less in the A/J CB vs. the DBA/2J CB. B: genes expressed more in the A/J CB vs. the DBA/2J CB. GO branches are Cellular Component, Molecular Function, and Biological Process. GO branches are subdivided into functional annotation terms. The GO Mining Tool can be found at the NetAffx Analysis Center (http://www.affymetrix.com). A detailed description of the terms and features is available in the GO Mining Tool Manual (https://www.affymetrix.com/support/technical/manual/go_manual.affx). The parameters established for analysis of our microarray profiles were as follows: 1) Ontology branch: Cellular Component, Molecular Function, and Biological Process; 2) GeneChip Array: 430_2 (i.e., Affymetrix Mouse Genome 430 2.0); 3) Probe Set Count Threshold: 1; 4) Probe Set Percentage Threshold: 0.05; 5) Probe Set Probability Threshold: 0.05 (i.e., P value threshold for significance). The P values calculated reflect significance of paired distribution ({chi}2 test). Functional annotation terms with P < 0.05 are presented.

 

View this table:
[in this window]
[in a new window]

 
Table 2. Differentially expressed genes that pertain to CB function

 

View this table:
[in this window]
[in a new window]

 
Table 3. Differentially expressed genes that may be related to CB function

 

View this table:
[in this window]
[in a new window]

 
Table 4. Novel, differentially expressed genes that may be related to CB morphology

 

Figure 2
View larger version (48K):
[in this window]
[in a new window]

 
Fig. 2. Interactions of genes. Most of the gene interactions were identified and described by Ingenuity Pathway Analysis (Ingenuity Systems). Gene interaction references can be found in Supplement 4. Genes: glial cell-derived neurotrophic factor (Gdnf); bone morphogenetic protein-2 (Bmp2); brain-derived neurotrophic factor (Bdnf); potassium large conductance calcium-activated channel, subfamily M, beta-member-2 (Kcnmb2); Gdnf signaling receptor tyrosine kinase c-RET (Ret); Gdnf family receptor-{alpha}1 (Gfra1); tyrosine hydroxylase (Th); nitric oxide synthase-1 (neuronal NO synthase; Nos1); choline acetyltransferase (Chat); Tachykinin-1 (Tac1); tryptophan hydroxylase-1 (tryptophan 5-monooxygenase; Tph1); single-minded homolog-1 (Sim1); hypoxia-inducible factor-1, {alpha}-subunit (Hif1a); aryl-hydrocarbon receptor nuclear translocator-2 (Arnt2); paired-like homeobox-2a (Phox2a); achaete-scute complex-like-1 (Ascl1); paired-like homeobox-2b (Phox2b); distal-less homeobox-2 (Dlx2); and msh homeobox homolog-2 (Msx2).

 
Real-time RT-PCR analysis. Real-time RT-PCR was performed to confirm the expressional differences of some genes in Fig. 2. The results are presented in Table 5. Statistical analysis of {Delta}Ct mostly corroborates observed differences in the microarray expression trends (i.e., increase or decrease). Also, analysis with the REST 2005 tool further corroborates differences in expression trends.


View this table:
[in this window]
[in a new window]

 
Table 5. Real-time RT-PCR analysis

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This study has demonstrated differences in steady-state gene expression in the CB between the DBA/2J and A/J mice. Previous studies, including ours, demonstrated that the DBA/2J mice have a larger, hypoxic-sensitive CB compared with the A/J mice (76, 102, 114). These differences are likely genetically regulated. Although genetic differences may also exist in the CB afferent system, the current study has focused on the CB. We have hypothesized that 1) genes related to CB function are expressed less in the A/J mice compared with the DBA/2J mice, and 2) gene expression levels of morphogenic and trophic factors of the CB are significantly lower in the A/J mice than the DBA/2J mice. Our microarray analysis, together with real-time RT-PCR, identified many genes that fit these categories and were differentially expressed between the CBs of the two strains. We believe that this work will lead to future studies to examine the role of each gene in CB function and morphology.

We demonstrated differences in multiple genes expressed in the CB between the A/J and DBA/2J strains at 4 wk of age (Fig. 1). GO classification and functional annotation terms are broad designations for genes. Many genes have multiple functions; therefore, significant redundancy is incorporated into this manner of gene expression sorting. In addition, the CB is heterogeneous in its cellular composition. It contains not only chemosensory glomus cells, but also supporting sheath cells, in addition to blood vessels, blood cells, nerve fibers, Schwann cells, and fibroblasts. Therefore, Fig. 1 may not represent the gene expression differences in components related to oxygen sensing. Thus we examined all significant, differentially expressed genes and categorized them based on our knowledge regarding the CB function and morphology (Tables 24). The adult A/J mice exhibits a smaller CB with fewer chemosensory cells (114), a decreased ventilatory (102) and cardiovascular response to hypoxia (16), as well as a decreased arousal response to hypoxia during sleep (93) compared with the DBA/2J mice. Therefore, one can suspect that transcriptional differences between CBs of the two strains may account for such phenotypic expression. We have particularly focused on the genes related to current theories of hypoxic chemotransduction to assign significance to the many differentially expressed genes in our study. Due to the large number of genes discovered, we have limited our discussion mainly to the genes that may well be involved in chemoreception/chemotransduction processes (Table 2), the genes involved in the development of neural crest-derived cells (Table 4), and the genes that are closely linked to each other (Fig. 2).

The presence of oxygen-sensitive potassium channels has been established in the CB (41). These include Ca2+-activated K+ channels (also called BK or maxi-K+ channels) (80) and background TASK-like K+ channels in the rat (111), Kv4.1 and Kv4.3 (Kv4.3–1 variant) in the rabbit (94), and Kv3 in the C57bl/6J mouse (82). Furthermore, oxygen-sensitive K+ channels in neuroepithelial bodies are Kv3.3 (Kcnc3) channels (107) and Kv1.2 (Kcna2) channels in PC12 cells (21). Since hypoxic sensitivity of the CB is compromised in the A/J mice, the expression of these channel subunits might well have been less in these mice. In our comparative gene expression data, however, no significant transcriptional differences are seen in Kv channel subunits or TASK channels. However, a specific beta-subunit of BK channels (beta2: Kcnmb2) was expressed significantly less in the A/J CB (Table 2). The result is consistent with our previous report (76). Inhibition of K+ channels by hypoxia is believed to depolarize glomus cells, leading to activation of voltage-dependent Ca2+ channels. In the case of rabbit glomus cells, contribution of voltage-dependent Na+ channels in the depolarizing processes has been shown (67). These channels appear to be expressed less in the CB of the A/J mice compared with that of the DBA/2J mice (Table 2). Decreased expression of these channels may contribute to the low hypoxic sensitivity of the A/J CB as well.

The response to changes in oxygen tension is mediated by neurotransmitters released from glomus cells. The CB contains numerous neurotransmitters and neuromodulators. It is believed that ACh (37) and ATP (118) are two primary, excitatory ligands in the chemotransduction of hypoxia by the CB. Substance P (84) and 5-HT (serotonin) (117) are positive modulators, whereas dopamine has complicated modulatory effects (49). Nitric oxide (NO) (84, 104) and GABA (34) appear to be negative modulators of chemoreceptor function. In the context of our present study, the following genes do not appear to be differentially expressed between the A/J and DBA/2J mice: 1) genes involved in ACh synthesis and in storage, Chat, and vesicular ACh transporter (VAChT); 2) tachykinin-1 (Tac1), which encodes for a precursor of both substance P and neurokinin A (56); 3) NO synthase-1 (neuronal NO synthase). On the other hand, the expression of Tph1, the rate-limiting enzyme in 5-HT synthesis, is significantly lower in the CB of the A/J mice than in that of the DBA/2J mice (Tables 3 and 5). In regard to dopamine, Th, a key enzyme of catecholamine synthesis, appears to be expressed less in the AJ CB compared with the DBA/2J CB (Table 2; 1.82-fold difference). However, analysis of real-time RT-PCR shows an insignificant difference in the expression of Th (Table 5). Furthermore, the expression of aromatic L-amino acid decarboxylase, an enzyme that converts dihydroxyphenylamine to dopamine, was not different between the two strains. Our data also demonstrate that GABA-A receptor subunit beta3 was expressed significantly less in the A/J CB compared with the DBA/2J CB (Table 3). In addition, several transcripts involved in neurotransmitter packaging (Slc6a4) and neurotransmitter vesicle trafficking (Snap25, Stx6, Syt1, and Syt7) were expressed less in the CB of the A/J mice compared with the CB of the DBA/2J mice (Tables 2 and 3). Thus the variable expression levels of genes associated with neurotransmitters and neuromodulators in the CB indicate the complexity associated with different CB chemotransduction between these two strains of mice.

The differences described above may not directly address the issue of the actual hypoxic sensor of the CB. Current theories suggest sensor mechanisms such as NADPH-oxidase (4, 45) and AMP-kinase (AMPK) (30, 31). NADPH-oxidase is thought to couple changes in oxygen tension to K+ channel activity via reactive oxygen species production (45). AMPK may function to link mitochondrial oxidative phosphorylation to calcium-mediated signaling in the CB (30). Another current theory suggests that heme oxygenase-2, with NADPH-cytochrome P-450 reductase, functions to alter BK channel dynamics under different oxygen tensions (112). However, we found no significant differences in the expression of these genes between the A/J CB and the DBA/2J CB. These findings suggest that other genes may regulate the differences in hypoxic sensitivity of the CB between the two mice. For example, Hif1a is expressed less in the A/J CB vs. the DBA/2J CB (Tables 2 and 5). Although precise mechanisms are not well understood, this gene is a candidate in the early steps of hypoxic chemotransduction or perhaps functions in the ability of the CB to adapt to long-term hypoxia or long-term intermittent hypoxia (54, 92).

An important step in the hypoxic chemotransduction paradigm is the integration of the CB output to the central nervous system via chemosensory afferent neurons of the petrosal ganglion. Gdnf and brain-derived neurotrophic factor (Bdnf) have been shown to play critical roles in the survival of chemosensory afferent neurons (29) as well as in the development of medullary neurons responsible for respiratory control (47). Our data showed no difference in Bdnf expression between the A/J and DBA/2J strains. On the other hand, Gdnf is expressed significantly less in the A/J CB compared with the DBA/2J CB (Tables 4 and 5). Thus the low expression of Gdnf in the A/J CB may result in the underdevelopment of chemosensory afferent neurons and may contribute to low hypoxic sensitivity of this strain. A question is whether Gdnf plays a role in the development of the CB. Leitner et al.(62) showed that the addition of Gdnf did not influence the survival of rat glomus cells in culture. However, their study is not conclusive regarding the role of Gdnf in glomus cell survival because Gdnf is naturally present in the rat CB. Their study confirmed that an excess amount of Gdnf did not influence rat glomus cell survival in culture. A/J mice may serve as a good model for studying the effect of Gdnf on respiratory control because of its low expression of Gdnf in the CB.

Previous studies have shown that the transcription factor Ascl1 (Mash1) (50) and Phox2b (79) function as regulators in the development of the autonomic nervous system including the CB. Ascl1 expression in the CB was similar between the A/J and DBA/2J mice, but Phox2b was expressed less in the A/J mice compared with the DBA/2J mice (Tables 2 and 5). In the Phox2b mutant, CB forms embryonically to some extent and then degenerates (24). In addition, Msx2, Dlx1, and Dlx2 were expressed less in the A/J CB vs. the DBA/2J CB. These genes function to control neural crest cell development (116). With respect to development of the CB, we found that the CB of the 1-day-old A/J mouse appears similar in volume and glomus cell quantity compared with its DBA/2J counterpart (10). However, the CB of A/J mice does not appear to grow through adulthood compared with the DBA/2J mice (10). Perhaps the genes mentioned here work together for maintaining morphology and function of the CB.

Thus it is important to know how the differentially expressed genes in the CB of the two strains of mice interact with each other. Figure 2 depicts an intriguing interrelationship among the genes described above. For example, Gdnf can be regulated by the dimerization of two basic helix-loop-helix/PAS transcription factors, Arnt2 (also called Hif2beta) and single-minded homolog-1 (Sim1) (65). Decreased expression of Arnt2 may lead to less expression of Gdnf in the A/J CB (Tables 3 and 5), although Sim1 expression is not significantly different between both strains. Gdnf has been shown to regulate BK channel expression in chick lumbar spinal motoneurons (69). We found that gene expression for BK channel beta2-subunit (Kcnmb2) is lower in the CB of A/J mice, suggesting a regulatory role of Gdnf in BK channel expression in the CB as well. Understanding the relationship among Arnt2, Gdnf, and Kcnmb2 in the CB may expand our knowledge of the hypoxic response of the CB.

Another intriguing gene found in our comparative CB gene expression analysis is Bmp2. Bmp2 was expressed less in the A/J CB compared with the DBA/2J CB. Bmp2 induces expression of paired-like (aristaless) homeobox-2a (Phox2a) (66), Gdnf signaling receptor tyrosine kinase c-RET (66), and Msx2 (96). Also, Bmp2 could regulate the neural differentiation into either cholinergic or noradrenergic neurons (110). The binding partner of Msx2, Dlx2, is involved in neuron differentiation into specific phenotypes, including Th positive cells (7). These genes could contribute to the phenotypic differences in morphology of the CB as well as neurotransmitter expression between the two strains of mice. A particular interest is an upregulatory effect of Bmp2 on Tph1 (33). Because Tph1 is the only significant, differentially expressed gene among genes that control known neurotransmitters in the CB, further studies would provide important information regarding Bmp2's respective role in CB development and serotonin production.

In this study, we used the microarray technique to profile differential gene expression in the CBs of two strains of mice that were not exposed to any environmental stressors. Ganfornina et al. (39) examined the CB gene expression profile in mice (female, C57BL/6J01aHsd) that were exposed to physiological hypoxia. They identified groups of genes that were either switched on or switched off under physiological hypoxia. We compared our results to theirs under the premise that the A/J mice are low responders to hypoxia and therefore may not express genes identified as being switched on during hypoxia in the above study. The switched on genes are: serine (or cysteine) proteinase inhibitor, clade A, member 3 N (Serpina3n); methyl-CpG binding domain protein-1 (Mbd1); a disintegrin and metalloproteinase domain 15 (metargidin; Adam15); cytochrome P-450, family 2, subfamily b, polypeptide 20 (Cyp2b20); heat shock protein 12B (Hspa12b); and protease, Ser19 (neuropsin; Prss19). We identified similar genes from our study that were expressed less in the CB of the A/J vs. the DBA/2J: Serpina1b, Serpina3c, Adam1a, Cyp1b1, Hspa2, and Prss12. However, Mbd1 was the only gene from our data set that matched with their data and it was expressed less in the A/J CB than in the DBA/2J CB. Cyp2j9, a cytochrome P-450 family member, was found to be expressed more in the A/J CB than in the DBA/2J CB. It may be premature to state that these genes are involved in the hypoxic response of the CB due to the strain and gender differences between these two studies. Future studies including both sexes, A/J and DBA/2J mice exposed to hypoxia, and other strains would shed light on the roles of these genes in the CB hypoxic chemotransduction process.

This study offers a comprehensive analysis of the differential gene expression in the CBs of the A/J and DBA/2J strains of mice. We consider this study a useful roadmap to investigate CB development and hypoxic chemosensing/chemotransduction processes because these two strains of mice show clear differences in CB morphology and hypoxic sensitivity. However, some cautions must be applied. For this study, we used 4-wk-old animals. In our previous study (10), we found that glomus cell density within the CB was similar in young animals of both strains. However, the glomus cell density begins to diverge at 4 wk of age (10). Thus genes differentially expressed at this age might be related to some morphological and functional differences of the CB between the two strains. In addition, other genes may well be responsible for the differences between the two strains at different ages, including sex-dependent genes. Some technical issues also must be considered for interpretation of our data. For example, the CB has heterogeneous cellular populations. Therefore, careful verification is required to determine whether a certain gene is specific for glomus cells. Also, due to the small size of the CB, possible contamination of other tissues cannot be totally ruled out. In the current study, the CB was dissected under a dissecting microscope (x30 magnification) with microsurgical instruments. It was technically impossible to identify the microstructures of the CB and eliminate all extraneous tissues. This is an inherent problem with use of the whole CB. Therefore, our gene expression data and the subsequent interpretation must be evaluated carefully for significance. Particularly, the relatively small fold changes in expression would have to be examined further to evaluate their significance in this model of gene expression. In future studies, additions of histological techniques (e.g., in situ hybridization and immunohistochemistry) and possibly gene expression profiling on the single-cell level would help us to determine the possible location of an mRNA transcript, its protein product, and its relative expression in the CBs between these two strains. Regarding the confirmation of microarray analysis, we interpreted real-time transcript expression relative to endogenous control (beta-actin) because initial RNA of the CBs was not quantified due to small sample size. We believe that beta-actin is justified as an endogenous control in this study for the following reasons. First, in our microarray analysis, in which the amount of RNA used was controlled, beta-actin expression in the CB did not differ between the two strains. Second, beta-actin expression per a pair of CBs in real-time RT-PCR analysis was similar between both strains (mean Ct ± SE: 27.04 ± 0.21 in DBA/2J, n = 28; and 27.07 ± 0.23 in A/J, n = 29). Third, beta-actin expression in the brain (used as a representative neural tissue) was indistinguishable between the two strains when normalized with total RNA quantity (A. Balbir, unpublished observation). Fourth, the utilization of the REST 2005 software added additional statistical power by analyzing the relationship between Ct values of the target genes and the beta-actin, endogenous control that was used for this study.

In summary, the transcriptional differences uncovered in this investigation generally support our hypotheses. The results implicate many genes involved in the development of the CB and associated hypoxic chemotransduction. All of these genes could function in concert, or independently, during the early stages of life to establish functional CB chemotransduction processes. However, further detailed studies will be required to determine the relationship between the genes depicted in Fig. 2. Also, future studies will need to address how these genes actually influence the function and development of the CB and its afferent system.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This work was supported by National Institutes of Health Grants HL-72293, P30-ES-038819, and AHA-0255358N. A. Balbir, R. S. Fitzgerald, and S. Biswal are supported by T32-HL-07534, HL-50712, and HL-081205, respectively.


    FOOTNOTES
 

Address for reprint requests and other correspondence: M. Shirahata, Division of Physiology, Dept. of Environmental Health Sciences, Johns Hopkins Bloomberg School of Public Health, E7610, 615 N. Wolfe St., Baltimore, MD 21205 (e-mail: mshiraha{at}jhsph.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Abeliovich A, Schmitz Y, Farinas I, Choi-Lundberg D, Ho WH, Castillo PE, Shinsky N, Verdugo JMG, Armanini M, Ryan A. Mice lacking {alpha}-synuclein display functional deficits in the nigrostriatal dopamine system. Neuron 25: 239–252, 2000.[CrossRef][ISI][Medline]
  2. Abudara V, Eyzaguirre C. Electrical coupling between cultured glomus cells of the rat carotid body: observations with current and voltage clamping. Brain Res 664: 257–265, 1994.[CrossRef][ISI][Medline]
  3. Abudara V, Garces G, Saez JC. Cells of the carotid body express connexin43 which is up-regulated by cAMP. Brain Res 849: 25–33, 1999.[CrossRef][ISI][Medline]
  4. Acker H, Dufau E, Huber J, Sylvester D. Indications to an NADPH oxidase as a possible pO2 sensor in the rat carotid body. FEBS Lett 256: 75–78, 1989.[CrossRef][ISI][Medline]
  5. Akakura N, Kobayashi M, Horiuchi I, Suzuki A, Wang J, Chen J, Niizeki H, Kawamura Ki Hosokawa M, Asaka M. Constitutive expression of hypoxia-inducible factor-1alpha renders pancreatic cancer cells resistant to apoptosis induced by hypoxia and nutrient deprivation. Cancer Res 61: 6548–6554, 2001.[Abstract/Free Full Text]
  6. Allen AM. Angiotensin AT1 receptor-mediated excitation of rat carotid body chemoreceptor afferent activity. J Physiol 510: 773–781, 1998.[Abstract/Free Full Text]
  7. Andrews GL, Yun K, Rubenstein JL, Mastick GS. Dlx transcription factors regulate differentiation of dopaminergic neurons of the ventral thalamus. Mol Cell Neurosci 23: 107–120, 2003.[CrossRef][ISI][Medline]
  8. Applied Biosystems. Amplification efficiency of TaqMan assays- on-demand gene expression oroducts [Online]. 2004. http://docs.appliedbiosystems.com/pebiodocs/00113186.pdf.
  9. Assadi AH, Zhang G, Beffert U, McNeil RS, Renfro AL, Niu S, Quattrocchi CC, Antalffy BA, Sheldon M, Armstrong DD, Wynshaw-Boris A, Herz J, D'Arcangelo G, Clark GD. Interaction of reelin signaling and Lis1 in brain development. Nat Genet 35: 270–276, 2003.[CrossRef][ISI][Medline]
  10. Balbir A, Okumura M, Schofield B, Coram J, Tankersley C, Fitzgerald RS, O'Donnell CP, Shirahata M. Genetic regulation of chemoreceptor development in DBA/2J and A/J strains of mice. Adv Exp Med Biol 580: 99–104, 2006.[ISI][Medline]
  11. Balbir A, Otsubo T, Okumura M, Fujii K, Shirahata M. Calcium-activated large conductance K (BK) channels in glomus cells. J Neurosci: Neuroscience 2005 Conference Abstracts/Annual Meeting Publications (program no. 846.11), 2005.
  12. Belmadani A, Tran PB, Ren D, Assimacopoulos S, Grove EA, Miller RJ. The chemokine stromal cell-derived factor-1 regulates the migration of sensory neuron progenitors. J Neurosci 25: 3995–4003, 2005.[Abstract/Free Full Text]
  13. Bisgard GE. Carotid body mechanisms in acclimatization to hypoxia. Respir Physiol 121: 237–246, 2000.[CrossRef][ISI][Medline]
  14. Buckler KJ, Vaughan-Jones RD, Peers C, Nye PC. Intracellular pH and its regulation in isolated type I carotid body cells of the neonatal rat. J Physiol 436: 107–129, 1991.[Abstract/Free Full Text]
  15. Buniel MC, Schilling WP, Kunze DL. Distribution of transient receptor potential channels in the rat carotid chemosensory pathway. J Comp Neurol 464: 404–413, 2003.[CrossRef][ISI][Medline]
  16. Campen MJ, Tagaito Y, Li J, Balbir A, Tankersley CG, Smith P, Schwartz A, O'Donnell CP. Phenotypic variation in cardiovascular responses to acute hypoxic and hypercapnic exposure in mice. Physiol Genomics 20: 15–20, 2004.[Abstract/Free Full Text]
  17. Ceradini DJ, Kulkarni AR, Callaghan MJ, Tepper OM, Bastidas N, Kleinman ME, Capla JM, Galiano RD, Levine JP, Gurtner GC. Progenitor cell trafficking is regulated by hypoxic gradients through HIF-1 induction of SDF-1. Nat Med 10: 858–864, 2004.[CrossRef][ISI][Medline]
  18. Chen H, Rubin E, Zhang H, Chung S, Jie CC, Garrett E, Biswal S, Sukumar S. Identification of transcriptional targets of HOXA5. J Biol Chem 280: 19373–19380, 2005.[Abstract/Free Full Text]
  19. Cole S, Lindenberg LB, Galioto FM Jr, Howe PE, DeGraff AC Jr, Davis JM, Lubka R, Gross EM. Ultrastructural abnormalities of the carotid body in sudden infant death syndrome. Pediatrics 63: 13–17, 1979.[Abstract/Free Full Text]
  20. Collins DD, Scoggin CH, Zwillich CW, Weil JV. Hereditary aspects of decreased hypoxic response. J Clin Invest 62: 105–110, 1978.[ISI][Medline]
  21. Conforti L, Bodi I, Nisbet JW, Millhorn DE. O2-sensitive K+ channels: role of the Kv1.2-subunit in mediating the hypoxic response. J Physiol 524: 783–793, 2000.[Abstract/Free Full Text]
  22. Couve A, Restituito S, Brandon JM, Charles KJ, Bawagan H, Freeman KB, Pangalos MN, Calver AR, Moss SJ. Marlin-1, a novel RNA-binding protein associates with GABA receptors. J Biol Chem 279: 13934–13943, 2004.[Abstract/Free Full Text]
  23. Cutz E, Ma TK, Perrin DG, Moore AM, Becker LE. Peripheral chemoreceptors in congenital central hypoventilation syndrome. Am J Respir Crit Care Med 155: 358–363, 1997.[Abstract]
  24. Dauger S, Pattyn A, Lofaso F, Gaultier C, Goridis C, Gallego J, Brunet JF. Phox2b controls the development of peripheral chemoreceptors and afferent visceral pathways. Development 130: 6635–6642, 2003.[Abstract/Free Full Text]
  25. Davenne M, Maconochie MK, Neun R, Pattyn A, Chambon P, Krumlauf R, Rijli FM. Hoxa2 and Hoxb2 control dorsoventral patterns of neuronal development in the rostral hindbrain. Neuron 22: 677–691, 1999.[CrossRef][ISI][Medline]
  26. Deans MR, Gibson JR, Sellitto C, Connors BW, Paul DL. Synchronous activity of inhibitory networks in neocortex requires electrical synapses containing connexin36. Neuron 31: 477–485, 2001.[CrossRef][ISI][Medline]
  27. Doyle DA, Cabral J, Pfuetzner RA, Kuo A, Gulbis JM, Cohen SL, Chait BT, MacKinnon R. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science 280: 69–77, 1998.[Abstract/Free Full Text]
  28. Eisele JH, Wuyam B, Savourey G, Eterradossi J, Bittel JH, Benchetrit G. Individuality of breathing patterns during hypoxia and exercise. J Appl Physiol 72: 2446–2453, 1992.[Abstract/Free Full Text]
  29. Erickson JT, Brosenitsch TA, Katz DM. Brain-derived neurotrophic factor and glial cell line-derived neurotrophic factor are required simultaneously for survival of dopaminergic primary sensory neurons in vivo. J Neurosci 21: 581–589, 2001.[Abstract/Free Full Text]
  30. Evans AM. AMP-activated protein kinase and the regulation of Ca2+ signaling in O2-sensing cells. J Physiol 574: 113–123, 2006.[Abstract/Free Full Text]
  31. Evans AM, Mustard KJW, Wyatt CN, Peers C, Dipp M, Kumar P, Kinnear NP, Hardie DG. Does AMP-activated protein kinase couple inhibition of mitochondrial oxidative phosphorylation by hypoxia to calcium signaling in O2-sensing cells? J Biol Chem 280: 41504–41511, 2005.[Abstract/Free Full Text]
  32. Eyzaguirre C, Fitzgerald RS, Lahiri S, Zapata P. Arterial chemoreceptors. In: Handbook of Physiology, Section 2: The Cardiovascular System, edited by Shepherd JT and Abboud FM. Bethesda, MD: American Physiological Society, 1983, p. 557–621.
  33. Fann MJ, Patterson PH. Depolarization differentially regulates the effects of bone morphogenetic protein (BMP)-2, BMP-6, and activin A on sympathetic neuronal phenotype. J Neurochem 63: 2074–2079, 1994.[ISI][Medline]
  34. Fearon IM, Zhang M, Vollmer C, Nurse CA. GABA mediates autoreceptor feedback inhibition in the rat carotid body via presynaptic GABAB receptors and TASK-1. J Physiol 553: 83–94, 2003.[Abstract/Free Full Text]
  35. Fitzgerald RS, Shirahata M. Systemic responses elicited by stimulating the carotid body: primary and secondary mechamisms. In: Carotid Body Chemoreceptors, edited by Gonzalez C. Berlin, Germany: Springer-Verlag, 1997, p. 171–202.
  36. Fitzgerald RS, Shirahata M, Lahiri S. Amiloride and carotid body chemoreception of hypercapnia and hypoxia. Respir Physiol 81: 337–347, 1990.[CrossRef][ISI][Medline]
  37. Fitzgerald RS. Oxygen and carotid body chemotransduction: the cholinergic hypothesis–a brief history and new evaluation. Respir Physiol 120: 89–104, 2000.[CrossRef][ISI][Medline]
  38. Fuller DD, Bavis RW, Vidruk EH, Wang ZY, Olson EB Jr, Bisgard GE, Mitchell GS. Life-long impairment of hypoxic phrenic responses in rats following 1 month of developmental hyperoxia. J Physiol 538: 947–955, 2002.[Abstract/Free Full Text]
  39. Ganfornina MD, Perez-Garcia MT, Gutierrez G, Miguel-Velado E, Lopez-Lopez JR, Marin A, Sanchez D, Gonzalez C. Comparative gene expression profile of mouse carotid body and adrenal medulla under physiological hypoxia. J Physiol 566: 491–503, 2005.[Abstract/Free Full Text]
  40. Glatt K, Glatt H, Lalande M. Structure and organization of GABRB3 and GABRA5. Genomics 41: 63–69, 1997.[CrossRef][ISI][Medline]
  41. Gonzalez C, Almaraz L, Obeso A, Rigual R. Carotid body chemoreceptors: from natural stimuli to sensory discharges. Physiol Rev 74: 829–898, 1994.[Free Full Text]
  42. Grenningloh G, Pribilla I, Prior P, Multhaup G, Beyreuther K, Taleb O, Betz H. Cloning and expression of the 58 kd beta-subunit of the inhibitory glycine receptor. Neuron 4: 963–970, 1990.[CrossRef][ISI][Medline]
  43. He L, Chen J, Dinger B, Stensaas L, Fidone S. Effect of chronic hypoxia on purinergic synaptic transmission in rat carotid body. J Appl Physiol 100: 157–162, 2006.[Abstract/Free Full Text]
  44. He L, Dinger B, Fidone S. Effect of chronic hypoxia on cholinergic chemotransmission in rat carotid body. J Appl Physiol 98: 614–619, 2005.[Abstract/Free Full Text]
  45. He L, Dinger B, Sanders K, Hoidal J, Obeso A, Stensaas L, Fidone S, Gonzalez C. Effect of p47phox gene deletion on ROS production and oxygen sensing in mouse carotid body chemoreceptor cells. Am J Physiol Lung Cell Mol Physiol 289: L916–L924, 2005.[Abstract/Free Full Text]
  46. Heath D, Smith P, Fitch R, Harris P. Comparative pathology of the enlarged carotid body. J Comp Pathol 95: 259–271, 1985.[CrossRef][ISI][Medline]
  47. Huang L, Guo H, Hellard DT, Katz DM. Glial cell line-derived neurotrophic factor (GDNF) is required for differentiation of pontine noradrenergic neurons and patterning of central respiratory output. Neuroscience 130: 95–105, 2005.[CrossRef][ISI][Medline]
  48. Iturriaga R, Rey S, Del Rio R. Cardiovascular and ventilatory acclimatization induced by chronic intermittent hypoxia: a role for the carotid body in the pathophysiology of sleep apnea. Biol Res 38: 335–340, 2005.[ISI][Medline]
  49. Iturriaga R, Alcayaga J. Neurotransmission in the carotid body: transmitters and modulators between glomus cells and petrosal ganglion nerve terminals. Brain Res Rev 47: 46–53, 2004.[CrossRef][Medline]
  50. Kameda Y. Mash1 is required for glomus cell formation in the mouse carotid body. Dev Biol 283: 128–139, 2005.[CrossRef][ISI][Medline]
  51. Katz DM. Neuronal growth factors and development of respiratory control. Respir Physiol Neurobiol 135: 155–165, 2003.[CrossRef][ISI][Medline]
  52. Katz DM. Regulation of respiratory neuron development by neurotrophic and transcriptional signaling mechanisms. Respir Physiol Neurobiol 149: 99–109, 2005.[CrossRef][ISI][Medline]
  53. Kawakami Y, Yoshikawa T, Shida A, Asanuma Y, Murao M. Control of breathing in young twins. J Appl Physiol 52: 537–542, 1982.[Abstract/Free Full Text]
  54. Kline DD, Peng YJ, Manalo DJ, Semenza GL, Prabhakar NR. Defective carotid body function and impaired ventilatory responses to chronic hypoxia in mice partially deficient for hypoxia-inducible factor 1alpha. Proc Natl Acad Sci USA 99: 821–826, 2002.[Abstract/Free Full Text]
  55. Koerner P, Hesslinger C, Schaefermeyer A, Prinz C, Gratzl M. Evidence for histamine as a transmitter in rat carotid body sensor cells. J Neurochem 91: 493–500, 2004.[CrossRef][ISI]
  56. Krause JE, Chirgwin JM, Carter MS, Xu ZS, Hershey AD. Three rat preprotachykinin mRNAs encode the neuropeptides substance P and neurokinin A. Proc Natl Acad Sci USA 84: 881–885, 1987.[Abstract/Free Full Text]
  57. Krutzik SR, Ochoa MT, Sieling PA, Uematsu S, Ng YW, Legaspi A, Liu PT, Cole ST, Godowski PJ, Maeda Y, Sarno EN, Norgard MV, Brennan PJ, Akira S, Rea TH, Modlin RL. Activation and regulation of Toll-like receptors 2 and 1 in human leprosy. Nat Med 9: 525–532, 2003.[CrossRef][ISI][Medline]
  58. Kugler P. In situ measurements of enzyme activities in the brain. Histochem J 25: 329–338, 1993.[CrossRef][ISI][Medline]
  59. Lahiri S, Rozanov C, Cherniack NS. Altered structure and function of the carotid body at high altitude and associated chemoreflexes. High Alt Med Biol 1: 63–74, 2000.[CrossRef][Medline]
  60. Lee N, Chen J, Sun L, Wu S, Gray KR, Rich A, Huang M, Lin JH, Feder JN, Janovitz EB, Levesque PC, Blanar MA. Expression and characterization of human transient receptor potential melastatin 3 (hTRPM3). J Biol Chem 278: 20890–20897, 2003.[Abstract/Free Full Text]
  61. Lehnardt S, Massillon L, Follett P, Jensen FE, Ratan R, Rosenberg PA, Volpe JJ, Vartanian T. Activation of innate immunity in the CNS triggers neurodegeneration through a Toll-like receptor 4-dependent pathway. Proc Natl Acad Sci USA 100: 8514–8519, 2003.[Abstract/Free Full Text]
  62. Leitner ML, Wang LH, Osborne PA, Golden JP, Milbrandt J, Johnson EM Jr. Expression and function of GDNF family ligands and receptors in the carotid body. Exp Neurol 191, Suppl 1: S68–S79, 2005.
  63. Lerche C, Scherer CR, Seebohm G, Derst C, Wei AD, Busch AE, Steinmeyer K. Molecular cloning and functional expression of KCNQ5, a potassium channel subunit that may contribute to neuronal M-current diversity. J Biol Chem 275: 22395–22400, 2000.[Abstract/Free Full Text]
  64. Leung PS, Lam SY, Fung ML. Chronic hypoxia upregulates the expression and function of AT(1) receptor in rat carotid body. J Endocrinol 167: 517–524, 2000.[Abstract]
  65. Liu C, Goshu E, Wells A, Fan CM. Identification of the downstream targets of SIM1 and ARNT2, a pair of transcription factors essential for neuroendocrine cell differentiation. J Biol Chem 278: 44857–44867, 2003.[Abstract/Free Full Text]
  66. Lo L, Tiveron MC, Anderson DJ. MASH1 activates expression of the paired homeodomain transcription factor Phox2a, and couples pan-neuronal and subtype-specific components of autonomic neuronal identity. Development 125: 609–620, 1998.[Abstract]
  67. Lopez-Lopez J, Gonzalez C, Urena J, Lopez-Barneo J. Low pO2 selectively inhibits K+ channel activity in chemoreceptor cells of the mammalian carotid body. J Gen Physiol 93: 1001–1015, 1989.[Abstract/Free Full Text]
  68. Maltepe E, Keith B, Arsham AM, Brorson JR, Simon MC. The role of ARNT2 in yumor angiogenesis and the neural response to hypoxia. Biochem Biophys Res Commun 273: 231–238, 2000.[CrossRef][ISI][Medline]
  69. Martin-Caraballo M, Dryer SE. Glial cell line-derived neurotrophic factor and target-dependent regulation of large-conductance KCa channels in developing chick lumbar motoneurons. J Neurosci 22: 10201–10208, 2002.[Abstract/Free Full Text]
  70. Martinez A, Saldise L, Ramirez MJ, Belzunegui S, Zudaire E, Luquin MR, Cuttitta F. Adrenomedullin expression and function in the rat carotid body. J Endocrinol 176: 95–102, 2003.[Abstract]
  71. McMahon HT, Missler M, Li C, Sudhof TC. Complexins: cytosolic proteins that regulate SNAP receptor function. Cell 83: 111–119, 1995.[CrossRef][ISI][Medline]
  72. Missler M, Zhang W, Rohlmann A, Kattenstroth G, Hammer RE, Gottmann K, Sudhof TC. {alpha}-Neurexins couple Ca2+ channels to synaptic vesicle exocytosis. Nature 423: 939–948, 2003.[CrossRef][Medline]
  73. Misura KM, Bock JB, Gonzalez LC Jr, Scheller RH, Weis WI. Three-dimensional structure of the amino-terminal domain of syntaxin 6, a SNAP-25 C homolog. Proc Natl Acad Sci USA 99: 9184–9189, 2002.[Abstract/Free Full Text]
  74. Nurse CA. Neurotransmission and neuromodulation in the chemosensory carotid body. Auton Neurosci 120: 1–9, 2005.[CrossRef][ISI][Medline]
  75. O'Connor L, Strasser A, O'Reilly LA, Hausmann G, Adams JM, Cory S, Huang DC. Bim: a novel member of the Bcl-2 family that promotes apoptosis. EMBO J 17: 384–395, 1998.[CrossRef][ISI][Medline]
  76. Otsubo T, Yamaguchi S, Okumura M, Shirahata M. Differential expression of oxygen sensitivity in voltage-dependent K+ channels in inbred strains of mice. Adv Exp Med Biol 580: 319–324, 2006.[ISI][Medline]
  77. Parmantier E, Lynn B, Lawson D, Turmaine M, Namini SS, Chakrabarti L, McMahon AP, Jessen KR, Mirsky R. Schwann cell-derived Desert hedgehog controls the development of peripheral nerve sheaths. Neuron 23: 713–724, 1999.[CrossRef][ISI][Medline]
  78. Paschen W, Hedreen JC, Ross CA. RNA editing of the glutamate receptor subunits GluR2 and GluR6 in human brain tissue. J Neurochem 63: 1596–1602, 1994.[ISI][Medline]
  79. Pattyn A, Morin X, Cremer H, Goridis C, Brunet JF. The homeobox gene Phox2b is essential for the development of autonomic neural crest derivatives. Nature 399: 366–370, 1999.[CrossRef][Medline]
  80. Peers C. Hypoxic suppression of K+ currents in type I carotid body cells: selective effect on the Ca2+-activated K+ current. Neurosci Lett 119: 253–256, 1990.[CrossRef][ISI][Medline]
  81. Peng YJ, Prabhakar NR. Effect of two paradigms of chronic intermittent hypoxia on carotid body sensory activity. J Appl Physiol 96: 1236–1242, 2004.[Abstract/Free Full Text]
  82. Perez-Garcia MT, Colinas O, Miguel-Velado E, Moreno-Dominguez A, Lopez-Lopez JR. Characterization of the Kv channels of mouse carotid body chemoreceptor cells and their role in oxygen sensing. J Physiol 557: 457–471, 2004.[Abstract/Free Full Text]
  83. Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res 30: e36, 2002.
  84. Prabhakar NR. Oxygen sensing by the carotid body chemoreceptors. J Appl Physiol 88: 2287–2295, 2000.[Abstract/Free Full Text]
  85. Prabhakar NR, Jacono FJ. Cellular and molecular mechanisms associated with carotid body adaptations to chronic hypoxia. High Alt Med Biol 6: 112–120, 2005.[CrossRef][Medline]
  86. Puthalakath H, Villunger A, O'Reilly LA, Beaumont JG, Coultas L, Cheney RE, Huang DCS, Strasser A. Bmf: a proapoptotic BH3-only protein regulated by interaction with the myosin V actin motor complex, activated by anoiki