|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
COMPARATIVE AND EVOLUTIONARY PHYSIOLOGY
1Department of Zoology, The University of British Columbia, Vancouver, British Columbia, Canada; and 2Department of Animal Science, University of California Davis, Davis, California
Submitted 16 March 2008 ; accepted in final form 23 June 2008
| ABSTRACT |
|---|
|
|
|---|
fG'). The activity of PDH was rapidly (within 3 h) decreased by
50% upon hypoxia exposure and remained depressed relative to normoxic samples throughout. Inactivation of PDH was primarily mediated via posttranslational modification following the accumulation of acetyl-CoA and subsequent activation of pyruvate dehydrogenase kinase (PDK). Estimated changes in cytoplasmic and mitochondrial [NAD+]/[NADH] did not parallel one another, suggesting the mitochondrial NADH shuttles do not function during hypoxia exposure. Large increases in the expression of PDK (PDK isoform 2) were consistent with decreased PDH activity; however, these changes in mRNA were not associated with changes in total PDK-2 protein content assessed using mammalian antibodies. No other changes in the expression of other known hypoxia-responsive genes (e.g., lactate dehydrogenase-A or -B) were observed in either muscle or liver. Gibbs free energy; pyruvate dehydrogenase kinase; energy charge; Gibbs free energy; muscle; fish
PDH is the rate-limiting enzyme that regulates the rate of entry of glycolytically derived pyruvate into the TCA cycle and mitochondrial oxidative metabolism (20) and is regulated by both product inhibition (NADH and acetyl-CoA) and reversible covalent modification (phosphorylation/dephosphorylation). The transformation of PDH between the active form (PDHa) and inactive form is regulated by the relative activities of pyruvate dehydrogenase kinase (PDK), which phosphorylates PDH to deactivate it, and PDH phosphatase, which activates PDH by dephosphorylation (42). PDH kinase is activated by metabolite concentrations signaling a high cellular energy status (high [acetyl-CoA]/[free-CoA], [ATP]/[ADPfree] and low pyruvate and [NAD+]/[NADH]) (20). Expression of PDK in mammalian cell lines, in particular the PDK-1 isoform, is hypoxia responsive and regulated by the transcription factor hypoxia inducible factor-1 (HIF-1) (10). Increases in PDK-1 mRNA during hypoxia exposure has been linked to the phosphorylation of PDH, reduced mitochondrial oxygen consumption, and ROS generation (33, 43). Expression of PDK-1 in HIF-1
null lymphoma cells resulted in increases in cellular [ATP], an attenuation of hypoxic ROS generation, and a lack of hypoxia-induced apoptosis during severe hypoxia exposure (24). Thus, the regulation of PDK appears to be an important aspect of hypoxia cell survival even in hypoxia-intolerant mammalian cells; however, few studies have examined the regulation of PDH and PDK in animals that have evolved to survive hypoxia exposure. Richards et al. (40) demonstrated a modest decrease in muscle PDH activity in the hypoxia-tolerant Amazonian cichlid (Astronotus ocellatus) but the detailed mechanisms regulating PDH during hypoxia are not known in any hypoxia-tolerant animals.
Here we characterize metabolic responses of the hypoxia-tolerant killifish (Fundulus heteroclitus) to 15-h severe hypoxia exposure followed by 12-h recovery in normoxic water. Killifish inhabit estuaries along the east coast of North America that undergo diel and seasonal fluctuations in dissolved oxygen. Diel fluctuations in oxygen are caused by plant and animal respiration at night when no photosynthesis occurs to replenish dissolved oxygen. As a result of the environmental fluctuations in dissolved oxygen, killifish have evolved an amazing suite of adaptations to enhance hypoxia survival including modifications to behavior (48) and biochemical and molecular responses (25, 37). In F. grandis, long-term hypoxia exposure led to tissue-specific phenotypes in metabolic capacity with an upregulation in glycolytic capacity (i.e., increases in enzyme activities) in liver and decreases in muscle (29). Taken together, these data suggest that muscle severely reduces ATP utilization during hypoxia exposure, in part through reductions in motor activity, and therefore we focused on the muscle in the present study to determine whether PDH modulation was associated with hypoxic survival. Specifically, we examined the activity of PDH in muscle of killifish during hypoxia exposure and examined the factors that regulate PDH activity, including changes in allosteric modulators (e.g., [ATP]/[ADPfree], [NAD+]/[NADH] ratios, [pyruvate]) and the patterns of PDK mRNA expression in an attempt to elucidate the importance of PDH regulation in hypoxia survival in the hypoxia-tolerant killifish.
| METHODS AND MATERIALS |
|---|
|
|
|---|
Adult killifish (F. heteroclitus heteroclitus L.; weighing 3.1 ± 0.13 g wet weight; mean ± SE) were captured from New Brunswick, Canada and transported to the University of British Columbia and housed indoors in aquaria with filtered synthetic brackish water (10 g/l; Deep Ocean, Carson, CA) made up in dechlorinated City of Vancouver tap water. Water temperature was maintained at 20°C, and oxygen partial pressure was
160 Torr. Fish were held under these conditions for at least 30 days before experimentation. Fish were fed commercial trout chow (PMI Nutrition International, Brentwood, MO) at
2% body wt/day. All experimental procedures were conducted according to the guidelines set out by the Canadian Council for Animal Care administered by the University of British Columbia Animal Care Committee and were approved by the University of British Columbia Animal Care Committee.
Gene Identification and Sequencing
Tissue sampling and RNA extraction. Adult killifish were taken directly from their holding tanks and killed by concussion; samples of brain, gill, intestine, liver, muscle, spleen, and heart were immediately dissected and frozen in liquid nitrogen. Samples were stored at –80°C until RNA extraction. Total RNA was extracted from tissues using the guanidinium thiocyanate method outlined by Chomczynski and Sacchi (9) using TriPure Isolation Reagent (Roche Diagnostics, Montreal, QC, Canada). Total RNA was quantified spectrophotometrically, and the quality was assessed using agarose gel electrophoresis to verify the presence of two ribosomal RNA bands.
cDNA was synthesized from 5 µg of total RNA isolated from the above tissues using oligo(dT18) primer and RevertAid H Minus M-MuLV reverse transcriptase following the manufacturers instructions (MBI Fermentas, Burlington, ON, Canada). Partial PDK isoform sequences were obtained using degenerate primers (PDK forward and reverse; Table 1) determined from conserved regions of available sequences in GenBank. Primers were designed with the assistance of GeneTool Lite software (version 1.0; DoubleTwist, Oakland, CA; http://www.biotool.com). PCRs were carried out in a PTC-200 MJ Research thermocycler using Taq DNA polymerase (MBI Fermentas), and cDNA was isolated from the above tissues. Each PCR reaction consisted of an initial treatment of 5 min at 94°C, followed by 40 cycles of 30 s at 94°C, 30 s at 49°C, and 90 s at 72°C. The resulting PCR products were electrophoresed along with a DNA ladder on 1.5% agarose gels containing ethidium bromide to confirm the presence of a product of the appropriate size (
900 bp). The PCR reactions yielded a single visible band during electrophoresis; therefore, the PCR product was directly ligated into pGEM-T easy vectors (Promega; Fisher Scientific, Nepean, ON, Canada) without gel purification. The ligated product was introduced into heat shock competent Escherichia coli (strain 109; Promega) and colonies grown on Luria Bertani agar plates containing 100 ng/ml ampicillin, 0.1 mM isopropyl β-D-1-thiogalactopyranoside, and 20 ng/ml 5-bromo-4-chloro-3-indolyl-b-D-galactoside (X-Cal). Several colonies containing ligated PCR products were selected and grown overnight in liquid culture. Plasmids were harvested from liquid culture using GenElute Plasmid Miniprep kit (Sigma-Aldrich, Oakville, ON, Canada) and sequenced using an ABI377 automated fluorescent sequencer (Applied Biosystems).
|
Phylogenetic Analysis
Amino acid sequences were deduced from the nucleotide sequence of each isoform using GeneTool Lite software. Protein sequences or deduced amino acid sequences of other PDK genes were obtained from GenBank or genome sequences available on Ensembl: human2 (Homo sapiens, NP002602); dog2 (Canis familiaris, XP548195); mouse2 (Mus musculus, AAH21764 [GenBank] ); rat2 (Rattus norvegicus, NP110499); opossum2 (Monodelphis domestica, XP001374336); chimpanzee2 (Pan troglodytes, XP001168889); frog2 (Xenopus leavis, NP001080097); frog4 (NP001006803); zebrafish3 (Danio rerio, XP00133259); zebrafish2 (NP957290); pufferfish (Tetraodon nigroviridis, CAF91901 [GenBank] ); rat1 (NP446278); chicken1 (Gallus gallus, NP001026523); mouse1 (NP766253); human1 (NP002601); cow (Bos Taurus, XP607085); chicken3 (NP001006259); mouse3 (NP663605); opossum3 (XP001362369); horse (Equus caballus, XP001493781); mouse4 (NP038771); dog3 (XP859771); human3 (NP005382); rat4 (NP446003); dog4 (XP539427); human4 (NP002603); chimpanzee 4 (XP527822); monkey4 (Macaca mulatta, XP001093255); platypus (Ornithorhynchus anatinus, XP001514685); dog1 (XP534032); horse (XP001495143); opossum1 (XP001377950); chimpanzee1 (XP001146941); cow (XP583960); monkey1 (XP001086316); bat4 (Rhinolophus ferrumequinum, DQ468386 [GenBank] ); squirrel4 (Spermophilus tridecemlineatus, AF020845 [GenBank] ); bee (Apis mellifera, XP393904); Stickleback2a (Gasterosteus aculeatus, Ensemble Gene ID ENSGACG00000003817); Stickleback2 (Ensemble Gene ID ENSGACG00000005722); stickleback4 (Ensemble Gene ID ENSGACG00000006683); stickleback1 (Ensemble Gene ID ENSGACG00000004843); stickleback3 (Ensemble Gene ID ENSGACG00000002603); fugu4 (Takifugu rubripes, Ensemble Gene ID ENSTRUG00000015972); fugu2 (Ensemble Gene ID ENSTRUG00000013086).
Phylogenetic analysis was performed using the neighbor-joining methods (41) with complete deletion of gaps using MEGA4 software (51). The support for each node was assessed using 500 bootstrap replicates. Isoforms were named according to their position in the phylogenetic tree.
PDK Isoform Tissues Distribution
Tissue distribution of the two identified killifish PDK isoforms was estimated using quantative PCR (qPCR) with isoform-specific primers. Primers were designed with the assistance of Primer Express software (primer sequences in Table 1; version 2.0.0; Applied Biosystems, Foster City, CA). The same cDNA samples used above for gene identification were used to determine the tissue distribution of each mRNA. qPCR was performed on an ABI 7000 sequence analyzer (Applied Biosytems) using the protocols outlined in Richards et al. (39). The absolute mRNA levels of each PDK isoform among tissues was estimated according to the formula, efficiency–Ct, where efficiency is the slope of a standard curve and Ct corresponds to the threshold cycle number. All results were standardized to total RNA and normalized to the levels of PDKFh-2 mRNA detected in the brain.
Experimental Protocols
Whole animal respirometry (Pcrit and
O2).
Routine oxygen consumption rate (
O2) and critical oxygen tension PO2 (Pcrit) were determined on individual fish by using closed respirometry. Briefly, individual fish were placed into 250-ml darkened glass respirometers and allowed to habituate for at least 12 h before respirometry. During this period, air-saturated brackish water (10 g/l) was circulated through the respirometer. At the beginning of a trial, water circulation was stopped, an oxygen probe was fixed into each respirometer, and the respirometer was sealed. Water temperature was maintained at 20°C by placing the respirometers in a temperature-regulated wet table. Continuous oxygen readings were obtained using a fiber optic oxygen probe (Foxy system; Ocean Optics, Dunedin, FL), while the water in the respirometer was continuously stirred in one corner with a magnetic stir bar. Oxygen partial pressure was monitored continuously until it decreased < 20 Torr or the fish showed signs of distress. Each trial took
1.5 h to complete, and the fish typically remained quiescent throughout. Any fish that struggled was removed from the data set. Water oxygen content was corrected for differential solubility in brackish water at 19°C by using solubility coefficients from Boutilier et al. (5) and
O2 was calculated over sequential 10-min periods using the slope of the oxygen trace (water O2 over time) corrected for fish weight and respirometer volume. Changes in
O2 were plotted against the mean water PO2 for each 10-min interval, and the inflection point at which
O2 transitioned from being dependent on environmental PO2, termed Pcrit (36), was determined using the BASIC program designed by Yeager and Ultsch (58).
To understand the temporal effects of hypoxia exposure on
O2, we measured
O2 on a separate group of fish over a 30-min period in normoxic water (124 Torr) and after 3-h exposure to 10 Torr. Briefly, fish were introduced to the same respirometers used for the Pcrit determination and were allowed to habituate overnight. To record resting
O2 under normoxic conditions, we stopped water flow to the respirometer and monitored water PO2 over
30 min or to the point at which water PO2 had decreased by 20%, whichever came first. At this point, water flow through the respirometer was reestablished, and the head tank supplying each respirometer was bubbled with nitrogen gas to achieve the desired level of hypoxia (10 Torr). Fish were maintained under these hypoxic flow-through conditions for 2 h 45 min at which point the respirometers were sealed, and
O2 was determined over a 30-min period.
Hypoxic Exposure for Metabolic Analysis
One day before experimentation, four fish were transferred into each of 14, 4-l plastic exposure chambers that had a mesh-covered top and mesh sides with a 1.2-l reservoir when the exposure chamber was emerged. The exposure chambers containing fish were placed into a 550-l glass aquaria with aeration, and filtered brackish water and maintained at 19°C. Water was circulated from the bulk tank into each chamber via a submersible pump to ensure a uniform oxygen tension throughout the exposure system. Fish were not fed during the 1-day transfer recovery period or during experimentation.
Following the 1-day recovery period, blood, white muscle, and liver were sampled from fish exposed to the normoxic water (127 ± 6 Torr; n = 8). To sample fish, the exposure chamber was gently removed from the aquaria, confining the fish to the bottom 1.2-l reservoir of the chamber. Fish were then anesthetized by injecting a concentrated solution of benzocaine dissolved in 95% ethanol (Sigma-Aldrich; final concentration 1 g/l) into the water. At complete anesthesia (<45 s) fish were removed from the water, and a section of trunk musculature was taken posterior to the visceral cavity and immediately freeze-clamped between two aluminum blocks cooled in liquid nitrogen. The liver was also dissected free of the animal and frozen in liquid nitrogen. Tissue sampling from all four fish took <2 min and two exposure chambers were sampled for a total of eight fish per sampling time.
Following normoxia sampling, hypoxia was induced by bubbling nitrogen into the large tank of water and covering the surface with plastic. Water PO2 decreased from normoxia to 23 ± 1 Torr over the first 1 h of nitrogen bubbling and was then maintained at 6 ± 1 Torr from 1.5 to 15 h. Water oxygen was monitored inside the exposure chambers using an YSI-85 dissolved oxygen meter. Eight fish were sampled at 3-, 6-, 9-, 12-, and 15-h exposure to hypoxia, as described above. Following the sampling at 15 h, nitrogen gas was replaced with air, and the fish were allowed to recover for 12 h in normoxia (107 ± 5.4 Torr; n = 6) and were then sampled.
Analytical Techniques
Frozen muscle was broken into pieces (50–100 mg) using an insulated mortar and pestle cooled in liquid nitrogen. Several aliquots of muscle were stored separately in liquid nitrogen for determination of the activity of the PDHa as previously described by Richards et al. (38). A portion of the frozen muscle (
50 mg) was ground into a fine powder under liquid nitrogen and used for intracellular pH (pHi) measurements as described by Pörtner et al. (35). The remaining muscle was lyophilized for 72 h, dissected free of connective tissue, powdered, and stored dry at –80°C for subsequent analysis.
For the determination of muscle glycogen,
20 mg of lyophilized muscle was digested in 1 ml 30% KOH at 100°C. Glycogen was isolated as described by Hassid and Abraham (17), and free glucose was determined after digestion with amyloglucosidase (2). For the extraction of metabolites from white muscle, aliquots of lyophilized muscle (
20 mg) were weighed into borosilicate tubes, 1 ml of ice-cold 1 M HClO4 was added, and then homogenized at the highest speed of a Syngene homogenizer for 20 s at 0°C. Homogenates were transferred to 2-ml microcentrifuge tubes, centrifuged at 4°C for 5 min at 20,000 g, and the supernatant was neutralized with 3 M K2CO3, except for samples for ammonia analysis, which were neutralized with 3M Tris-base. Neutralized extracts were assayed spectrophotometrically for ATP, CrP, creatine, lactate, pyruvate, glutamate,
-ketoglutarate, and glucose using the spectrophotometric methods described in Bergmeyer (2). Muscle ammonia was analyzed using a commercial kit (Sigma). Muscle acetyl-CoA was determined on neutralized extracts by radiometric methods previously described in Richards et al. (38).
The expression patterns of killifish PDK mRNA were estimated using qPCR and the isoform-specific primers used above and given in Table 1. The expression pattern of lactate dehydrogenase LDH-A (muscle) and LDH-B (liver) mRNA were examined using qRT-PCR primers (Table 1) designed from LDH-A and -B sequences from GenBank (LDH-A accession no. L43525; LDH-B accession no. M33969). Briefly, total RNA was extracted from
20 mg of muscle or liver and 5 µg of total RNA was reverse transcribed into cDNA using the protocols outlined above. qPCR was performed as described in Richards et al. (39). qPCR reactions contained 2 µl of cDNA, 4 pmol of each primer, and Universal SYBR Green master mix (Applied Biosystems) in a total volume of 20 µl. All qPCR reactions were performed as follows: 2 min at 50°C, 10 min at 95°C, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Melt curve analysis was performed following each reaction to confirm the presence of only a single product of the reaction. To control for genomic contamination, negative control samples were performed using RNA that had not be reversed transcribed. Genomic DNA contamination was present for all genes examined, but never constituted more than 1:1024 copies (<0.1%) for any gene. One randomly selected control sample was used to develop a standard curve relating threshold cycle to cDNA amount for each primer set. All results were expressed relative to these standard curves, and mRNA amounts were normalized to the expression of elongation factor (EF-1
), which was measured using primers described by Scott et al. (42). There was no effect of hypoxia exposure on the expression of EF-1
when expressed as a function of total RNA (data not shown); thus EF-1
appears to be a good control gene for hypoxia studies. Expression levels in hypoxia-exposed animals were expressed relative to the mean expression levels in the normoxia control samples. All samples were run in duplicate and the coefficient of variation between duplicate samples was always <15%.
PDH kinase 2 protein content was assessed using immunoblot analysis according to the protocols described in Richards et al. (40). Briefly, muscle samples (
20 mg) were homogenized in a buffer containing 100 mM TRIS·HCl, 1% SDS, 5 mM EDTA, 1 µg/ml aprotinin, 1 µg/ml pepstatin A, 1 µg/ml leupeptin, and 20 µg/ml PMSF, pH 7.5. Homogenates were centrifuged at 5,000 g for 10 min at 4°C, the supernatant was assayed for total protein using the methods of Bradford (7), and a portion of the supernatant was denatured by boiling for 3 min in SDS-sample buffer (26). Denaturing SDS-polyacrylamide gels were loaded with denatured muscle homogenates at a protein concentration of 20 µg protein/lane, and electrophoresed for 15 min at 75 V followed by 75 min at 150 V. An identical control sample was included on each gel to control for gel-to-gel variation. Following the electrophoresis, proteins were transferred to nitrocellulose membranes (Bio-Rad Laboratories) using a Trans-Blot semidry transfer cell (Bio-Rad). Blots were incubated for 1 h with 1 µg/ml of primary antibody (raised against the COOH-terminal end of rabbit PDK-2; peptide sequence 5'-VPSTEPKNTSTYRVS-'3; Abgent, San Diego, CA) and then subsequently incubated for 1 h with alkaline phosphatase-conjugated goat-anti-rabbit IgG secondary antibody (Stressgen, Ann Arbor, MI). Immunoblots were developed in alkaline phosphatase buffer containing 5-bromo-4-chloro-3-indolylphosphate and nitroblue tetrazolium, and the band intensity was quantified using a FluorChem 8800 image analyzer (Alpha Innotech, San Leandro, CA) and AlphaEase FC software (version 3.1.2). Individual samples were expressed relative to total protein, and subsequently each sample was normalized to the mean of the normoxia control samples.
Calculations and Statistical Analysis
Free cytosolic [ADP] and [AMP] were calculated from measured values of [ATP], [CrP], [creatine] ([Cr]), and pHi according to published protocols (13, 52), assuming equilibrium of the creatine kinase (K'CK) and adenylate kinase (K'AK) reactions. The equilibrium constants for K'CK and K'AK were corrected for experimental temperature, pH, and free Mg2+ (assumed to be 1 mM; see Refs. 13 and 52). Cellular phosphorylation potential was calculated according to the following equation
![]() |
Cytosolic free [Pi] was estimated assuming the inverse of CrP hydrolysis and starting with a resting tissue level of 1 µmol/g wet weight (16). The Gibbs free energy of ATP hydrolysis (
fG'; kJ/mol) was determined using the following equation
![]() |
where R is the universal gas constant (8.3145 J·K–1·mol–1), T is temperature in K, and
fGATP'0 is the standard transformed Gibbs energy of ATP hydrolysis (
fGATP'0 = –RTlnK'ATP) at measured pH and temperature and estimated free [Mg2+]. Cytoplasmic [NAD+]/[NADH] was estimated using the apparent equilibrium of the lactate dehydrogenase reaction [Keq from Wang et al., (55)] and measured pHi, [lactate], and [pyruvate]. Mitochondrial [NAD+]/[NADH] was estimated from the glutamate dehydrogenase (GDH) reaction using whole cell measurements of [NH3], [glutamate], [
-ketoglutarate], Keq from Williamson et al. (57) and estimated mitochondrial pH from Moyes et al. (31). All data are presented as means ± SE (n). All metabolite concentrations determined on lyophilized tissues were converted back to wet weights by taking into account a wet-to-dry ratio of 4:1 (55). Statistical analysis consisted of a one-way ANOVA followed by the Holm-Sidak post hoc method of pairwise multiple comparisons. All data were tested for normality (Kolmogorov-Smirnov test) and homogeneity of variance (Levene median test), and in cases where data sets did not meet these assumptions, data were log transformed and statistical analyses were repeated. For a small number of cases where log transformation did not yield data of normal distributed or equal variance, Kruskal Wallis one-way ANOVA on ranks was used for statistical analysis. Results were considered significant at P < 0.05.
| RESULTS |
|---|
|
|
|---|
Two full-length cDNA sequences coding for PDK were identified in killifish. Alignment of these killifish isoforms with other vertebrate PDK isoforms indicated that our killifish isoforms were 52–75% similar to other known fish and mammalian PDK isoforms. Phylogenetic analysis of the deduced amino acid sequence of these two killifish isoforms indicated that one isoform grouped more closely with other known PDK-4 isoforms (designated PDKFh-4), and the other isoform grouped closely with other known PDK-2 isoforms (designated PDKFh-2; Fig. 1). No isoforms matching PDK-1 or -3 were detected. PDKFh-2 and PDKFh-4 cDNA sequences have coding regions of 1224 and 1230 nucleotides, yielding an open reading frame of 409 and 410 amino acids, respectively.
|
(data not shown), the tissue distribution patterns of PDKFh-2 and -4 were qualitatively similar to that observed when normalized to total RNA.
|
Oxygen consumption rate and Pcrit.
During progressive decreases in environmental PO2,
O2 remained constant from water PO2 levels of 153 ± 3 Torr to 66 ± 1 Torr, below which
O2 decreased in a near-linear fashion with environmental PO2 (Fig. 3A). The point at which killifish transition their
O2 from being independent of environmental PO2 to being dependent on environmental PO2 was 63.9 ± 5.6 Torr (n = 6). Killifish exposed to hypoxia for 3 h (water PO2 10 ± 4 Torr) and tested under hypoxia had approximately sevenfold lower
O2 than killifish exposed to normoxic water (water PO2 124 ± 6 Torr; Fig. 3B).
|
Water oxygen decreased over the first 1.5-h exposure from normoxia (127 ± 6 Torr; n = 8) to hypoxia and then stabilized at 6 ± 1 Torr (n = 131) for the 15-h duration of the hypoxia exposure (Fig. 4A). During the decline in environmental oxygen, there was an initial period where the fish appeared agitated and were observed swimming around the top of each chamber, but by 2 h after the initiation of hypoxia, all fish were quiescent at the bottom of their individual chambers and remained in this state throughout the hypoxia exposure. Killifish that experienced 15-h exposure to hypoxia followed by 12-h recovery in normoxic water remained quiescent throughout the hypoxia and recovery period. No mortality was observed as a result of the hypoxia exposure or recovery.
|
White muscle [ATP] was not affected by exposure to hypoxia or hypoxia followed by 12-h recovery in normoxic water (recovery group; Fig. 4B). Exposure to hypoxia caused an initial 13% decrease in white muscle [CrP] over the first 3 h, which continued to decrease to 11% of the normoxic [CrP] values at 15-h exposure to hypoxia (Fig. 4C). White muscle [CrP] in the recovery group was similar to that observed in the initial normoxic control. Changes in white muscle free [Cr] were reciprocal and roughly stiochiometric to the changes observed in white muscle [CrP] (Table 2). White muscle pHi decreased within the first 3-h exposure to hypoxia and continued to decrease until 9 h, where white muscle pHi stabilized at
6.8 until 15 h (Table 2). pHi was similar between the initial normoxic controls and in the recovery group. Estimated [ADPfree] and [AMPfree] both increased in response to hypoxia exposure, although [ADPfree] was elevated by 6-h exposure to hypoxia, whereas [AMPfree] was only elevated after 15-h exposure to hypoxia (Table 2). Muscle [AMPfree] returned to normoxic levels in the recovery group, while [ADPfree] remained elevated compared with the normoxic controls. As a result of these changes in free adenylates, cellular energy charge ([ATP]/[ADPfree]) and [ATP]/[AMPfree] were significantly depressed after 6- and 3-h exposure to hypoxia, respectively, and remained low for the full 15-h exposure to hypoxia (Fig. 4D and Table 2). In the recovery group, both [ATP]/[ADPfree] and [ATP]/[AMPfree] were similar to the normoxic control values. Muscle phosphorylation potential was significantly depressed by 3-h exposure to hypoxia and continued to decrease throughout the hypoxia exposure returning to values that were not significantly different from the normoxic controls in the recovery group (Table 2). The
fG' decreased significantly by 6-h exposure to hypoxia and continued to decline until 15-h hypoxia exposure (Table 2). In the recovery group,
fG' was not different from that observed in the normoxic controls.
|
|
|
-ketoglutarate (Table 4). Estimated mitochondrial [NAD+]/[NADH] did not change significantly in response to hypoxia exposure or recovery, but there was an overall trend for a higher [NAD+]/[NADH] during hypoxia exposure compared with normoxic controls or recovery.
|
50% within 3-h exposure to hypoxia and remained at this low level for the full 15-h exposure, returning to values that were not significantly different from those in the recovery group (Fig. 7A). Total PDK-2 protein showed a general decrease of nearly 40%, between 6 and 15 h of hypoxia exposure (Fig. 7B). A full recovery of PDK-2 protein was observed in the recovery group. PDKFh-2 mRNA in white muscle increased
3.8-fold over the first 9-h exposure to hypoxia and continued to increase to 10.2 times the normoxic values at 15 h (Fig. 7C). In the recovery group, the level of PDKFh-2 expression was similar to the initial normoxic values. Exposure to hypoxia had no effect on PDKFh-4 mRNA expression in white muscle (Fig. 7D). White muscle LDH-A mRNA expression did not change in response to hypoxia exposure to recovery (Table 3). Liver
Hypoxia exposure did not affect the expression pattern of LDH-B or PDKFh-4 in liver (Table 5). In contrast, liver PDKFh-2 expression transiently decreased, by
38%, over the first 3-h exposure to hypoxia, returned to values that were not significantly different from the controls between 6 and 12 h, and significantly decreased again at 15-h exposure (Table 5). PDKFh-2 expression in the recovery group was similar to that observed in the initial normoxic control.
|
| DISCUSSION |
|---|
|
|
|---|
O2. In killifish muscle, these responses occurred in a biphasic fashion and showed a strong dependence on muscle CrP hydrolysis (Fig. 4C) and glycolytic production of lactate (Fig. 5B) to support ATP turnover during the first 6-h exposure to hypoxia. This was followed by a period of substantially reduced utilization of CrP and lower lactate accumulation at > 9-h exposure. Fish white muscle is generally considered a "closed" system, with little exchange of metabolites, in particular, lactate, with the extracellular fluid (54); therefore, these data indicate that during the first 6 h of exposure to hypoxia, increases in substrate-level phosphorylation are necessary to make up for the reduced capacity for ATP generation. After the first 6 h of exposure to hypoxia, the smaller changes in muscle metabolite concentrations coincide with stable muscle [ATP] (Fig. 4B) and suggest dramatic reductions in ATP consumption and tissue metabolic rate depression.
Although maintenance of a constant cellular [ATP] during hypoxia exposure (Fig. 4B) has been termed the "hallmark" of a hypoxia-tolerant animal, there was still a substantial disruption of cellular energetics and a loss of phosphorylation potential (Table 2) during hypoxia exposure. The loss of phosphorylation potential would have dramatic effects on the rates of ATP utilization as well as rates of cellular ATP production and substrate oxidation. In particular, hypoxia exposure was associated with a significant drop in
fG' from –66.8 kJ/mol to –50.9 kJ/mol over 15 h (Table 2). Estimates of the critical limit of
fG' for the maintenance of cellular function (16, 21) suggest that below a threshold of –52 kJ/mol, cellular processes such as ion pumping can no longer derive sufficient energy from ATP hydrolysis to be maintained. Thus, although total tissue [ATP] is not affected by hypoxia exposure, significant disruption in cellular energy balance is still observed in response to hypoxia, and this disruption affects rates of cellular ATP use and production.
The control of metabolic rate depression is thought to be primarily regulated by the kinetics of mitochondrial substrate oxidation (4). Consistent with this notion is the rapid (3-h) and sustained reduction in muscle PDHa activity during hypoxia exposure (Fig. 6A), which limits the rate of mitochondrial pyruvate oxidation, facilitating pyruvate reduction to lactate via cytoplasmic LDH (Fig. 5B). The conversion of pyruvate to lactate and the large accumulation of muscle lactate (Fig. 5B) during hypoxia exposure is meant to aid in the maintenance of cytoplasmic [NAD+]/[NADH] for sustained glycolytic flux; however, our estimates of cytoplasmic [NAD+]/[NADH] indicate that the cytoplasm rapidly becomes reduced during hypoxia exposure (Fig. 5C).
|
40 times more PDKFh-2 mRNA transcripts than PDKFh-4, suggesting that PDKFh-2 is the primary isoform present and therefore regulated to a greater degree by the increase in [acetyl-CoA] (Table 3).
Considerable controversy surrounds the use of the LDH and GDH equilibrium for estimating cytoplasmic and mitochondrial [NAD+]/[NADH] in tissues (e.g., see Ref. 22). For example, although GDH is restricted to the mitochondrial matrix and is potentially useful for estimating mitochondrial [NAD+]/[NADH], the total activity of GDH is low, in part due to a low mitochondrial volume density in muscle; and there is uncertainty to what extent cytoplasmic levels of NH3, glutamate and
-ketoglutarate influence the enzyme equilibrium. The usefulness of the LDH reaction for estimating cytoplasmic [NAD+]/[NADH] has also been drawn into question because there is evidence that LDH is not isolated to the cytoplasm but is also located in the mitochondrial fraction of fish white muscle (30). It is unlikely, however, that the presence of mitochondrial LDH influences our estimates of cytoplasmic [NAD+]/[NADH] to a great degree because of the low mitochondrial content characteristic of fish white muscle. Overall, our estimates of cytoplasmic and mitochondrial [NAD+]/[NADH] in killifish muscle must be considered tentative and await confirmation by other methods.
We observed a rapid inhibition of PDH in killifish muscle during hypoxia exposure that was consistent with work in mammalian ischemic/hypoxia cell models. These studies suggest that modulation of PDHa activity during hypoxia exposure is essential to cell survival leading to a stabilization of cellular [ATP] (43), an attenuation of hypoxia-related ROS generation (24), and inhibition of mitochondrial-induced apoptotic pathways (e.g., see Ref. 12). Important to the inhibition of PDH in mammalian cell culture systems during hypoxia (24) is the HIF-1-mediated increases in the expression of PDK-1 (23, 24). Four PDK isoforms have been identified in mammals (Fig. 1), but our gene identification approach using degenerate primes (Table 1) identified only two isoforms that loosely group with other known PDK-2 and -4 isoforms, each with a specific tissue distribution (Fig. 2). Fish often possess more isoforms than mammals because of genome duplication events (e.g., stickleback have two PDK-2 isoforms; Fig. 1), but our degenerate primers and PCR protocols only amplified a single representative of each of PDK-2 and -4. Of the isoforms identified in killifish, only PDKFh-2 was shown to be hypoxia responsive in muscle with mRNA levels increasing > 10-fold during the 15-h hypoxia exposure (Fig. 6C). The qRT-PCR primers used to amplify PDKFh-2 were designed in a highly variable portion of the gene and therefore can be considered highly gene specific and unlikely to yield amplification of other unidentified PDK-2 isoforms that may be present in the killifish genome. The hypoxia-induced increases in PDKFh-2 expression in killifish muscle suggest that this gene could be regulated by HIF via a hypoxic responsive element in the promoter region of PDKFh-2. However, the expression of other genes that are known to have associated hypoxic responsive element (37, 44), including muscle LDH-A (Table 3) and liver LDH-B (Table 5), were not affected by the hypoxia exposure in killifish, thus bringing into question whether other signal transduction pathways may be involved in regulating the PDKFh-2 hypoxic response. Other signal transduction cascades (e.g., MAPK; 27) are known to be activated during hypoxia and could contribute to the regulation of PDKFh-2 expression during hypoxia exposure in killifish muscle. Additional work characterizing the transcription factor binding sites in the promoter regions of PDKFh-2 is warranted. No changes in muscle PDKFh-4 mRNA were observed (Fig. 6D), but the heart-specific expression pattern of PDKFh-4 (Fig. 2) suggests that this isoform may only be important in heart tissues.
In spite of the hypoxia sensitivity of PDKFh-2 expression in muscle (Fig. 6C), PDKFh-2 expression was not hypoxia sensitive in the liver of killifish. In fact, the expression of PDKFh-2 significantly decreased within the first 3-h exposure to hypoxia and remained somewhat depressed compared with both the normoxia and recovery sample (Table 3). Similarly, the expression of both PDKFh-4 and LDH-B in the liver was not affected by the hypoxia exposure (Table 3). Clearly, patterns of gene expression in response to hypoxia exposure are specific to the tissues examined, which is consistent with the tissue-specific enzyme responses observed in the gulf killifish (F. grandis) acclimated to moderate hypoxia (
3.6 kPa or 27 Torr) for 1 mo (29). Martinez et al. (29) observed significant increases in the activities of five out of the eleven glycolytic enzymes that were analyzed in the liver (including LDH), which suggests a strong activation of substrate-level phosphorylation and oxygen-independent ATP production. In the same study, the results observed in liver differed dramatically from those observed in muscle where there were clear decreases in glycolytic enzyme activities, suggesting, unlike the liver, the muscle undergoes a substantial metabolic rate depression. Combined, these results suggest that the liver of killifish may make up for the hypoxia-imposed energy deficit through a strong activation of glycolysis, whereas in the muscle, strong inhibition of mitochondrial oxygen consumption and overall tissue metabolic rate depression may play a more important role in insuring hypoxia survival.
There is a lack of association between the hypoxia-induced changes in muscle PDKFh-2 mRNA expression and PDK-2 protein content (c.f. Fig. 6, B and C), such that the changes in mRNA were not translated into increased protein content during hypoxia exposure. Why then are there changes in mRNA transcription during hypoxia exposure without changes in total protein content? There are several possible explanations. First, the use of mammalian antibodies may not be specific for the PDK-2 isoforms identified in fish and analyzed using qPCR. Second, although the time course for hypoxia exposures in these various studies were chosen to replicate natural durations of hypoxia exposure, it is possible that these hypoxia exposures were too short for translation of the unregulated mRNAs to occur. Thus, it is possible that had the hypoxia exposure period been longer we might have observed increased protein content. Third, increases in mRNA transcript levels for select proteins may be required to support targeting of these transcripts to the limited capacity translation machinery (45) to maintain total protein levels during hypoxia. This may be increasingly important for nuclear-encoded mitochondrial proteins where hypoxia may affect rates of mitochondrial protein degradation and import. Fourth, and possibly most intriguing, is the notion that the high levels of mRNA transcripts accumulate during hypoxia to facilitate rapid turnover of select proteins during the early part of recovery from hypoxia. Although, Lewis et al. (28) demonstrated that there was no hyperactivation of protein synthesis during recovery from hypoxia exposure in the Amazonian cichlid (Astronatus ocellatus), increased levels of mRNA transcripts could mean preferential synthesis of metabolically important proteins without the need to increase total protein synthesis capacity. Our experimental design included a recovery group, but the only sample taken was 12 h after hypoxia exposure ended, and therefore any rapid protein turnover would be completed. Indeed, almost all parameters measured in this recovery group, except cytoplasmic [NAD+]/[NADH] and [ADPfree], had returned to control normoxic levels. Further research into the turnover rates of metabolically important proteins during hypoxia exposure and recovery is necessary to understand the counterintuitive increase in mRNA expression without changing total protein.
Killifish are well known for their hypoxia tolerance (34, 48) and have an amazing capacity to activate oxygen-independent means of ATP production to survive. It is surprising however, that killifish have relatively high Pcrit values (63.9 ± 5.6 Torr) and are oxygen conformers across a very large range of environmental oxygen tensions. The high Pcrit value is also surprising considering their very low whole red cell Hb-P50 of 3.8 to 5.0 Torr (11). Our estimates of Pcrit differ substantially from those measured in a closely related species, the gulf killifish (F. grandis), which had an estimated Pcrit of
34 Torr. The disparity in Pcrit values could be due to species differences or differences in their physiological condition. Henriksson et al. (18) demonstrated in two species of sculpin from the family Cottidae (Pacific staghorn sculpins, Leptocottus armatus and Prickly sculpin, Cottus asper) that measures of Pcrit are subject to change from any parameter that affects whole animal oxygen uptake (e.g., gill morphology). It is possible that differences in gill morphology or modifications to other aspects of the oxygen cascade in our killifish could impact their Pcrit, but despite this high Pcrit value, our fish were able to survive 15-h exposure to an enviromental O2 tension of
10% of their Pcrit value. At this level of hypoxia, an animal can either attempt to sustain metabolic rate through the activation of substrate level phosphorylation or reduce metabolic rate via metabolic rate depression, and it appears that killifish employ both strategies.
Perspectives and Significance
Hypoxia survival involves a dramatic reorganization of metabolic processes to reduce ATP demands to match the limited capacity for oxygen-independent ATP production. The role of PDH in hypoxic survival has been examined extensively in carcinoma cell lines and ischemic mammalian models; however, the present study is the first to examine the complex regulation of PDH in an organism that has evolved strategies for hypoxia survival. Inactivation of PDH in killifish appears central to hypoxic survival and reaffirms the importance of limiting mitochondrial pyruvate oxidation during hypoxia to limit mitochondrial respiration and prevent ROS production and mitochondrial-mediated initiation of apoptosis. The parallels observed between killifish PDH regulation during hypoxia exposure and the responses observed in many carcinoma cell lines opens the possibility of using organisms that have evolved to survive hypoxia exposure to understand the basic biology of cancer and other hypoxia related heath issues.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
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 |
|---|
|
|
|---|
from the estuarine teleost, Fundulus heteroclitus: interaction of HIF-2
with two ARNT2 splice variants. J Exp Zool 294: 17–29, 2002.[CrossRef][Web of Science][Medline]
-isoform switching in gills of rainbow trout (Oncorhynchus mykiss) during salinity transfer. J Exp Biol 206: 4475–4486, 2003.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |