|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Departamento de Bioquímica y Fisiología y Genética Molecular-IBGM, Universidad de Valladolid-CSIC, 47005 Valladolid, Spain
2 Departamento de Genética, Universidad de Sevilla, 41080 Sevilla, Spain
| Abstract |
|---|
|
|
|---|
(Received 18 April 2005;
accepted after revision 2 May 2005;
first published online 12 May 2005)
Corresponding author D. Sánchez: Departamento de Bioquímica y Biología Molecular y Fisiología, Facultad de Medicina, C/Ramón y Cajal, 7, Universidad de Valladolid, 47005 Valladolid, Spain. Email: lazarill{at}ibgm.uva.es
| Introduction |
|---|
|
|
|---|
We are interested in the process of oxygen sensing by CB chemoreceptor cells. This process is turned on after birth in mammals, and accounts for the respiratory and cardiovascular compensatory adaptations organized in response to both acute (lasting seconds to minutes) hypoxia (Prabhakar, 2000) and the sustained hypoxia experienced during acclimatization to high altitudes (Bisgard, 2000). In spite of the detailed knowledge of the electrophysiological and neurochemical changes in CB function generated by hypoxia (Gonzalez et al. 2003), the molecular nature of the components of the signalling cascade responsible for the hypoxic response is partially known (Prabhakar, 2001; Paulding et al. 2002).
In this work we have undertaken a global survey of hypoxia-dependent transcriptional changes to search for the components of the functional module of proteins participating in the physiological response to low PO2. In order to uncover genes differentially expressed in CB cells, we use the hybridization of CB mRNA with oligonucleotide microarrays of mouse genes. The transcripts are obtained from mice that were exposed to either normoxia (21% O2) or to physiological hypoxia (10% O2) for 24 h. The transcriptional profiles obtained for CB cells under these environmental conditions are compared and subtracted from the profile of control (non-chemoreceptor) AM cells exposed to the same PO2 conditions. The subsequent bioinformatical analysis reveals genes involved in the hypoxic response.
| Methods |
|---|
|
|
|---|
Three experiments were performed. In each, two groups of five mice were subjected to normoxia (21%) or hypoxia (10%), for 24 h in a modified plethysmography cage. The groups subjected to normoxia or hypoxia will hereafter be referred to as control and experimental groups. All animals were maintained on a 12 h-light cycle, with food and water ad libitum during the experimental period.
After the experiment, mice were killed by cervical dislocation and the adrenal medulla and carotid bodies were dissected and stored in RNAlaterTM (Ambion) at 4°C prior to the RNA purification step carried out 24 h later. The organs from each experiment were pooled according to each experimental condition. All the carotid body organs (a total of 10 per experimental group) were used for the RNA preparation. Only one adrenal medulla from each mouse (a total of 5 per group) was used, to compensate for the tissue size difference. The remainder of the adrenal medulla samples were stored at 20°C in RNAlater for future use.
RNA purification, amplification and labelling
RNA was purified from the tissue samples using the RNAqueous kit (Ambion), which included a DNAseI treatment to avoid genomic DNA contamination. The quality and purity of the adrenal medulla RNA was assayed by OD measurement at 260 and 280 nm and by electrophoresis on agarose gels. The RNA yields from carotid body samples were too low for detection on gels (approximately 200 ng/10 carotid bodies). An assumption was made that the quality would be similar to the adrenal medulla sample processed in parallel. The whole carotid body total RNA preparation proceeded directly to the amplification phase. Five hundred nanograms of adrenal medulla total RNA was used for amplification.
Two rounds of linear amplification by in vitro transcription were performed following a protocol based on the Eberwine method (Eberwine et al. 1992). Briefly, double stranded cDNA was synthesized using the total RNA as template primed with an oligonucleotide containing the T7 RNA polymerase recognition site plus an oligo-dT segment in order to select for the polyA-mRNA population. After purification, this cDNA was in vitro transcribed resulting in the first amplification step. A second round of amplification followed, using 500 ng of the amplified RNA from each sample. Then, a first strand of cDNA was synthesized from the amplified RNA using random hexamers as primers. This resulted in an estimated average size of cDNA fragments of 600 bp. The second cDNA strand was primed with the T7-oligodT primer and the fragments produced were subjected to in vitro transcription using a biotinylated mix of nucleotides. Therefore, we were linearly amplifying fragments of similar size from the 3' ends of the mRNAs present in the original sample. The biotinylated RNA was ready for hybridization to the chips.
We used the MessageAmp kit (Ambion) for the first round of amplification and the generation and purification of cDNA from the amplified RNA. The Enzo kit (Affymetrix) was used for the second in vitro transcription RNA amplification and biotin labelling.
Oligonucleotide array hybridization and detection
Fifteen micrograms of the biotinylated and fragmented probes were hybridized to Affymetrix GeneChip Mouse MOE430A arrays at 45°C for 16 h, following the manufacturer's protocols. After washes, the arrays were labelled with anti-biotinstreptavidinphycoerythrin antibody and scanned with an Agilent Gene Scanner G2500A.
Microarray data analysis
The analysis of single array gene expression and the comparative expression between samples was performed with the statistical algorithms implemented in the Affymetrix Microarray Analysis Suite v5.0 (MAS 5.0). Global scaling was carried out for each pair of probe arrays using an arbitrary target intensity of 100. The MAS 5.0 statistically determined the presence of a transcript on a given chip (detection P-value cut-offs
1 and
2 were 0.05 and 0.065, respectively). The detection call (presence, P, versus absence, A) was used to assemble the gene expression profile of each tissue under study.
After robust normalization with a perturbation factor of 1.1, the signal intensity values for each pair of control/experimental arrays in a given experiment were used to estimate a change P-value and its associated change call in gene expression. A quantitative estimation of the change in gene expression was derived from the signal-log ratio value calculated in MAS 5.0. Only genes that showed consistent change calls in their expression in the three experiments were selected for further study.
In parallel to the above MAS 5.0-based analysis, we performed an automatic study in GenePublisher (http://www.cbs.dtu.dk/services/GenePublisher) (Knudsen et al. 2003) in order to compare both types of analysis and extract the most significant changes in gene expression.
Descriptions for genes selected for analysis were obtained from the Affymetrix Analysis Center. Data mining with gene ontology classification of the selected transcripts, and gene ontology comparisons between tissues were carried out in the FatiGO server (http://fatigo.bioinfo.cnio.es/) (Al-Shahrour et al. 2004).
The software PATHWAYASSIST and the ResNetTM database (Stratagene) were used to explore the networks of interactions where the physiological hypoxia-regulated genes are potentially involved.
Real time RT-PCR analysis
Total RNA was extracted from normoxic and hypoxic adrenal medulla with the RNaqueous micro kit (Ambion). The RNA samples were treated with deoxyribonuclease I (Ambion) to avoid contaminating genomic DNA (given the intronless nature of the mouse GAPDH gene). cDNA synthesis was performed with MuLV-RT from 1 µg total RNA. Primer sets and TaqMan probes for two selected control genes were: mGAPDH: 5'- TGTGTCCGTCGTGGATCTG-3', 5'-GATGCCTGCTTCACCACCTT-3' and 5'-FAM-TGGAGAAACCTGCCAAGTATGATGACATCA-BHQ2-3'; mGus: 5'-CAATGGTACCGGCAGCC-3', 5'-AAGCTAGAAGGGACAGGCATGT-3' and 5'-FAM-TACGGGAGTCGGGCCCAGTCTTG-BHQ2-3'. Real-time amplifications were carried out in a Rotor-Gene RG3000 (Corbett Research) thermocycler.
Primer sets and TaqMan probes for experimental genes (Agouti-related protein, Peroxisomal
3, and Serpinh1) were obtained as Assays-on-demand kits Mm00475829 g1, Mm00478725 m1, and Mm438056 m1 from ABI.
Amplifications were performed in 25 µl final volume, using 12.5 µl of Absolute QPCR dUTP mix (ABgene), 0.5 µl of UNG (1 µ
µl1, ABI), 1 µl of the corresponding primer/probe mix and 1 µl of template DNA (70 ng µl1). Cycling conditions were: 2 min at 50°C, 15 min at 95°C, and 40 cycles of 95°C for 15 s followed by 60°C for 1 min. Amplification efficiencies for all the genes studied were determined by serial dilutions of template cDNA. All samples were amplified in triplicates and the mean critical threshold (CT) value for each transcript was normalized to m-Gus because of its constant expression in our microarray experiments, and the lack of a reported regulation by changes in environmental PO2. As reaction efficiencies were similar (differences < 0.1), the 
CT method (Livak & Schmittgen, 2001) was used to study variations in gene expression in the two conditions tested (normoxia and hypoxia).
| Results |
|---|
|
|
|---|
The transcriptional profiles of tissues and cells exposed to low PO2 levels have been already documented using microarray chip technologies (e.g. Beitner-Johnson et al. 2001; Gilbert et al. 2003; Hoshikawa et al. 2003; Koike et al. 2004; Aravindan et al. 2005). However, the low PO2 levels (< 20 mmHg), sometimes strict anoxia, used in some studies trigger a general cellular stress response in every tissue.
We are interested in the transcriptional changes generated by levels of hypoxia that would only generate a response in chemosensory organs. This range of PO2, about 10% (70 mmHg) for CB chemoreceptor cells, is considered a physiological hypoxia (Gonzalez, 1998), and mimics the oxygen levels found by an organism at high altitude. Thus, the selection of this PO2 will allow us to differentially identify the molecular components of the signalling cascade related to the sensitivity to physiological hypoxia.
We have selected for this study a hypoxia level of 10% PO2 that produces compensatory hyperventilation reflexes without generating a broad hypoxic cellular stress. The respiratory response to this PO2 level on C57BL/6J01aHsd mice was verified by plethysmography (results not shown). We selected a sustained hypoxia protocol because it mimics the PO2 stimulus that an organism experiences upon exposure to high altitude (at around 5000 meters above sea level). The selection of 24 h as the time of exposure for the experiments was due to the reported decline in transcription of hypoxia-regulated genes after 24 h of hypoxia treatment (Czyzyk-Krzeska et al. 1994).
The selection of the adrenal medulla as a control tissue for subtractive filtering of the CB array data is based on their common developmental origin, their neural and catecholamine-releasing phenotype, and their differential sensitivity to changes in environmental PO2. During embryonic life, the AM chromaffin cells respond to low PO2 to secrete catecholamines (Slotkin & Seidler, 1988). However, in adult life the AM cells are not sensitive to hypoxia levels in the physiological range for a chemoreceptor (Thompson & Nurse, 2000). Moreover, this hypoxia level fails to induce the transcription of tyrosine hydroxylase (TH), the limiting enzyme of the catecholamine synthesis pathway, in AM and superior cervical ganglion cells, while triggering a consistent increase in transcription in CB chemoreceptor cells (Czyzyk-Krzeska et al. 1992).
General description of microarray parameters
The two rounds of linear amplification rendered around 30 µg of double stranded cDNA, with 260/280 absorbance ratios of 2.022.08 indicating a high level of RNA purity. The preparation of total RNA showed tight and clear bands at the expected size for the ribosomal RNA with no signs of degradation.
The quality for the hybridization signals obtained in the microarrays was evaluated through the standard Affymetrix guidelines. Table 1 shows these parameters. The M/3' values measured in the mouse GAPDH control probe sets point to a good amplification of the transcripts. The 5'/3' ratio was not used since we are amplifying fragments from the 3' end of each mRNA. The percentage of present (P) calls (average 50.2%) is also in the accepted range by Affymetrix. The reliability index for triplicates is 0.98, which indicates an adequate level of reproducibility. This is also evidenced by the number of concordant calls in the triplicates, which averages 84.2%.
|
The detection call assigned by MAS 5.0 revealed 9286 transcripts of the MOE430A chip expressed in the mouse CB, and 9975 transcripts expressed in AM.
The gene expression profiles of the two organs (genes with P calls in the three chips) were studied by gene ontology (GO) classification (as defined by the GO consortium, http://www.geneontology.org; Ashburner et al. 2000). The overall profiles show a similar GO distribution in the three general ontology classes. Figure 1A shows this distribution for CB as an example.
|
Of special interest for this work was an enrichment of CB gene expression in genes related to cell communication and response to stimulus (biological process ontology), and metal-ion binding, protein binding, and receptor activity (molecular function ontology) (Fig. 1B).
Analysis of carotid body gene expression in response to physiological hypoxia
We selected probe sets that consistently showed a change in expression in the three experiments based on the change calls (increase, I; decrease, D; or no change, NC) generated by MAS 5.0 (see Methods). The change calls defined as medium increment (MI) or medium decrement (MD) were considered as NC. Entries with three equivalent change calls were considered for further analysis.
The selection of CB genes showing changes in gene expression under hypoxia treatment provided a list of 404 probe sets that increased their expression, and 391 probe sets that decreased their expression.
Because of the non-chemosensory nature of several cell types present in the whole CB tissue used for our array experiments, we designed a subtractive analysis of the transcriptional changes obtained in response to physiological hypoxia. We selected the probe sets in AM chips with consistent changes in gene expression upon physiological hypoxia. Sixty-three probe sets present in AM cells showed an increased gene expression, while 177 probe sets reduced their expression in response to hypoxia.
A fundamental, though expected, conclusion can be deduced from these numbers when compared to those obtained for CB (Fig. 2A): the adult carotid body is more sensitive than the adrenal medulla to physiological hypoxia, corroborating the adequate stimuli classical theory (Sherrington, 1947), and this sensitivity is reflected in assessable transcriptional changes.
|
The GO classification of the genes regulated specifically in the CB reveals genes mainly related to protein and nucleotide binding, hydrolase activity, and interestingly in the context of oxygen sensing, genes related to metal ion binding and oxidoreductase activity (Fig. 2B).
A comparison of the ontology terms that define the hypoxia up- and down-regulated CB genes demonstrated significant differences (asterisks in Fig. 2B) in genes participating in nucleotide binding, in the structural constituent of the ribosome (both more abundant in the up-regulated gene list), and in transferase activity (more abundant in the down-regulated gene list). A catalogue of hypoxia-regulated CB genes with consistent change calls, change fold > 2.5, and change P-values < 0.001 is shown in Table 2.
|
|
|
Three genes were selected for qRT-PCR experiments based on their consistent change in expression under hypoxia in both organs studied, as determined by the array experiments. Two control genes were selected for standardization of gene expression: glyceraldehyde 3 phosphate dehydrogenase (GAPDH) and ß-glucuronidase (Gus), based on their presumed lack of transcriptional regulation under the experimental conditions tested in our experiments. However, GAPDH, though widely used as a standard for quantitative gene expression studies, is known to change its transcriptional levels under hypoxia (Zhong & Simons, 1999). Our array experiments showed a mild increase in GAPDH expression (1.2-fold) in response to physiological hypoxia.
The lack of change of Gus in response to physiological hypoxia in our arrays led us to select this gene as a control for normalization.
Figure 4 shows the qRT-PCR results, after normalization, in comparison with the transcriptional changes evidenced by the chip arrays. The fold-changes obtained in these experiments are comparable to those obtained in the array experiments for each gene studied. This demonstrates the reliability of the arrays hybridization, and supports the biological information extracted from the data.
|
| Discussion |
|---|
|
|
|---|
We have to take into account that the tissues used for generating the hybridizing probes for our arrays are composed of several cell types. Thus, the transcriptional profiles we have obtained correspond in the CB to those of type I chemoreceptor cells, type II glia-like sustentacular cells, and blood vessel cells. In the AM they represent chromaffin cells, blood vessel cells, and some sympathetic neurones.
Both the CB type II cells and the chromaffin AM cells do not respond to physiological hypoxia (Lopez-Barneo et al. 1988; Thompson & Nurse, 2000). Also, the sympathetic neurones of the AM have not been reported to be oxygen sensors in the hypoxia range applied in our experiments. Thus, only blood and blood vessel cells, some of them known to respond to physiological hypoxia (Weir et al. 2002), might be obscuring the transcriptional changes of CB chemosensory cells. However, the analytical substraction of hypoxia-regulated genes from AM is expected to enrich the search by revealing genes involved in hypoxic transduction in CB cells.
Several laboratories have established that the effects on gene regulation upon environmental oxygen changes are associated with changes in the redox state of the cell. As an example, exposure to 10% O2 of a hypoxia-tolerant fish revealed, in several tissues, up-regulated genes such as lactate dehydrogenase A, enolase, triosephosphate isomerase, haeme-oxygenase-1, ferritin, transferrin, glucose-6-phosphatase, IGFBP-1, Erb-B2, and the MAP-kinase MKP-1 (Gracey et al. 2001). Similarly, both transcriptional and post-transcriptional regulation of several genes have been demonstrated in oxygen-sensing cells (reviewed by Seta & Millhorn, 2004). Genes such as the angiogenic factor VEGF, the dopamine D2 receptor, the potassium channel Kv
1.2, the vascular endothelium EPAS-1 (HIF-2
), erythropoietin, c-fos, Bnip-3 and serpine1 are reported to be up-regulated by very low (12%) O2 levels (Premkumar et al. 2000; Seta & Millhorn, 2004). Many of these changes in expression are linked to the stabilization of the transcription factor HIF-1
, which organizes the early response to extreme decreases in O2 (reviewed by Sharp & Bernaudin, 2004). Also, NF-
B transcription and function is affected by oxygen, and thus by the intracellular redox state (reviewed by Haddad, 2004). A twofold change in TH expression is reached after 1 h exposure to 5% O2 in PC12 cells (Czyzyk-Krzeska et al. 1994). Also, moderate hypoxia levels (10% O2) up-regulate the CB transcription of TH (reviewed by Millhorn et al. 2000; Paulding et al. 2002).
Because of these studies, we searched our microarray CB data for changes in genes known to change their expression under hypoxia. We found consistent changes in the expression level under physiological hypoxia of: TH (1.84-fold increase), ferritin (2.16-fold increase), and triosephosphate isomerase (2.01-fold increase). A group of genes known to increase their transcription in response to hypoxia showed analogous but non-consistent changes in our CB microarray experiments. The genes c-fos, Bnip-3, serpine1, lactate dehydrogenase A, and MKP-1 have two I and one NC change calls in our experiments. However, the transcriptional regulation of these genes occurs under more extreme hypoxic conditions.
The ventilatory acclimation occurring under maintained and moderate hypoxia is associated with a sensitization of CB chemosensory cells. These processes participate in the physiological adaptation of organisms to high-altitude environments. Under these conditions, the CB experiences a hyperplasia due to an increase in the number and size of chemosensory cells (Arias-Stella & Valcarcel, 1976; Laidler & Kay, 1978). Moreover, an increase in expression of Na+ channels and connexin 43, as well as in neurotransmitter release and afferent discharges, are part of the response to sustained hypoxia (Gonzalez-Guerrero et al. 1993; Stea et al. 1995; Chen et al. 2002a).
In this context, and though our experiments are not designed to reveal transcriptional changes related to chronic hypoxia and CB hyperplasia, we looked for regulated genes of potential interest for understanding the process of hypoxic chemotransduction in the CB cells. We focused on two main groups: (i) neurotransmitters and neurotransmitter receptors that have been reported to contribute to the low- PO2 secretory response of the organ, and (ii) membrane ion channels that could participate in the early events of the hypoxic transduction cascade.
Given the prominent role of dopamine in the chemosensory process (Gonzalez et al. 1995), and the reported hypoxia regulation of the enzyme TH (see above), we searched for genes involved in the catecholamine metabolism. As mentioned, TH consistently increased its expression in our arrays, and we found a moderate decrease in expression of the noradrenaline synthesizing enzyme dopamine ß-hydroxylase. This suggests that physiological hypoxia favours the sustained presence of dopamine in the synaptic environment, where it functions as a neurotransmitter on afferent nerve endings. However, the dopamine D2 receptor, which has been reported to be up-regulated by hypoxia (Huey & Powell, 2000; Bairam et al. 2003), did not show expression changes under our hypoxia protocol.
The gene endothelin I (ET-1) enhances the excitability of and has a mitogenic effect on CB glomus cells under chronic hypoxia (10% O2), roles that have been ascribed to an up-regulation of the ET receptor A and ET-1 itself (Chen et al. 2002b,c). Our results show the presence in mouse normoxic CB of ET-1, and of ET-2 and the ET-receptor B at a lower level. However, none of these genes showed regulation upon 24 h of exposure to hypoxia. Notwithstanding, the ET converting enzyme (ECE-1), a zinc metalloprotease involved in the processing of endothelins to biologically active peptides, increases its expression (1.73-fold) under physiological hypoxia. We found a similar effect of hypoxia on angiotensin. This is in agreement with the local angiotensin functional system that has been reported in the chemosensory cells of the CB (reviewed by Leung et al. 2003), which is associated with up-regulation of angiotensin II receptors (AT1 type) upon chronic (728 days) hypoxia (Leung et al. 2000). Angiotensinogen and the angiotensin-converting enzyme (ACE) are expressed in CB type I cells, and the ACE activity also increases under chronic hypoxia (Lam et al. 2004). Our microarray data reveal that ACE increases 1.87-fold specifically in mouse CB in response to 24 h of physiological hypoxia. Also, it is important to mention that these two genes (ECE-1 and ACE) are, in fact, switched-ON in our experiments (see Table 3).
The up-regulation by hypoxia of the enzymes carrying out the proteolytic processing of both ET-1 and angiotensin I suggests a hypoxia-mediated control process that would result in local high levels of active ET-1and angiotensin II in chemosensory cells. Of special interest for chemosensory transduction pathways is the common signalling pathway shared by ET-1 and angiotensin II, involving an IP3-mediated increase of [Ca2+]i that could enhance CB chemosensory cells excitability.
Serotonin generates a higher excitability in CB cells (reviewed by Bisgard, 2000) and is possibly involved in the transduction of the hypoxia stimulus. Our expression profile data reveal that serotonin receptors 1a, 3a, and 2c are expressed in mouse CB, and that the 5HT3a receptor expression is consistently down-regulated (1.55-fold) under 10% O2.
Adenosine receptor A2a is expressed in rat CB chemoreceptor cells (Gauda et al. 2000), and is up-regulated by chronic hypoxia in PC12 cells (Kobayashi et al. 2000). Our data show A2b to be present in normoxic mouse CB, in agreement with our experiments using A2b blockers (S. Conde, A. Obeso & C. González, unpublished observations). However, this receptor appears moderately down-regulated by 10% PO2 in our mouse arrays (two D and one NC change calls).
Also, CB cells show an up-regulation of VEGF and its receptor VEGFR-2 under chronic hypoxia in rat (Chen et al. 2003). Our array data show the presence of VEGF-C and the receptors VEGFR-1, 2, and 4 in normoxic CB. However, only the VEGFR-1 is up-regulated upon 24 h of physiological hypoxia, as is also the case for CB under chronic hypoxia (Tipoe & Fung, 2003).
Regarding the hypoxic regulation of membrane ion channels, our microarray data reveal the up-regulation of the ATP-sensitive potassium channel Kir6.1 and the intermediate-conductance and Ca2+-activated K+ channel Kcnn4.
The up-regulation (1.65-fold) of Kir6.1 is in agreement with knockout mice experiments that have implicated this channel in responses to metabolic stress by hypoxia and glucose deprivation (Seino & Miki, 2004), its transcription being induced by hypoxia (Melamed-Frank et al. 2001). Kcnn4 has been reported as a marker for cellular proliferation in vascular smooth muscle cells (Kohler et al. 2003). This role agrees with the consistent up-regulation (2.16-fold) we detect for Kcnn4 in CB subjected to hypoxia, which in turn could underlie the potent angiogenic effect produced by hypoxia on CB (reviewed by Fung, 2003).
Finally, we also found a consistent decrease (1.56-fold) in expression levels of the potassium channel modulatory subunit Kvß1 (Kcnab1) by 10% O2. Kvß associates with Kv1
subunits to induce fast inactivation in non-inactivating Kv1.1 and Kv1.5 channels, or to accelerate the fast inactivation of Kv1.4. Although these interactions were originally thought to be Kv1 subfamily-specific in heterologous systems, Kvß subunits also interact with Kv
subunits in other subfamilies (Aimond et al. 2005). The data currently available regarding the functional role of Kcnab1 obtained from knock-out mice are contradictory: while Kcnab1-deficient mice show alterations in the firing discharges of hipoccampal neurones due to a reduction in the proportion of fast inactivating current (Giese et al. 1998), targeted disruption of Kcnab1 in myocardial cells results in a decreased membrane expression of Kv4.3 channels and increased membrane expression of Kv2.1 channels without altering Kv1 subunit-encoded currents (Aimond et al. 2005). Clearly, in order to infer the functional role of its hypoxic down-regulation in our preparation, we will need to study the partners of Kcnab1 in CB cells and their cellular location. However, this down-regulation by sustained hypoxia, together with the fact that Kvb1.2 has been reported to confer oxygen sensitivity to Kv4.2 channels in heterologous systems (Perez-Garcia et al. 1999), could reflect an important role of this gene in the regulation of hypoxic sensitivity.
These findings differ from those of our recent report (Perez-Garcia et al. 2004) demonstrating the expression of Kv2 and Kv3 subfamilies in mouse CB chemoreceptor cells, and the modulation of Kv3 by acute hypoxic exposure. Several Kv2 and Kv3 subfamily members were detected in our CB arrays in normoxia, but we found no changes in their expression in response to physiological sustained hypoxia. These results suggest that different sets of channels are transducing the response of CB chemoreceptor cells upon acute versus sustained hypoxia, as we have recently documented for rabbit CB glomus cells (Kääb et al. 2005).
In summary, we present the first reported transcriptional profile of mouse CB in response to physiological sustained hypoxia. An interesting finding of our work is the fact that many of the genes we have identified as hypoxia regulated in CB are, in fact, functionally related to molecules known to be released by type I cells upon low O2 stimulation (see Fig. 4 and Supplemental material). These results also suggest that given our whole-organ approach and the long-term hypoxia exposure used, our microarray experiments are identifying not only the genes directly regulated by oxygen, but other genes regulated by molecules released as a consequence of CB activation. Obviously, these genes can be present both in type I cells and vascular vessels, but the subtractive approach followed in the analysis suggests that their regulation is CB specific and that those genes present in vascular cells are changing their expression as a consequence of paracrine factors released by type I cells in response to physiological hypoxia. Additionally, our results identify a new group of CB ion channels experiencing transcriptional changes under hypoxia. The expression of these channels needs to be demonstrated in chemosensory cells, but their regulation after the analytical subtraction from other hypoxia-exposed tissues suggests they are present in CB neural cells.
The data generated in this analysis will lead to future testing of the role of the candidate genes using neurochemical and electrophysiological experiments at the cellular level. In addition, these experiments could exploit the potential of the mouse model system by using mice carrying mutations that will alter their expression or function.
| Supplemental material |
|---|
|
|
|---|
DOI: 10.1113/jphysiol.2005.088815
http://jp.physoc.org/cgi/content/full/jphysiol.2005.088815/DC1
and contains supplemental material showing other interaction networks among CB genes regulated under hypoxia.
This material can also be found as part of the full-text HTML version available from http://www.blackwell-synergy.com
| References |
|---|
|
|
|---|
Al-Shahrour
F, Diaz-Uriarte
R
&
Dopazo
J (2004). FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics
20, 578580.
Aravindan N, Williams MT, Riedel BJ & Shaw AD (2005). Transcriptional responses of rat skeletal muscle following hypoxia-reoxygenation and near ischaemia-reperfusion. Acta Physiol Scand 183, 367377.[CrossRef][Medline]
Arias-Stella J & Valcarcel J (1976). Chief cell hyperplasia in the human carotid body at high altitudes; physiologic and pathologic significance. Hum Pathol 7, 361373.[Medline]
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM & Sherlock G (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25, 2529.[CrossRef][Medline]
Bairam A, Lajeunesse Y, Joseph V & Labelle Y (2003). Time dependent regulation of dopamine D1- and D2-receptor gene expression in the carotid body of developing rabbits by hypoxia. Adv Exp Med Biol 536, 541547.[Medline]
Beitner-Johnson D, Seta K, Yuan Y, Kim H, Rust RT, Conrad PW, Kobayashi S & Millhorn DE (2001). Identification of hypoxia-responsive genes in a dopaminergic cell line by subtractive cDNA libraries and microarray analysis. Parkinsonism Relat Disord 7, 273281.[CrossRef][Medline]
Bisgard GE (2000). Carotid body mechanisms in acclimatization to hypoxia. Respir Physiol 121, 237246.[CrossRef][Medline]
Chen J, Dinger B, Jyung R, Stensaas L & Fidone S (2003). Altered expression of vascular endothelial growth factor and FLK-1 receptor in chronically hypoxic rat carotid body. Adv Exp Med Biol 536, 583591.[Medline]
Chen
J, He
L, Dinger
B, Stensaas
L
&
Fidone
S (2002a). Chronic hypoxia upregulates connexin43 expression in rat carotid body and petrosal ganglion. J Appl Physiol
92, 14801486.
Chen
J, He
L, Dinger
B, Stensaas
L
&
Fidone
S (2002b). Role of endothelin and endothelin A-type receptor in adaptation of the carotid body to chronic hypoxia. Am J Physiol Lung Cell Mol Physiol
282, L1314L1323.
Chen Y, Tipoe GL, Liong E, Leung S, Lam SY, Iwase R, Tjong YW & Fung ML (2002c). Chronic hypoxia enhances endothelin-1-induced intracellular calcium elevation in rat carotid body chemoreceptors and up-regulates ETA receptor expression. Pflugers Arch 443, 565573.[CrossRef][Medline]
Czyzyk-Krzeska MF, Bayliss DA, Lawson EE & Millhorn DE (1992). Regulation of tyrosine hydroxylase gene expression in the rat carotid body by hypoxia. J Neurochem 58, 15381546.[CrossRef][Medline]
Czyzyk-Krzeska
MF, Furnari
BA, Lawson
EE
&
Millhorn
DE (1994). Hypoxia increases rate of transcription and stability of tyrosine hydroxylase mRNA in pheochromocytoma (PC12) cells. J Biol Chem
269, 760764.
Eberwine
J, Yeh
H, Miyashiro
K, Cao
Y, Nair
S, Finnell
R, Zettel
M
&
Coleman
P (1992). Analysis of gene expression in single live neurons. Proc Natl Acad Sci U S A
89, 30103014.
Fung ML (2003). Hypoxia-inducible factor-1: a molecular hint of physiological changes in the carotid body during long-term hypoxemia? Curr Drug Targets Cardiovasc Haematol Disord 3, 254259.[CrossRef][Medline]
Gauda EB, Northington FJ, Linden J & Rosin DL (2000). Differential expression of A2a, A1-adenosine and D2-dopamine receptor genes in rat peripheral arterial chemoreceptors during postnatal development. Brain Res 872, 110.[CrossRef][Medline]
Giese
KP, Storm
JF, Reuter
D, Fedorov
NB, Shao
LR, Leicher
T, Pongs
O
&
Silva
AJ (1998). Reduced K+ channel inactivation, spike broadening, and after-hyperpolarization in Kvß1.1-deficient mice with impaired learning. Learn Mem
5, 257273.
Gilbert RW, Costain WJ, Blanchard ME, Mullen KL, Currie RW & Robertson HA (2003). DNA microarray analysis of hippocampal gene expression measured twelve hours after hypoxia-ischemia in the mouse. J Cereb Blood Flow Metab 23, 11951211.[Medline]
Gonzalez C (1998). Sensitivity to physiologic hypoxia. In Oxygen Regulation of Ion Channels and Gene Expression, ed. Lopez-Barneo J & Weir EK, pp. 321336. Futura Publishing Co., Armonk, NY, USA.
Gonzalez
C, Almaraz
L, Obeso
A
&
Rigual
R (1994). Carotid body chemoreceptors: from natural stimuli to sensory discharges. Physiol Rev
74, 829898.
Gonzalez C, Lopez-Lopez JR, Obeso A, Perez-Garcia MT & Rocher A (1995). Cellular mechanisms of oxygen chemoreception in the carotid body. Respir Physiol 102, 137147.[CrossRef][Medline]
Gonzalez C, Rocher A & Zapata P (2003). Arterial chemoreceptors: cellular and molecular mechanisms in the adaptative and homeostatic function of the carotid body. Rev Neurol 36, 239254.[CrossRef][Medline]
Gonzalez-Guerrero PR, Rigual R & Gonzalez C (1993). Effects of chronic hypoxia on opioid peptide and catecholamine levels and on the release of dopamine in the rabbit carotid body. J Neurochem 60, 17691776.[CrossRef][Medline]
Gracey
AY, Troll
JV
&
Somero
GN (2001). Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis. Proc Natl Acad Sci U S A
98, 19931998.
Haddad JJ (2004). Oxygen sensing and oxidant/redox-related pathways. Biochem Biophys Res Commun 316, 969977.[CrossRef][Medline]
Hoshikawa
Y, Nana-Sinkam
P, Moore
MD, Sotto-Santiago
S, Phang
T, Keith
RL, Morris
KG, Kondo
T, Tuder
RM, Voelkel
NF
&
Geraci
MW (2003). Hypoxia induces different genes in the lungs of rats compared with mice. Physiol Genomics
12, 209219.
Huey KA & Powell FL (2000). Time-dependent changes in dopamine D2-receptor mRNA in the arterial chemoreflex pathway with chronic hypoxia. Brain Res Mol Brain Res 75, 264270.[Medline]
Kääb S, Miguel-Velado E, Lopez-Lopez JR & Perez-Gaiúa MT (2005). Down regulation of kv 3.4 channels by chronic hypoxia increases acute oxygen sensitivity in rabbit carotid body. J. Physiol. (in press).
Knudsen
S, Workman
C, Sicheritz-Ponten
T
&
Friis
C (2003). GenePublisher: Automated analysis of DNA microarray data. Nucl Acids Res
31, 34713476.
Kobayashi
S, Conforti
L
&
Millhorn
DE (2000). Gene expression and function of adenosine A2A receptor in the rat carotid body. Am J Physiol Lung Cell Mol Physiol
279, L273L282.
Kohler
R, Wulff
H, Eichler
I, Kneifel
M, Neumann
D, Knorr
A, Grgic
I, Kampfe
D, Si
H, Wibawa
J, Real
R, Borner
K, Brakemeier
S, Orzechowski
HD, Reusch
HP, Paul
M, Chandy
KG
&
Hoyer
J (2003). Blockade of the intermediate-conductance calcium-activated potassium channel as a new therapeutic strategy for restenosis. Circulation
108, 11191125.
Koike
T, Kimura
N, Miyazaki
K, Yabuta
T, Kumamoto
K, Takenoshita
S, Chen
J, Kobayashi
M, Hosokawa
M, Taniguchi
A, Kojima
T, Ishida
N, Kawakita
M, Yamamoto
H, Takematsu
H, Suzuki
A, Kozutsumi
Y
&
Kannagi
R (2004). Hypoxia induces adhesion molecules on cancer cells: a missing link between Warburg effect and induction of selectin-ligand carbohydrates. Proc Natl Acad Sci U S A
101, 81328137.
Laidler P & Kay JM (1978). A quantitative study of some ultrastructural features of the type I cells in the carotid bodies of rats living at a simulated altitude of 4300 metres. J Neurocytol 7, 183192.[CrossRef][Medline]
Lam
SY, Fung
ML
&
Leung
PS (2004). Regulation of the angiotensin-converting enzyme activity by a time-course hypoxia in the carotid body. J Appl Physiol
96, 809813.
Leung PS, Fung ML & Tam MS (2003). Renin-angiotensin system in the carotid body. Int J Biochem Cell Biol 35, 847854.[CrossRef][Medline]
Leung PS, Lam SY & Fung ML (2000). Chronic hypoxia upregulates the expression and function of AT1 receptor in rat carotid body. J Endocrinol 167, 517524.[Abstract]
Livak
KJ
&
Schmittgen
TD (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2
CT method. Methods
25, 402408.[CrossRef][Medline]
Lopez-Barneo
J, Lopez-Lopez
JR, Urena
J
&
Gonzalez
C (1988). Chemotransduction in the carotid body: K+ current modulated by PO2 in type I chemoreceptor cells. Science
241, 580582.
Melamed-Frank M, Terzic A, Carrasco AJ, Nevo E, Avivi A & Levy AP (2001). Reciprocal regulation of expression of pore-forming KATP channel genes by hypoxia. Mol Cell Biochem 225, 145150.[CrossRef][Medline]
Millhorn DE, Beitner-Johnson D, Conforti L, Conrad PW, Kobayashi S, Yuan Y & Rust R (2000). Gene regulation during hypoxia in excitable oxygen-sensing cells: depolarization-transcription coupling. Adv Exp Med Biol 475, 131142.[Medline]
Paulding WR, Schnell PO, Bauer AL, Striet JB, Nash JA, Kuznetsova AV & Czyzyk-Krzeska MF (2002). Regulation of gene expression for neurotransmitters during adaptation to hypoxia in oxygen-sensitive neuroendocrine cells. Hypoxia-induced regulation of mRNA stability. Microsc Res Tech 59, 178187.[CrossRef][Medline]
Peers
C (2004). Interactions of chemostimuli at the single cell level: studies in a model system. Exp Physiol
89, 6065.
Perez-Garcia
MT, Colinas
O, Miguel-Velado
E, Moreno-Dominguez
A
&
Lopez-Lopez
JR (2004). Characterization of the Kv channels of mouse carotid body chemoreceptor cells and their role in oxygen sensing. J Physiol
557, 457471.
Perez-Garcia
MT, Lopez-Lopez
JR
&
Gonzalez
C (1999). Kvß1.2 subunit coexpression in HEK293 cells confers O2 sensitivity to kv4.2 but not to Shaker channels. J General Physiol
113, 897907.
Prabhakar
NR (2000). Oxygen sensing by the carotid body chemoreceptors. J Appl Physiol
88, 22872295.
Prabhakar
NR (2001). Oxygen sensing during intermittent hypoxia: cellular and molecular mechanisms. J Appl Physiol
90, 19861994.
Premkumar DR, Adhikary G, Overholt JL, Simonson MS, Cherniack NS & Prabhakar NR (2000). Intracellular pathways linking hypoxia to activation of c-fos and AP-1. Adv Exp Med Biol 475, 101109.[Medline]
Reiner
A, Yekutieli
D
&
Benjamini
Y (2003). Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics
19, 368375.
Seino
S
&
Miki
T (2004). Gene targeting approach to clarification of ion channel function: studies of Kir6.x null mice. J Physiol
554, 295300.
Seta
KA
&
Millhorn
DE (2004). Functional genomics approach to hypoxia signaling. J Appl Physiol
96, 765773.
Sharp FR & Bernaudin M (2004). HIF1 and oxygen sensing in the brain. Nat Rev Neurosci 5, 437448.[CrossRef][Medline]
Sherrington C (1947). The Integrative Action of the Nervous System. Yale University Press, New Haven.
Slotkin TA & Seidler FJ (1988). Adrenomedullary catecholamine release in the fetus and newborn: secretory mechanisms and their role in stress and survival. J Dev Physiol 10, 116.[Medline]
Stea A, Jackson A, Macintyre L & Nurse CA (1995). Long-term modulation of inward currents in O2 chemoreceptors by chronic hypoxia and cyclic AMP in vitro. J Neurosci 15, 21922202.[Abstract]
Thompson RJ & Nurse CA (2000). O2-chemosensitivity in developing rat adrenal chromaffin cells. Adv Exp Med Biol 475, 601609.[Medline]
Tipoe
GL
&
Fung
ML (2003). Expression of HIF-1
, VEGF and VEGF receptors in the carotid body of chronically hypoxic rat. Respir Physiol Neurobiol
138, 143154.[CrossRef][Medline]
Weir EK, Hong Z, Porter VA & Reeve HL (2002). Redox signaling in oxygen sensing by vessels. Respir Physiol Neurobiol 132, 121130.[CrossRef][Medline]
Zhong H & Simons JW (1999). Direct comparison of GAPDH, ß-actin, cyclophilin, and 28S rRNA as internal standards for quantifying RNA levels under hypoxia. Biochem Biophys Res Commun 259, 523526.[CrossRef][Medline]
| Acknowledgements |
|---|
This article has been cited by other articles:
![]() |
A. Balbir, H. Lee, M. Okumura, S. Biswal, R. S. Fitzgerald, and M. Shirahata A search for genes that may confer divergent morphology and function in the carotid body between two strains of mice Am J Physiol Lung Cell Mol Physiol, March 1, 2007; 292(3): L704 - L715. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Acker, J. Fandrey, and H. Acker The good, the bad and the ugly in oxygen-sensing: ROS, cytochromes and prolyl-hydroxylases Cardiovasc Res, July 15, 2006; 71(2): 195 - 207. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Tajima, K. Schonherr, S. Niedling, M. Kaatz, H. Kanno, R. Schonherr, and S. H. Heinemann Ca2+-activated K+ channels in human melanoma cells are up-regulated by hypoxia involving hypoxia-inducible factor-1{alpha} and the von Hippel-Lindau protein J. Physiol., March 1, 2006; 571(2): 349 - 359. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||