Am J Physiol Regul Integr Comp Physiol 289: R1512-R1519, 2005.
First published June 30, 2005; doi:10.1152/ajpregu.00089.2005
0363-6119/05 $8.00
COMPARATIVE AND EVOLUTIONARY PHYSIOLOGY
Gene expression profiling of the long-term adaptive response to hypoxia in the gills of adult zebrafish
David L. M. van der Meer,1
Guido E. E. J. M. van den Thillart,1
Frans Witte,1
Merijn A. G. de Bakker,1
Jaya Besser,2
Michael K. Richardson,1
Herman P. Spaink,2
Jelani T. D. Leito,1 and
Christoph P. Bagowski1
1Departments of Integrative Zoology and 2Molecular Cellular Biology, Institute of Biology, University of Leiden, Leiden, The Netherlands
Submitted 10 February 2005
; accepted in final form 16 June 2005
 |
ABSTRACT
|
|---|
Low oxygen levels (hypoxia) play a role in clinical conditions such as stroke, chronic ischemia, and cancer. To better understand these diseases, it is crucial to study the responses of vertebrates to hypoxia. Among vertebrates, some teleosts have developed the ability to adapt to extremely low oxygen levels. We have studied long-term adaptive responses to hypoxia in adult zebrafish. We used zebrafish that survived severe hypoxic conditions for 3 wk and showed adaptive behavioral and phenotypic changes. We used cDNA microarrays to investigate hypoxia-induced changes in expression of 15,532 genes in the respiratory organs (the gills). We have identified 367 differentially expressed genes of which 117 showed hypoxia-induced and 250 hypoxia-reduced expressions. Metabolic depression was indicated by repression of genes in the TCA cycle in the electron transport chain and of genes involved in protein biosynthesis. We observed enhanced expression of the monocarboxylate transporter and of the oxygen transporter myoglobin. The hypoxia-induced group further included the genes for Niemann-Pick C disease and for Wolman disease [lysosomal acid lipase (LAL)]. Both diseases lead to a similar intra- and extracellular accumulation of cholesterol and glycolipids. The Niemann-Pick C protein binds to cholesterol from internal lysosomal membranes and is involved in cholesterol trafficking. LAL is responsible for lysosomal cholesterol degradation. Our data suggest a novel adaptive mechanism to hypoxia, the induction of genes for lysosomal lipid trafficking and degradation. Studying physiological responses to hypoxia in species tolerant for extremely low oxygen levels can help identify novel regulatory genes, which may have important clinical implications.
microarray analysis; Niemann-Pick disease; lysosomal lipase
FOR MOST ORGANISMS, OXYGEN is essential, and oxygen uptake and metabolism are tightly controlled processes. However, changes in the environment, as well as pathophysiological conditions, can lead to reduced oxygen availability. A number of animal species are able to tolerate severe hypoxic stress. Examples are various fish, frog, and turtle species that can often tolerate situations of anoxia and some snakes and insects that can endure severe hypoxia (16). All situations are evocative of ischemia/reperfusion events that are highly damaging to most mammals and birds. In general, hypoxia even for brief periods, can be detrimental or fatal to humans and most mammals. Insufficient blood supply leads to tissue hypoxia in the brain during stroke, in the heart during acute infarction, and to chronic ischemia (reviewed in Ref. 36). Insufficient oxygen supply can lead to severe necrosis and/or apoptosis in tissues of different organs (13, 19).
In mammals and fish, it was shown that hypoxia, to facilitate better oxygen uptake, leads to an increase in ventilation rate (10, 28). For fish species of two different families, it has also been demonstrated that chronic hypoxia leads to morphological changes of the gills bringing about an increase in gill surface area. In several cichlid species, chronic hypoxia caused elongation of gill filaments and increased the size of the secondary lamellae (6, 46).
In the cyprinid, Crucian carp (Carassius carassius), which is especially well suited to withstand periods of anoxia, hypoxia leads to dramatic and reversible gross morphological changes (39). The secondary lamellae are the primary site of oxygen uptake in the gills, and in Crucian carp, they are embedded in a cell mass. Hypoxia in these fish reduces the cell mass through a combination of increasing apoptosis and decreasing proliferation, and thereby leads to an increase in the lamella surface exposed to water, and thus to an increase in gill surface area (39). However, in other systems, adaptations to hypoxia have been shown to lead to inhibition of apoptosis and to stimulation of the antioxidant defense and of stress responses (3, 5, and 14). These hypoxia-induced phenotypic changes and adaptations are due to changes in gene expression.
Two earlier studies have shed light on gene expression changes induced by hypoxia in the goby fish (Gillichthys mirabilis) and in zebrafish embryos [Danio (Brachydanino) rerio], respectively (15, 42). Ton et al. (43) have profiled immediate early gene expression changes of 4,512 genes in zebrafish embryos exposed to extreme hypoxia. In their study, the conditions used were lethal, and zebrafish embryos cannot survive longer than 24 h under these low oxygen conditions (Ton C, personal communication). Thus a struggle for survival was induced in these embryos, which was intended by the investigators, because they were studying the gene-level mechanism of embryonic survival. However, the conditions also led to an arrest of growth and development. It had been shown (30) earlier that zebrafish embryos under hypoxia are developmentally arrested and under anoxic conditions, embryos arrest in S and G2 cell cycle phases. These earlier studies have clearly shown how gene expression pattern can change in the zebrafish embryo in the immediate early response to hypoxia. However, in our study, we used adult zebrafish instead of embryos and nonlethal instead of lethal conditions and focused our work on a single organ, the zebrafish gills. Thereby, for the first time, we have studied long-term adaptive responses to severely low oxygen levels.
Understanding changes in gene expression in fish exposed to hypoxia could reveal new mechanisms of hypoxia tolerance and shed light on the evolution of this adaptive response in vertebrates. This can have important clinical implications.
 |
MATERIALS AND METHODS
|
|---|
Animal handling and hypoxia treatment.
Adult wild-type zebrafish (Danio rerio) were obtained from a local pet store and handled in compliance with local animal care regulations and standard protocols. Our protocol was approved by the review board of Leiden University in accordance with animal protocols of the government of The Netherlands. Fish were kept at 28°C in aquaria with 12:12-h light-dark cycles and were fed twice daily with commercial flake food. Oxygen levels were gradually decreased in 4 days from 8090% to 60, 40, 20, and the final 10% air saturation. After day 4 the fish were kept for an additional 21 days at 10% air saturation. (At 100% air saturation and 28°C the O2 concentration is 8 mg/l). In parallel, a control group was kept at 8090% air-saturated water. Both groups were kept in identical aquaria of 100 liters. The oxygen level on the hypoxia group was kept constant by a controller (Applikon) connected to an oxygen electrode and solenoid valve inline with an air diffuser. The oxygen level in the tank was kept constant by adding oxygen via the diffuser and thereby compensating the oxygen consumption of the fish.
Gill dissection.
The fish were killed with an overdose of anesthetic (MS-222 tricaine methanesulfonate; Argent Chemical Laboratories). Gill covers were removed, and the branchial arches I to IV were dissected from both sides of each fish. For RNA preparation, the filaments of both hemibranches were removed from the bony elements of each gill arch. For scanning electron microscopy, the gill arches of three fish from each group were left intact and fixed immediately in Karnovsky fixative (4% paraformaldehyde and 2.5% glutaraldehyde in 0.1 M phosphate buffer, pH 7.2) for 4 h at 4°C. After three washes in 0.1 M phosphate buffer, pH 7.2, they were transferred to 70% ethanol.
RNA preparation.
After dissection, gills were homogenized in a Dounce homogenizer using 1 ml Trizol solution (GIBCO-BRL, Life Technologies). The whole lamella was used and freed from the pharyngular bone. For each individual sample, a pool of gills from three different animals was used. After Trizol extraction, total RNA was further purified using RNAeasy columns (Quiagen). RNA samples were analyzed for quality control by Lab-on-a-Chip analysis (Agilent) and on agarose gels.
Microarray construction.
Microarrays of the MWG Biotech 14k Zebrafish Oligonucleotide Set were produced using the facilities of the Leiden Genome Technology Center (http://www.lgtc.nl). Oligonucleotides were dissolved in 150 mM phosphate buffer (pH 8.5) to a concentration of 20 µM. With Omnigrid 100 (Genemachines), oligonucleotides were spotted on Codelink activated slides (Amersham Biosciences) according to the Amersham Codelink protocol and as described elsewhere (26, 40a). Microarrays were spotted according to the MIAME guidelines with 15,532 unique zebrafish oligonucleotides from the MWG oligonucleotide set. Ambions array control oligonucleotides (8 sense oligos) and three custom-designed oligonucleotides were spotted together with the MWG Oligonucleotide Set. Spotted Codelink activated slides were treated with blocking solution and washed according to the manufacturers instructions.
Gene expression analysis.
cDNA was synthesized from 2 µg total RNA and then amplified using Ambions Message Amp kit by ServiceXS (Leiden, the Netherlands). During in vitro transcription from the cDNA template, 5-(3-aminoallyl)-UTP (Ambion) was incorporated. Subsequently, 5 µg of cRNA was used for coupling of Amershams monoreactive Cye 3 or Cye 5 dyes and 1.5 µg of the fluorescent-labeled cRNA was used for the hybridizations as previously described (40a). The Ambion array control RNA spikes (amount per spike varying from 0200 pg) were added during the labeling reaction. Labeled cRNA was purified over RNAeasy columns (Quiagen), and the quality was checked by Agilents Lab-on-a-Chip total RNA nanobiosizing assay. Prehybridization, hybridization (55°C), and scanning were carried out as previously described (26, 40a).
Data analysis.
The local background was subtracted from the fluorescent value of each spot. Feature intensities were extracted from scanned microarray images using GenePix Pro 5.1. (Axon Instruments). Images were visually inspected, and spots from low-quality areas of the array were flagged and excluded from further analysis. Spots were also excluded from analysis if both the combined fluorescent intensity for both channels was < 1.4 times that of the local background and the pixel-by-pixel correlation coefficient of the spot was < 0.4. Fluorescence ratios were normalized as previously described (44). A hierarchical clustering algorithm (11) was applied to those genes that fulfilled all of the following conditions: 1) they were not flagged, 2) they were differentially expressed (equal or higher than 1.5-fold induced or equal or < 2-fold repressed), 3) they had P values
0.02; and they were consistent on at least four of the five arrays. In total, 367 genes matched all of these criteria. The average linkage clustered data were visualized using Treeview (11). Data were submitted to GenBank (NCBI Gene Expression Omnibus Series Accession No. GSE2069; Sample Accession Nos. GSM37552, GSM37577, GSM37578, GSM37579, GSM37580; Platform Accession No. GPL1743; http://www.ncbi.nlm.nih.gov/projects/geo).
Real-time quantitative RT-PCR.
For verification of gene expression, we used quantitative real-time RT-PCR. The Roche Master SYBR Green kit was used for the RT-PCR reactions. The annealing and synthesis temperature was 55°C alternating with 96°C for 45 cycles. Dissociation protocols were used to measure melting curves and control for unspecific signals from the primers. We used 100 ng of total RNA per reaction. A standard curve for
-actin using 1, 5, 10, 100, and 500 ng of total RNA was used for normalization. Samples were measured in the Roche LightCycler. The Primer3 software (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) was used to design primers for short amplicons between 50 and 100 bases. The primers used are shown in Table 1.
Scanning electron microscopy.
Gills stored in 70% ethanol were dehydrated to prepare them for scanning electron microscopy. Gills were dehydrated through an ascending acetone series until reaching 100%. Acetone in the tissue was replaced with liquid CO2 cooled to 10°C in a Critical Point Dryer (Baltec CPD 030). The temperature was raised to 40°C after several rinses with CO2 to saturate the gills completely with CO2. During this temperature rise, the critical point was reached and the CO2 liquid was instantaneously converted to gas (35). Silver glue was used to attach the gills to stubs, and they were then sputter coated (Polaron Equipment, scanning electron microscopy coating unit E5100) with gold for 3 min in a low vacuum/argon environment. In the scanning electron microscopy (Jeol JSM-6400) it is possible to tilt and turn the specimen and inspect and photograph the gills from every angle (35). The acceleration voltage was set to 80 kV and the working distance was between 19 and 22 mm. Pictures were taken of the whole gill arch and the primary and secondary lamellae.
 |
RESULTS AND DISCUSSION
|
|---|
Cellular and organismal responses to hypoxia come in two kinds, those that help to survive short hypoxia exposures and those that are necessary to continue living under hypoxia. While the first kind is based on slowing down or even completely shutting down metabolic processes and also involves induction of cell stress responses (29), the second kind provides the key adaptive mechanisms of long-term survival. Among vertebrates, some teleost fish have developed the ability to adapt to extreme low oxygen levels. Earlier studies have measured the survival time during lethal hypoxia in adult zebrafish (34). It was shown that zebrafish can survive a PO2 of 15 Torr (10% air saturation; 0.8 mg/l) for 2 days, whereas immediate exposure to PO2 of 8 Torr (5% air saturation; 0.4 mg/l) was lethal within 12 h (34). It was further shown that preconditioning to hypoxia lowered the mortality rate (34).
For our study, we have used adult zebrafish and exposed them to hypoxia. We were interested in the adaptive response to hypoxia and gradually lowered the PO2 to a nonlethal dose of 10% air saturation (0.8 mg/l). At this oxygen level, zebrafish are still able to grow and gain weight and they can survive these conditions for longer than 6 mo (data not shown). We observed a decrease in activities, a complete lack of rapid turns and movements, and an increased ventilation rate (Supplemental data for this article may be found at http://ajpregu.physiology.org/cgi/content/full/00089.2005/DCI); shown are control fish (normoxia) or zebrafish, which have been under hypoxic conditions for 14 days (hypoxia). We have focused our initial studies on the respiratory organs, the zebrafish gills.
We were particularly interested in gene expression changes associated with the adaptive response to hypoxia in vertebrates. We used cDNA microarrays to profile gene expression changes in response to hypoxia in gills of adult zebrafish. We have identified 367 genes that were differentially expressed under hypoxic conditions. The majority of differentially regulated genes showed a decrease in gene expression (68.1% or 250 genes in total) compared with genes that showed increased expression levels (31.9% or 117 genes in total). Criteria used for differential expression were
1.5 induced and twofold reduced. If criteria for both cases were twofold induced, 286 genes showed differential expression, including 36 induced genes. In our study, we used the same criteria employed in earlier studies on zebrafish embryos (43). All relevant data are submitted (NCBI Gene Expression Omnibus).
Earlier studies (6, 46) have shown that hypoxia causes rapid phenotypic changes in gills of species of two other fish families leading to an increase in surface area. Using scanning electron microscopy, we found different phenotypic changes in zebrafish gills, which could potentially have a similar effect. We observed an increase in cell number and a decrease in cell size in epithelial cells of the secondary lamellae (Fig. 1). The observed cells also appeared to be more ruffled under hypoxic conditions compared with gills isolated from the normoxic control group (Fig. 1). This could lead to an increase in effective surface area of the gills and potentially increase the oxygen extraction capacity at low PO2.

View larger version (171K):
[in this window]
[in a new window]
|
Fig. 1. Scanning electron microscopy pictures of gills from normoxia control zebrafish (left) and hypoxia-treated fish (right). Magnification for the 2 pictures on top, for pictures in the middle, and for pictures on the bottom was comparable. The scale indicated is only valid for the two bottom pictures. The top pictures show the primary lamellae, whereas higher magnifications only show secondary lamellae (middle and bottom). Under hypoxic conditions, cells seen on the surface of the secondary lamellae appear smaller and the surface appears to be more ruffled. Similar data were seen in gills of 3 independent preparations. Several gill arches were investigated for each preparation.
|
|
The observed increase in cell number in the secondary lamellae in the gills could be explained by changes in the gene expression we observed for genes important for cell growth and cell division (Fig. 2G). The protooncogene neuroblastoma-ras is a membrane-bound small G protein often coupled to growth factor receptors. Ras serves as a molecular switch, converting signals from the cell membrane to the nucleus leading to protein synthesis and cell proliferation (1). Hypoxia not only induces neuroblastoma-ras gene expression in the gills (Fig. 2G), it also leads to a decrease of expression of the ras inhibitor p70 ste5 (Ref. 25; Fig. 2G). Furthermore, hypoxia increases expression of cyclin ania-6a (cyclin L; Fig. 2G), which plays a role in promoting cell cycle entry and is commonly overexpressed in primary tumors (44a). These changes in gene expression indicate cell proliferation in the zebrafish gills under long-term hypoxia.

View larger version (73K):
[in this window]
[in a new window]
|
Fig. 2. Selected genes from a total of 367 genes found to be differentially expressed under hypoxic conditions were grouped into categories based on their biological function or their biological properties. The cluster analysis of the following categories is shown: metabolism (A), translation (B), disease related (C), stress response (D), cytoskeleton and cell architecture (E), apoptosis (F), growth regulation (G), phosphatases (H), immune response (I), Degradation (J), chaperones (K), and transport (L). Quantitative changes in gene expression are induced genes (red), and repressed genes (green). Results of the individual slides are in vertical columns 1-5 (starting at the left of each cluster). Vertical column 6 (far right column of the cluster) represents the mean. Gene expression changes had to be consistent on at least 4 of 5 arrays. Each slide represents a comparison of normoxia vs. hypoxia samples. Hypoxia samples were labeled with Cye 3 and normoxia samples with Cye 5 (1, 3, 4, 5) except for number 2, which was a dye swap experiment (corresponding to number 5) (the color indication was later changed so that it is consistent between all 6 vertical columns). The complete annotated data set is available at NCBI Geo Gene Expression Omnibus (see MATERIALS AND METHODS for details) and a complete list of the 367 differential expressed genes and an expression profile (supplemental data for this article may be found at http://ajpregu.physiology.org/cgi/content/full/00089.2005/DCI).
|
|
Fish respond to low oxygen levels by activation of anaerobic glycolysis and, in most cases, by metabolic depression and energy conservation (7, 41, 44b, 48). This is reflected by gene expression changes observed earlier in zebrafish embryos exposed to hypoxia (43). We observed similar gene expression changes in the gills of the adult fish. Genes of the TCA cycle and the electron transfer chain that are involved in mitochondrial energy production were found to be repressed. No genes involved in these processes were found to be induced (Fig. 2A). In contrast, genes involved in glycolysis, like glycogen phosphorylase and aldolase, showed increased expression under hypoxic conditions (Fig. 2A). These findings are not novel, but indicate that some gene expression changes are comparable between the long-term hypoxia response in adult zebrafish gills and the immediate response in embryos.
We discovered in the gills that hypoxic conditions reduced gene expression of components of fat metabolism, uptake, and transport. Examples are acyl-CoA dehydrogenase, intestinal fatty acid binding protein, and oocyte-type fatty acid binding protein, respectively (Fig. 2A).
-oxidation consumes a lot of oxygen, and decreasing
-oxidation could help to save oxygen. It is well known that ischemia, due to the accumulation of intermediates of the
-oxidation, causes major cell damage and leads to heart failure (27).
One of the most energy-consuming cellular processes is protein biosynthesis. In the goby fish (Gillichthys mirabilis), and in zebrafish embryos, hypoxia leads to repression of genes encoding components of the protein translation machinery (15, 43). Likewise, in our studies, ribosomal proteins and their subunits are the largest group of hypoxia-repressed genes. No ribosomal protein was found to be induced by hypoxia (Fig. 2B). Conversely, we found that hypoxia stimulated gene expression of two translational elongation factors, eef-1a and eef-2 (Fig. 2B). Roles other than translational control have been assigned to eef-1a, and it is possible that this is true for eef-2 as well (23). We further found that hypoxia decreased gene expression of the eukaryotic initiation factor-2
subunit kinase (PEK) (Fig. 2B). In response to environmental stresses, PEK has been shown to reduce protein synthesis by phosphorylation of the
-subunit of eukaryotic translation initiation factor-2
(37). It is possible that PEK is not directly involved in hypoxia-induced repression of translation. PEK had also been shown to be involved in the cellular stress response (40). Another possible explanation for this observation is that the reduction in PEK expression does not correlate with a decrease in PEK activity.
An organism uses many different strategies to adjust to hypoxic conditions. We found several novel hypoxia-dependent gene expression changes that point to physiological adaptations to low oxygen levels previously not described. An example is the monocarboxylate transporter (mct4), which transports metabolites like pyruvate and lactate. In our study, we found that hypoxia increased gene expression of mct4 (Fig. 2L). If oxygen supplies do not meet demands, enhanced anaerobic metabolism leads to accumulation of lactate, a phenomenon best known in muscle tissue. The monocarboxylate transporter can help to remove accumulated lactate and transport it to the bloodstream, particularly to the heart and other aerobic tissues. Further studies are necessary to show that increased transport of monocarboxylates is indeed an important adaptive mechanism to chronic hypoxia.
Another important physiological adaptation to hypoxia is indicated by our findings of enhanced expression of the oxygen transporter myoglobin (Fig. 2L). Myoglobin increases the oxygen diffusion rate through tissues (38). Thus by increasing myoglobin levels, the organism could compensate for the gradient decline and stabilize the oxygen flux to its tissues.
A novel mechanism of adapting to low oxygen levels is suggested by our findings that genes causative for two genetic metabolic disorders affecting cholesterol trafficking and degradation showed enhanced expression under hypoxia (Fig. 2C). One is the Niemann-Pick disease gene C (NPC) and the other is lysosomal acid lipase (LAL; cholesterol esterase), which is linked to Wolman disease. The Niemann-Pick disease is an inherited metabolic disorder. NPC is a neurovisceral disorder characterized by accumulation of cholesterol and glycolipids in the lysosomal/late endosomal system, which is due to mutations on either the NPC1 or the NPC2 genes (31). Although NPC1 and NPC2 proteins appear essential for proper cellular cholesterol trafficking, their precise functions and relationship have remained elusive (31).
LAL is essential for the hydrolysis of cholesteryl esters and triglycerides in the lysosome (49). Defective LAL activity in Wolman disease results in large amounts of lipids, particularly cholesteryl esters and glycerides, to accumulate in the brain, liver, spleen, lymph nodes, and other tissues (36). Strikingly, both Niemann-Pick and Wolman disease lead to accumulation of cholesteryl esters and glycerides in several tissues, and both are connected to lysosomes. Lysosomes represent the machinery for cellular turnover. We speculate that both NPC and LAL are involved in biomembrane turnover, transport, or generation of nonmetabolizable lipids, such as cholesterol and sphingosine. Thus upregulation may be expected in active tissues and tissues exposed to wear and tear, especially in their membranes. Hypoxia would lead to a higher need for biomembranes possibly due to higher wear and tear because of the increase in ventilation rate. This hypothesis could explain our observation of ruffled cell membranes on the surface of the secondary lamellae (Fig. 1). Both genes have so far not been implicated in a response to hypoxia. We have confirmed our microarray results for both genes by quantitative PCR (Fig. 3). The upregulation of these genes corroborates clearly with the above-mentioned downregulation of genes for lipid oxidation. Additional studies are necessary to further clarify the role of the Niemann-Pick disease gene and of LAL in adaptation to hypoxia.
Hypoxic conditions represent high stress levels that can lead to the induction of genes involved in the cellular stress response, like heat shock proteins (HSPs). Ischemia-induced expression of HSP70 has been demonstrated in tissues of intact animals (8) and in mammalian cell culture cells (2). HSP70 functions as a chaperone and is also known to protect cells against apoptosis (18). We observed increased expression of genes for HSP70 and HSP40 (Fig. 2K). We also observed decreased expression for other chaperones like chaperonin 10 (Fig. 2K), which is involved in the folding of newly synthesized proteins (12). Because protein biosynthesis appears to be lowered, it fits that chaperonin 10 expression is also downregulated. In contrast to HSP70 and HSP40, other genes involved in the cellular stress response showed decreased expression (Fig. 2D). This could be due to the time point we have chosen; different stress proteins are switched on and off at different times during exposure to stress. We investigated the long-term response to hypoxia. A future study including earlier time points could shed light on the temporal induction of stress response and other genes.
As pointed out in the introduction, anoxia induces apoptosis in the gills of the Crucian carp leading to an increase in gill surface area. We have not found any indication of hypoxia-induced cell death neither by apoptosis nor necrosis in zebrafish gills (data not shown). We have discovered expression changes for several genes involved in apoptosis (Fig. 2F). However, the majority of the observed gene regulations tend to be antiapoptotic.
We discovered that hypoxia induced expression levels of four different phosphatases, whereas no phosphatase was found with repressed expression (Fig. 2H). This was similar for the clustering of disease-related genes, which exclusively showed enhanced gene expression (Fig. 2C and Table 2).
View this table:
[in this window]
[in a new window]
|
Table 2. Different disease-related genes with enhanced expression upon long-term hypoxia treatment in the zebrafish gills
|
|
Clustering for genes of the immune response revealed that this group solely showed a decrease in gene expression under hypoxic conditions (Fig. 2I). However, these observations should not lead to general conclusions, because the numbers of phosphatases, disease-related genes, and immune response genes found are not large enough.
Several of the genes we found differentially expressed in the zebrafish had been described as hypoxia-responsive genes in other systems. To further verify our results, we used quantitative real-time PCR for nine transcripts. We confirmed the gene expression changes of these nine transcripts found in the microarray studies by this independent method (Fig. 3). The fold induction values were not directly comparable, but in all cases, induction or reduction was confirmed. Quantitative differences between array data and qPCR results have been reported previously (26, 43).
Our study has shown that cDNA microarray analysis is well suited to studying the responses to physiological stresses, such as low oxygen levels. In this study, we have identified specific changes in gene expression in the gills of adult zebrafish and identified possible novel mechanisms for the long-term adaptive response to hypoxia. We identified several genes that currently have not been linked to hypoxia. It remains to be determined which of these changes are tissue specific. We are currently investigating adaptive gene expression changes to long-term exposure to hypoxia in zebrafish heart and muscle tissues. Future studies will be very useful to obtain a global view of tissue- and species-specific gene expression changes and to delineate the role of the individual genes found. This will lead to a better understanding of the adaptive responses to hypoxia and findings in zebrafish might have clinical implications in humans.
 |
GRANT
|
|---|
This work was supported by European Commission Sixth Framework Prorgamme Grant (contract LSHG-CT-2003503496, ZF-Models).
 |
ACKNOWLEDGMENTS
|
|---|
We thank Tim Gant (Medical Research Council Toxicology Unit, Leicester UK) for his help with data analysis and software, Robert Geisler (Max-Planck-Institute, Tuebingen) for help with the supply of the oligonucleotide library, Andre Wijfjes for advice on the microarrays, and Duddy Oyib for critical reading of the manuscript.
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: C. P. Bagowski, Institute of Biology, Univ. of Leiden, Wassenaarseweg 64, 2333 AL Leiden, The Netherlands (e-mail: bagowski{at}rulbim.leidenuniv.nl)
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
|
|---|
- Adjei AA. Blocking oncogenic Ras signaling for cancer therapy. J Natl Cancer Inst 93: 10621074, 2001.[Abstract/Free Full Text]
- Benjamin IJ, Kroger B, and Williams S. Activation of the heat shock transcription factor by hypoxia in mammalian cells. Proc Natl Acad Sci USA 87: 62636267, 1990.[Abstract/Free Full Text]
- Blagosklonny MV. Hypoxia-inducible factor: Achilles heel of antiangiogenic cancer therapy. Int J Oncol 19: 257262, 2001.[ISI][Medline]
- Brown SC, Torelli S, Brockington M, Yuva Y, Jimenez C, Feng L, Anderson L, Ugo I, Kroger S, Bushby K, Voit T, Sewry C, and Muntoni F. Abnormalities in
-dystroglycan expression in MDC1C and LGMD2I muscular dystrophies. Am J Pathol 164: 727737, 2004.[Abstract/Free Full Text]
- Brusselmans K, Bono F, Collen D, Herbert JM, Carmeliet P, and Dewerchin M. A novel role for VEGF as an autocrine survival factor for embryonic stem cells during hypoxia. J Biol Chem 280: 34933499, 2005.[Abstract/Free Full Text]
- Chapman LJ, Galis F, And Shinn J. Phenotypic plasticity and the possible role of genetic assimilation: hypoxia-induced trade-offs in the morphological traits of an African cichlid. Ecol Lett 3: 387393, 2000.[CrossRef]
- Dalla Via J, van den Thillart GEEJM, Cattani O, and de Zwaan A. Influence of long-term hypoxia exposure on the energy metabolism of Solea solea. II. Intermediary metabolism in blood, liver, and muscle. Mar Ecol 111: 1727, 1994.
- Dillmann WH, Mehta HB, Barrieux A, Guth BD, Neeley WE, and Ross J Jr. Ischemia of the dog heart induces the appearance of a cardiac mRNA coding for a protein with migration characteristics similar to heat-shock/stress protein 70. Circ Res 59: 110114, 1986.[Abstract/Free Full Text]
- Du H, Sheriff S, Bezerra J, Leonova T, and Grabowski GA. Molecular and enzymatic analyses of lysosomal acid lipase in cholesteryl ester storage disease. Mol Genet Metab 64: 126134, 1998.[CrossRef][ISI][Medline]
- Duffin J and Mahamed S. Adaptation in the respiratory control system. Can J Physiol Pharmacol 81: 765773, 2003.[CrossRef][ISI][Medline]
- Eisen MB, Spellman PT, Brown PO, and Botstein D. Cluster analysis and display of genome-wide expression patterns. Genetics 95: 1486314868, 1998.
- Fenton WA and Horwich AL. Chaperonin-mediated protein folding: fate of substrate polypeptide. Q Rev Biophys 36: 229256, 2003.[CrossRef][ISI][Medline]
- Geng YJ. Molecular mechanisms for cardiovascular stem cell apoptosis and growth in the hearts with atherosclerotic coronary disease and ischemic heart failure. Ann NY Acad Sci 1010: 687697, 2003.[Abstract/Free Full Text]
- Giaccia AJ. Hypoxic stress proteins: survival of the fittest. Semin Radiat Oncol 6: 4658, 1996.[CrossRef][ISI][Medline]
- Gracey AY, Troll JV, and Somero GN. Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis. Proc Natl Acad Sci USA 98: 19931998, 2001.[Abstract/Free Full Text]
- Hermes-Lima M and Zenteno-Savin T. Animal response to drastic changes in oxygen availability and physiological oxidative stress. Comp Biochem Physiol C Toxicol Pharmacol 133: 537556, 2002.[CrossRef][ISI][Medline]
- Hohfeld J. Regulation of the heat shock conjugate Hsc70 in the mammalian cell: the characterization of the anti-apoptotic protein BAG-1 provides novel insights. J Biol Chem 3: 269274, 1998.
- Ikonomidou C and Turski L. Excitotoxicity and neurodegenerative diseases. Curr Opin Neurol 8: 487497, 1995.[CrossRef][ISI][Medline]
- Jaaskelainen P, Miettinen R, Karkkainen P, Toivonen L, Laakso M, and Kuusisto J. Genetics of hypertrophic cardiomyopathy in Eastern Finland: few founder mutations with benign or intermediary phenotypes. Ann Med 36: 2332, 2004.[CrossRef][ISI][Medline]
- Janckila AJ, Nakasato YR, Neustadt DH, and Yam LT. Disease-specific expression of tartrate-resistant acid phosphatase isoforms. J Bone Miner Res 18: 19161919, 2003.[CrossRef][ISI][Medline]
- Kato M and Dobyns WB. Lissencephaly and the molecular basis of neuronal migration. Hum Mol Genet 12: R89R96, 2003.[Abstract/Free Full Text]
- Lamberti A, Caraglia M, Longo O, Marra M, Abbruzzese A, and Arcari P. The translation elongation factor 1A in tumorigenesis, signal transduction and apoptosis. Amino Acids 26: 443448, 2004.[ISI][Medline]
- Lombard-Platet G, Savary S, Sarde CO, Mandel JL, and Chimini G. A close relative of the adrenoleukodystrophy (ALD) gene codes for a peroxisomal protein with a specific expression pattern. Proc Natl Acad Sci USA 93: 12651269, 1996.[Abstract/Free Full Text]
- Majidi M, Gutkind JS, and Lichy JH. Deletion of the COOH terminus converts the ST5 p70 protein from an inhibitor of RAS signaling to an activator with transforming activity in NIH-3T3 cells. J Biol Chem 275: 65606565, 2000.[Abstract/Free Full Text]
- Meijer AH, Verbeek FJ, Salas-Vidal E, Corredor-Adamez M, Bussman J, van der Saar A, Otto GW, Geisler R, and Spaink HP. Transcriptome profiling of adult zebrafish at the late stage of chronic tuberculosis due to Mycobacterium marinum infection. Mol Immunol 42: 185203, 2005.
- Mitchell GA, Kassovska-Bratinova S, Boukaftane Y, Robert MF, Wang SP, Ashmarina L, Lambert M, Lapierre P, and Potier E. Medical aspects of ketone body metabolism. Clin Invest Med 18: 193216, 1995.[ISI][Medline]
- Mortola JP. Implications of hypoxic hypometabolism during mammalian ontogenesis. Respir Physiol Neurobiol 141: 345356, 2004.[CrossRef][ISI][Medline]
- Ning W, Chu TJ, Li CJ, Choi AM, and Peters DG. Genome-wide analysis of the endothelial transcriptome under short-term chronic hypoxia. Physiol Genomics 18: 7078, 2004.[Abstract/Free Full Text]
- Padilla PA and Roth MB. Oxygen deprivation causes suspended animation in the zebrafish embryo. Proc Natl Acad Sci USA 98: 73317335, 2001.[Abstract/Free Full Text]
- Patterson MC and Platt F. Therapy of Niemann-Pick disease, type C. Biochim Biophys Acta 1685: 7782, 2004.[Medline]
- Redon R, Hussenet T, Bour G, Caulee K, Jost B, Muller D, Abecassis J, and du Manoir S. Amplicon mapping and transcriptional analysis pinpoint cyclin L as a candidate oncogene in head and neck cancer. Cancer Res 62: 62116217, 2002.[Abstract/Free Full Text]
- Rees BB, Sudradjat FA, and Love JW. Acclimation to hypoxia increases survival time of zebrafish, Danio rerio, during lethal hypoxia. J Exp Zool 289: 266272, 2001.[CrossRef][ISI][Medline]
- Robinson G and Gray T. Electron Microscopy 3: Theory and Practice of Histology Techniques (3rd ed.). Edited by Bancroft JD and Stevens A. London: Churchill Livingstone, 1990, p. 563581.
- Semenza GL, Agani F, Feldser D, Iyer N, Kotch L, Laughner E, and Yu A. Hypoxia, HIF-1, and the pathophysiology of common human diseases. Adv Exp Med Biol 475: 123130, 2000.[ISI][Medline]
- Shi Y, Vattem KM, Sood R, An J, Liang J, Stramm L, and Wek RC. Identification and characterization of pancreatic eukaryotic initiation factor-2
subunit kinase, PEK, involved in translational control. Mol Cell Biol 18: 74997509, 1998.[Abstract/Free Full Text]
- Sidell BD. Intracellular oxygen diffusion: the roles of myoglobin and lipid at cold body temperature. J Exp Biol 201: 11191128, 1998.[Abstract]
- Sollid J, De Angelis P, Gundersen K, and Nilsson GE. Hypoxia induces adaptive and reversible gross morphological changes in crucian carp gills. J Exp Biol 206: 36673673, 2003.[Abstract/Free Full Text]
- Sood R, Porter AC, Ma K, Quilliam LA, and Wek RC. Pancreatic eukaryotic initiation factor-2
kinase (PEK) homologues in humans, Drosophila melanogaster and Caenorhabditis elegans that mediate translational control in response to endoplasmic reticulum stress. Biochem J 346: 281293, 2000.
- t Hoen PA, de Kort F, van Ommen GJ, and den Dunnen JT. Fluorescent labelling of cRNA for microarray applications. Nucleic Acids Res 31: e20e30, 2003.[Abstract/Free Full Text]
- Ton C, Stamatiou D, Dzau VJ, and Liew CC. Construction of a zebrafish cDNA microarray: gene expression profiling of the zebrafish during development. Biochem Biophys Res Commun 296: 11341142, 2002.[CrossRef][ISI][Medline]
- Ton C, Stamatiou D, and Liew CC. Gene expression profile of zebrafish exposed to hypoxia during development. Physiol Genomics 13: 97106, 2003.[Abstract/Free Full Text]
- Turton NJ, Judah DJ, Riley J, Davies R, Lipson D, Styles JA, Smith AG, and Gant TW. Gene expression and amplification in breast carcinoma cells with intrinsic and acquired doxorubicin resistance. Oncogene 20: 13001318, 2001.[CrossRef][ISI][Medline]
- Van den Thillart GEEJM, Dalla Via J, Vitali G, and Cortesi P. Influence of long-term hypoxia exposure on the energy metabolism of Solea solea. I. Critical O2 levels for aerobic and anaerobic metabolism. Mar Ecol 104: 109117, 1994.
- Van Raaij MT, Pit DS, Balm PH, Steffens AB, and van den Thillart GEEJM. Behavioral strategy and the physiological stress response in rainbow trout exposed to severe hypoxia. Horm Behav 30: 8592, 1996.[CrossRef][Medline]
- Wallgren-Pettersson C, Donner K, Sewry C, Bijlsma E, Lammens M, Bushby K, Giovannucci Uzielli ML, Lapi E, Odent S, Akcoren Z, Topaloglu H, and Pelin K. Mutations in the nebulin gene can cause severe congenital nemaline myopathy. Neuromuscul Disord 12: 674679, 2002.[CrossRef][ISI][Medline]
- Witte F, Wanink JH, Rutjes HA, van der Meer HJ, and van den Thillart GEEJM. Eutrophication and the fish fauna of Lake Victoria. In: Tropical Eutrophic Lakes: Their Restoration and Management. Edited by Vikram Reddy M. New Delhi: Oxford, 2005.
- Zani VJ, Asou N, Jadayel D, Heward JM, Shipley J, Nacheva E, Takasuki K, Catovsky D, and Dyer MJ. Molecular cloning of complex chromosomal translocation t(8,14,12)(q24.1;q323;q241) in a Burkitt lymphoma cell line defines a new gene (BCL7A) with homology to caldesmon. Blood 87: 31243134, 1996.[Abstract/Free Full Text]
- Zhou BS, Wu RS, Randall DJ, and Lam PK. Bioenergetics and RNA/DNA ratios in the common carp (Cyprinus carpio) under hypoxia. J Comp Physiol [B] 171: 4957, 2001.[CrossRef][Medline]
- Zschenker O, Jung N, Rethmeier J, Trautwein S, Hertel S, Zeigler M, and Ameis D. Characterization of lysosomal acid lipase mutations in the signal peptide and mature polypeptide region causing Wolman disease. J Lipid Res 42: 10331040, 2001.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
S. Kalujnaia, I. S. McWilliam, V. A. Zaguinaiko, A. L. Feilen, J. Nicholson, N. Hazon, C. P. Cutler, and G. Cramb
Transcriptomic approach to the study of osmoregulation in the European eel Anguilla anguilla
Physiol Genomics,
November 14, 2007;
31(3):
385 - 401.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
R. N. Morrison, G. A. Cooper, B. F. Koop, M. L. Rise, A. R. Bridle, M. B. Adams, and B. F. Nowak
Transcriptome profiling the gills of amoebic gill disease (AGD)-affected Atlantic salmon (Salmo salar L.): a role for tumor suppressor p53 in AGD pathogenesis?
Physiol Genomics,
September 14, 2006;
26(1):
15 - 34.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
A. F. Riggs and T. A. Gorr
A globin in every cell?
PNAS,
February 21, 2006;
103(8):
2469 - 2470.
[Full Text]
[PDF]
|
 |
|
Copyright © 2005 by the American Physiological Society.