Overexpression and Knockdown of Hypoxia-Inducible Factor 1 Disrupt the Expression of Steroidogenic Enzyme Genes and Early Embryonic Development in Zebrafish.

Hypoxia is an important environmental stressor leading to endocrine disruption and reproductive impairment in fish. Although the hypoxia-inducible factor 1 (HIF-1) is known to regulate the transcription of various genes mediating oxygen homeostasis, its role in modulating steroidogenesis-related gene expression remains poorly understood. In this study, the regulatory effect of HIF-1 on the expression of 9 steroidogenic enzyme genes was investigated in zebrafish embryos using a "gain-of-function and loss-of-function" approach. Eight of the genes, CYP11a, CYP11b2, 3β-HSD, HMGCR, CYP17a1, 17β-HSD2, CYP19a, and CYP19b, were found to be differentially upregulated at 24 and 48 hpf following zHIF-1α-ΔODD overexpression (a mutant zebrafish HIF-1α protein with proline-414 and proline-557 deleted). Knockdown of zHIF-1α also affected the expression pattern of the steroidogenic enzyme genes. Overexpression of zHIF-1α and hypoxia exposure resulted in downregulated StAR expression but upregulated CYP11a and 3β-HSD expression in zebrafish embryos. Conversely, the expression patterns of these 3 genes were reversed in embryos in which zHIF-1α was knocked down under normoxia, suggesting that these 3 genes are regulated by HIF-1. Overall, the findings from this study indicate that HIF-1-mediated mechanisms are likely involved in the regulation of specific steroidogenic genes.

Introduction

Hypoxia caused by aquatic ecosystem enrichment with nutrients and organic matter is a major global threat that is predicted to worsen as climate change progresses. The adverse impacts of aquatic hypoxia include significant reductions in fisheries production, fish growth, alteration of species composition, and mass fish mortality.1,2 Hypoxia impairs fish reproduction by inhibiting testicular and ovarian development, reducing sperm and egg quality, affecting fertilization, and hatching and influencing the survival of larvae and juvenile fitness.3 Ovarian and testicular growth, gametogenesis, and oocyte and sperm production were significantly reduced in Atlantic croakers in hypoxic regions of the Mississippi River4,5 and in the northern Gulf of Mexico.6

In fish, the inhibition of key reproductive processes by hypoxia is associated not only with reduced metabolism but also with the repression of specific hormones and hormone receptors along the hypothalamus-pituitary-gonad (HPG) axis.3,7 Numerous in vivo studies have shown that fish reproduction is affected by hypoxia via alteration of the steroidogenesis pathway.7–11 Steroid hormones biosynthesis is a fundamental process for reproduction in vertebrates, and all steroid hormones are derived from cholesterol. The rate-limiting step in steroid production is the delivery of cholesterol from the outer to the inner mitochondrial membrane, which is mediated first by the steroidogenic acute regulatory (StAR) protein12 followed by the conversion of cholesterol to pregnenolone by the cytochrome P450 cholesterol side-chain cleavage (P450scc) enzyme encoded by CYP11a1.13 Pregnenolone is then metabolized to progesterone by 3β-hydroxysteroid dehydrogenase/Δ5-Δ4 isomerase (encoded by 3β-HSD)14 or hydroxylated to 17α-hydroxyprogesterone by cytochrome P450 17α-hydroxylase (encoded by CYP17a1). Progesterone and 17α-hydroxyprogesterone are then converted through a series of steps to androgen and estrogen by enzymes such as 17β-hydroxysteroid dehydrogenase (17β-HSD) and P450 aromatase (CYP19).15 Estradiol and testosterone are steroid hormones that are important for sexual and reproductive development in animals.16,17 Molecular oxygen is essential for the synthesis of sex steroid hormones,18 and recent studies revealed that the expression of certain steroidogenic enzyme genes, such as StAR, CYP19, and CYP17A1 (which specifies 17α-monooxygenase), may be regulated by hypoxia.19

Hypoxia-inducible factor 1 is a heterodimeric protein consisting of an oxygen-labile HIF-1α subunit and a constitutively expressed HIF-1β subunit.20 Due to the oxygen-sensitive oxygen-dependent degradation domain (ODD) in the HIF-1α subunit, HIF-1 is rapidly degraded under normoxia.21 HIF-1α is stabilized under hypoxia and is transported into the nucleus where it dimerizes with HIF-1β22 to form a heterodimeric complex that binds to the hypoxia-responsive elements (HREs) in the promoters of HIF-target genes.23 HIF-1 has been shown to negatively regulate StAR expression and steroidogenesis in granulosa cells under hypoxia.24

Although more than 100 different genes are now known to be controlled by HIF-1, whether steroidogenic enzyme genes are regulated under hypoxia through HIF-1 remains poorly understood. Zebrafish (Danio rerio) have 2 zHIF-1α paralogs: zHIF-1aA and zHIF-1aB. Recent studies found that zHIF-1aB was more sensitive to oxygen than zHIF-1aA in zebrafish.25,26 Using a gain-of-function and loss-of-function approach, we describe here the effects of zHIF-1α (zHIF-1aB) on the expression of 9 selected steroidogenic genes and embryonic development in zebrafish (D. rerio).

Materials and Methods

Zebrafish maintenance

Wild-type adult zebrafish were maintained under a constant 14 hours:10 hours light:dark cycle at 28°C. Flow through systems were set up in the laboratory to provide a constant temperature and normoxic (7.0 ± 0.2 mg O2 L−1) or hypoxic (1.0 ± 0.2 mg O2 L−1) environments during experimental periods as described previously by Shang et al.10 Embryos were obtained by natural spawning and cultured in zebrafish E3 embryo medium (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, 0.33 mM MgSO4, and 5% methylene blue) at 28.5°C and staged according to Kimmel et al.27

Plasmid construction

The full-length zHIF-1α complementary DNA (cDNA) was amplified by reverse transcription polymerase chain reaction (RT-PCR) from zebrafish embryos and cloned into the pGEM-T Easy Vector (Promega, Madison, WI, USA) to form pGEMT-zHIF-1α. Due to the short half-life of the wild-type HIF-1α protein under normoxia,28 a zHIF-1α mutant construct (pGEMT-zHIF-1α-ΔODD) was synthesized to increase the stability of the zHIF-1α protein during its overexpression in normoxic embryos. To form pGEMT-zHIF-1α-ΔODD, 2 sites containing proline residues (proline-414 and proline-557) in wild-type ODD (at nucleotide positions 1234–1251: TTG GCA CCT GCA GGC and at nucleotide positions 1659–1680: ATG CTG GCT CCT TAC ATC CCA) were deleted from pGEMT-zHIF-1α using the Gene Tailor Site-Directed Mutagenesis Kit (Invitrogen, Carlsbad, CA USA).

To monitor the overexpression of zHIF-1α in developing zebrafish embryos, a zHIF-1α-eGFP fusion expression vector was constructed. The open reading frame (ORF) of enhanced green fluorescent protein (eGFP) was amplified from the pEGFP-N1 (Clontech Laboratories, Mountain View, CA, USA) plasmid using the gene-specific primers EGFP-F-Sma: 5′-TCC CCC GGG GGA ATG GTG AGC AAG GGC GA-3′ and EGFP-R-Not: 5′-TTG CGG CCG CAA TTA CTT GTA CAG CTC-3′, which contained built-in SmaI and NotI sites (underlined), respectively. SmaI-digested and NotI-digested eGFP PCR products and the pCMV-TNT vector (Promega) were purified and ligated overnight at 16°C. The plasmid for which the sequence was verified, pCMV-eGFP, was subsequently used in the subcloning experiment. The zHIF-1α-ΔODD ORF was amplified from pGEMT-zHIF-1α-ΔODD with the specific primers HIF-1α-F-Kpn: 5′-CGG GGT ACC CCG ATG GAT ACT GGA GTT GTC AC-3′ and HIF-1α-R-Sal: 5′-ACG CGT CGA CGT CGG TTG ACT TGG TCC AGA GCA C-3′, which contained built-in KpnI and SalI sites, respectively. The reverse primer was designed to eliminate the stop codon by converting TGA (stop codon) to CGA (arginine). The KpnI-digested and SalI-digested pCMV-eGFP plasmid and zHIF-1α-ΔODD PCR products were ligated with T4 ligase (Invitrogen) to produce pCMV-zHIF-1α-ΔODD-eGFP. The sequence of this construct was verified by sequencing the DNA.

Hypoxia exposure experiment

Zebrafish embryos were divided into 2 groups: 1 group was reared in a hypoxic system and the other in a normoxic system. Tanks used for hypoxia exposure were set up as previously described.10 Dissolved O2 concentrations were monitored continuously using dissolved oxygen meters (Cole-Parmer 01972-00, Vernon Hills, IL, USA) and polarographic probes (Cole-Parmer 5643-00), and adjustment was made using a dissolved oxygen controller (Cole-Parmer 01972-00). Fertilized embryos (4 hours postfertilization [hpf ]) were placed in net cages and allowed to develop to desired stages of development (24, 48, or 72 hpf). Each treatment consisted of 5 replicates of 60 embryos each. At the end of the experiments, the embryos were immersed in 200 µL of TRI Reagent (Molecular Research Center, Cincinnati, OH, USA) and stored at −80°C until RNA was extracted. The embryos used for hormone analysis were snap-frozen in liquid nitrogen and stored at −80°C. Animal care and experiments were conducted in accordance with the City University of Hong Kong animal care guidelines.

RNA isolation and first-strand cDNA synthesis

Total RNA was extracted using TRI Reagent (Molecular Research Center) according to the manufacturer’s instructions. First-strand cDNA was synthesized in a 25-µL reaction that contained 0.25 µg (0.5 µg/µL) of random primer, 5 µL of 5 × reaction buffer, 1.25 µL of 10 mM deoxynucleotide triphosphates, 0.65 µL of 40 U ribonuclease (RNase inhibitor) (Invitrogen), 1 µL of 200 U M-MLV Reverse Transcriptase (Promega), and 1 µg of total RNA, which was pretreated with RNase-free deoxyribonuclease I (Invitrogen) to eliminate genomic DNA contamination. cDNA synthesis was performed at 37°C for 1 hour. The cDNAs were kept at −20°C until they were analyzed by real-time PCR.

Real-time PCR

Real-time PCR, which was used to quantify gene expression, was performed using a StepOnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) as described previously.11 PCR assays were conducted using the SYBR Green–based detection method (#4367659; Applied Biosystems, Warrington, UK) according to the manufacturer’s instructions. Sequences of the primers used for real-time PCR are shown in Table 1. Melting curve analysis was performed after the PCR was completed to assess the amplification specificity. The identity of the PCR amplicons was confirmed by DNA sequencing. The cycling conditions were as follows: initial denaturation at 95°C for 10 minutes, followed by 40 cycles performed at 95°C for 15 seconds, and 0°C for 30 seconds. β-Actin was employed as an endogenous control for normalization. All PCR reactions were performed in triplicate. Fold change was calculated according to the following formula: fold , where ΔCT = CT(target) − CT(β-actin) and Δ(ΔCT) = ΔCT(stimulated) − ΔCT(control). All data were expressed as the mean ± SD and P ⩽ .05 was considered statistically significant.

Table 1

Sequences of real-time PCR primers.

Estradiol (E2) and testosterone (T) measurement

In total, 60 zebrafish embryos from each replicate were collected and immediately stored at −80°C for hormone extraction (E2 and T) as described by Yu et al.11 E2 and T levels were measured using a commercial competitive enzyme-linked immunosorbent assay (ELISA) kit (Cayman Chemical, Ann Arbor, MI, USA).

Overexpression of zHIF-1α

Capped eGFP and zHIF1α-ΔODD-eGFP messenger RNAs (mRNAs) were synthesized from linearized pCMV-eGFP and pCMV-zHIF-1α-ΔODD-eGFP, respectively, using the T7 mMESSAGE mMACHINE Kit (Ambion, Austin, TX, USA) according to the manufacturer’s instructions. mRNAs were microinjected at a concentration of 0.5 µg/µL in Hanks buffer containing 0.15% phenol red, and each embryo was injected with 2 nL mRNA at a 1 to 2 cell stage using a FemtoJet microinjection system (Eppendorf AG, Hamburg, Germany). Five replicates of 60 fertilized embryos each were used for both treatments. Injected embryos were incubated under normoxic conditions and sampled at various developmental stages for analyses.

zHIF-1α knockdown

Antisense morpholino oligonucleotides (MOs) were obtained from Gene Tools (Philomath, OR, USA). An MO targeted against the translation initiation region of zHIF-1α (zHIF-1α MO: 5′-CAG TGA CAA CTC CAG TAT CCA TTC C-3′) was used to knockdown the zHIF-1α protein. The standard control MO (5′-CCT CTT ACC TCA GTT ACA ATT TAT A-3′) was used as an injection control. The optimized morpholino amount for zHIF-1α knockdown experiments was 8 ng/embryo. Injected embryos were reared to the desired developmental stages under hypoxic or normoxic conditions as described above.

Statistical analysis

The Student t test was used to test the null hypothesis that there is no significant difference in mean gene expression between the treatment and concurrent control groups. Significance was set at P ⩽ .05. All statistical analyses were conducted using Prism 3.02 (GraphPad, San Diego, CA, USA).

Results

Hypoxia alters the expression pattern of steroidogenic enzyme genes

To confirm that HIF-1 signaling was activated in the zebrafish embryos under the experimental hypoxic conditions employed in this study, real-time PCR was performed to measure the expression of 2 HIF-1 target genes, IGFBP-1a29 and CITED-2,30 in embryos at 24, 48, and 72 hpf following exposure to hypoxia. Both the IGFBP-1a and CITED-2 mRNAs were significantly upregulated in embryos under hypoxia at all detected stages (Figure 1A), which indicated that the hypoxic conditions used in the study were appropriate.

Figure 1

Effects of hypoxia on the expression of hypoxia markers and steroidogenic genes. Effect of hypoxia on the expression of (A) the HIF-1 target genes, IGFBP-1a and CITED-2 and (B) 9 steroidogenic enzyme genes. Fertilized zebrafish embryos were exposed to normoxic (7.0 ± 0.2 mgO2L−1) or hypoxic (1.0 ± 0.2 mg O2L−1) conditions for 24, 48, or 72 hours. Gene expression was quantified using real-time PCR and normalized against β-actin mRNA. Data are presented as the mean relative fold change ± SD with respect to the normoxia control (its expression level was arbitrarily set to 1). Expression levels significantly different from the control are indicated by asterisks (t test, n ⩾ 4, *P ⩽ .05, **P ⩽ .01, ***P ⩽ .001). HIF-1 indicates hypoxia-inducible factor 1; mRNA, messenger RNA.

The mRNA expression levels for 9 selected steroidogenic genes in embryos exposed to hypoxia were differentially affected (Figure 1B). Five steroidogenic genes were significantly downregulated. StAR expression was downregulated by 2.1 ± 0.05 (P ⩽ .01) and 2.6 ± 0.2 (P ⩽ .001) fold at 24 and 72 hpf, respectively; HMGCR was reduced by 1.4 ± 0.02 (P ⩽ .01) and 7.1 ± 2.78 (P ⩽ .001) fold at 24 and 72 hpf, respectively; 17β-HSD2 was reduced by 1.8 ± 0.19 (P ⩽ .01) and 1.9 ± 0.2 (P ⩽ .01) fold at 48 and 72 hpf, respectively; CYP19a was downregulated by 1.5 ± 0.01 (P ⩽ .01) and 2.1 ± 0.02 (P ⩽ .05) fold at 48 and 72 hpf, respectively; and CYP19b was downregulated by 0.5 ± 0.01 (P ⩽ .001) fold at 72 hpf. In contrast, 3 genes were significantly upregulated in response to hypoxia. CYP11a was upregulated by 3.5 ± 0.62 (P ⩽ .01), 3.6 ± 0.8 (P ⩽ .01), and 1.3 ± 0.19 (P ⩽ .05) fold at 24, 48, and 72 hpf, respectively; CYP11b2 was upregulated by 1.6 ± 0.11 (P ⩽ .05) and 1.4 ± 0.09 (P ⩽ .01) fold at 24 and 48 hpf, respectively; and 3β-HSD was increased by 1.6 ± 0.11 (P ⩽ .05) and 5.4 ± 1.04 (P ⩽ .001) fold at 48 and 72 hpf, respectively. CYP17a1 mRNA expression showed no significant changes in hypoxic embryos in any of the 3 developmental stages.

Overexpression of zHIF-1α mRNA delays embryonic development and the expression pattern of steroidogenic genes

To investigate the relationship between zHIF-1α expression and steroidogenesis in zebrafish, overexpression and knockdown of zHIF-1α were conducted by microinjection of the in vitro–transcribed zHIF-1α-ΔODD-eGFP mRNAs and zHIF-1α translation-blocking morpholinos (MO) into embryos, respectively. Five treatment replicates were used in each experiment. The total numbers of embryos injected with the eGFP mRNA (control) and zHIF-1α-ΔODD-eGFP mRNA were 719 and 1081, respectively. Based on the embryo staging criteria described by Westerfield,31 the development of embryos microinjected with zHIF-1α-ΔODD-eGFP mRNA was found to be delayed by at least 6 hours when compared with that of the eGFP-injected and uninjected controls at 24 hpf (Figure 2A).

Figure 2

Effects of zHIF-1α overexpression on embryonic development and the expression of hypoxia markers and steroidogenic genes. (A) Development of zebrafish embryos at 24 hpf. a, uninjected control embryo; b, embryo microinjected with eGFP mRNA shows normal development (Prim-6 stage); c, embryo microinjected with zHIF-1α-ΔODD-eGFP mRNA shows developmental retardation (at approximately 18-somite stage). (B) Percentage of mortality and morphological abnormality in zebrafish embryos (24 hpf) microinjected with zHIF-1α mRNA. The effect of zHIF-1α overexpression on steroidogenic gene expression was examined by microinjecting zHIF-1α mRNA into 1- to 2-cell stage zebrafish embryos and maintained under normoxic conditions. Zebrafish embryos microinjected with eGFP mRNA and maintained under normoxia were used as the control. Because no statistically significant difference was observed between the data sets of the microinjected control and the no-injection control (One-way analysis of variance with a P < .05 threshold), the latter data set is presented here. Control: uninjected embryos (n = 373, dead = 31, abnormal = 4); eGFP: embryos microinjected with eGFP mRNA (n = 719, dead = 48, abnormal = 16); and zHIF-1α: embryos microinjected with zHIF-1α-ΔODD-eGFP mRNA (n = 1081, dead = 193, abnormal = 297). (C) Effects of zHIF-1α overexpression on IGFBP-1a and CITED-2 expression. (D) Effects of zHIF-1α overexpression on steroidogenic gene expression. Gene expression was quantified using real-time PCR and normalized with β-actin mRNA. Data are presented as the mean relative fold change ± SD with respect to the gene expression level in the control (its expression level was arbitrarily set to 1) for each experiment. Expression levels significantly different from the control are indicated by asterisks (t test, n ⩾ 4, *P ⩽ .05, **P ⩽ .01, ***P ⩽ .001). eGFP indicates enhanced green fluorescent protein; mRNA, messenger RNA.

At 24 hpf, 17.8% of the embryos microinjected with zHIF-1α-ΔODD-eGFP mRNA were dead compared with the 6.1% dead embryos observed with eGFP mRNA microinjection. The percentage of mortality of the embryos microinjected with zHIF-1α mRNAs was approximately 3-fold greater than that of the eGFP control (Figure 2B). In addition, the percentage of abnormal embryos observed in the zHIF-1α microinjection experiments was 27.5% (297/1081), whereas only 2.2% (16/719) of abnormal embryos were observed in the eGFP control group. Following microinjection with the zHIF-1α-ΔODD-eGFP mRNAs, IGFBP-1a and CITED-2 were significantly upregulated in embryos at 24 and 48 hpf (Figure 2C), which suggested successful overexpression of zHIF-1α in the embryos up until 48 hpf.

Following microinjection of the zHIF-1α mRNAs, the expression patterns of all 9 steroidogenic genes were found to be affected in zebrafish embryos at 24 and 48 hpf (Figure 2D). With the exception of StAR which was downregulated by 1.82 ± 0.12 (P ⩽ .01) and 1.43 ± 0.12 (P ⩽ .01) fold at 24 and 48 hpf, respectively, all other steroidogenic genes were upregulated by the overexpression of zHIF-1α. CYP11a was upregulated by 2.4 ± 0.04 (P ⩽ .01) and 7.0 ± 1.2 (P ⩽ .01) fold at 24 and 48 hpf, respectively; CYP11b2 was upregulated by 2.1 ± 0.6 (P ⩽ .01) and 3.4 ± 0.13 (P ⩽ .001) fold at 24 and 48 hpf, respectively; CYP17a1 was upregulated by 2.0 ± 0.3 (P ⩽ .01) and 3.2 ± 0.06 (P ⩽ .001) fold at 24 and 48 hpf, respectively; 17β-HSD2 was upregulated by 2.98 ± 0.1 (P ⩽ .05) and 1.6 ± 0.39 (P ⩽ .05) fold at 24 and 48 hpf, respectively; CYP19a was upregulated by 1.34 ± 0.05 (P ⩽ .05) and 2.8 ± 0.27 (P ⩽ .01) fold at 24 and 48 hpf, respectively; CYP19b was upregulated by 1.6 ± 0.11 (P ⩽ .01) and 4.9 ± 0.21 (P ⩽ .001) fold at 24 and 48 hpf, respectively; and 3β-HSD and HMGCR were upregulated by 1.9 ± 0.4 (P ⩽ .05) and 3.5 ± 0.4 (P ⩽ .001) fold, respectively, at 48 hpf.

Knockdown of zHIF-1α alters expression of steroidogenic genes

The experiment in which zHIF-1α was knocked down with a morpholino was conducted to complement the study in which zHIF-1α was overexpressed. Following zHIF-1α knockdown, morphological abnormalities were observed in morphant embryos at 48 hpf under normoxia (Figure 3A). The observed deformities included small head circumference, bent notochords, deformed hearts, distorted abdomens, and short curved tails. Knocking down zHIF-1α also resulted in significantly shorter zebrafish embryo body length at the same developmental stage (Figure 3A). In morphant embryos at 24 hpf, the expression levels of IGFBP-1a and CITED-2 were significantly reduced by 2.7 ± 0.48 (P ⩽ .001) and 1.9 ± 0.14 (P ⩽ .01) fold, respectively, compared with those of the control. However, IGFBP-1a was elevated by 1.6 ± 0.18 (P ⩽ .001) fold in morphants at 48 hpf, whereas CITED-2 expression was unaffected.

Figure 3

Effects of zHIF-1α knockdown on zebrafish embryonic development and the expression of hypoxia markers and steroidogenic genes under normoxia. (A) Morphological abnormalities of the zHIF-1α knockdown embryos at 48 hpf under normoxia. The effect of zHIF-1α knockdown on steroidogenic gene expression was examined by microinjecting zHIF-1α MOs into 1- to 2-cell stage zebrafish embryos and maintained under either normoxia or hypoxia. Zebrafish embryos microinjected with a standard control MO were used as controls. All morphants were imaged at the left lateral view. a, Control: embryo injected with control MO; b-j, embryos microinjected with zHIF-1α MO; b, c, and d, morphants with decreased body length, bended notochords (black arrows), and distorted abdomens (blue arrows); e-j, morphants with bended notochords, small head circles (red circles), and short curved tail (yellow arrows). (B) Effects of zHIF-1α knockdown on IGFBP-1a and CITED-2 expression under normoxia. (C) Effects of zHIF-1α knockdown on steroidogenic gene expression under normoxia. Gene expression was quantified by real-time PCR and normalized against β-actin mRNA. Data are presented as the mean relative fold change ± SD with respect to the gene expression level in the control (its expression level was arbitrarily set to 1) for each experiment. Expression levels significantly different from the control are indicated by asterisks (t test, n ⩾ 4, *P ⩽ .05, **P ⩽ .01, ***P ⩽ .001). mRNA indicates messenger RNA; MO, morpholino oligonucleotide.

The expression patterns of the steroidogenic genes in the zHIF-1α morphants were different from those of the control embryos (microinjected with control MO) under normoxia (Figure 3C). Three steroidogenic genes were specifically upregulated in the zHIF-1α morphants. StAR was significantly increased by 1.7 ± 0.24 (P ⩽ .01) fold, CYP11b2 was upregulated by 2.1 ± 0.5 (P ⩽ .05) fold, and CYP17a1 was upregulated by 1.3 ± 0.04 (P ⩽ .05) fold at 24 hpf. In contrast, CYP11a was downregulated by 4.8 ± 0.01 (P ⩽ .01) at 24 hpf, whereas 3β-HSD and 17β-HSD2 were downregulated by 1.7 ± 0.02 (P ⩽ .01) and 1.8 ± 0.1 (P ⩽ .05) fold, respectively, at 24 hpf. No significant changes in HMGCR or CYP19b expression were observed in embryos in the zHIF-1α morphants at 24 hpf compared with embryos microinjected with control MO. Among the 9 steroidogenic genes examined, only StAR, CYP11a, 3β-HSD, and 17β-HSD2 in the zHIF-1α morphants exhibited opposite patterns of expression to those of the embryos overexpressing zHIF-1α. These observations suggest that expressions of these genes might be directly and/or indirectly regulated by HIF-1.

To confirm whether changes in the expression of steroidogenic genes under hypoxia (Figure 1B) were due to regulation by zHIF-1α, embryos in which zHIF-1α had been knocked down were exposed to hypoxia. Comparing embryos in which zHIF-1α had been knocked down with the control MO-injected embryos under hypoxia, IGFBP-1a was downregulated by 4 ± 0.04 (P ⩽ .001) and 1.47 ± 0.02 (P ⩽ .01) fold at 24 and 48 hpf, respectively, whereas CITED-2 was unaffected at both stages (Figure 4A). These observations suggest that zHIF-1α was successfully downregulated in the knockdown treatment.

Figure 4

Effects of zHIF-1α knockdown on the expression of hypoxia markers and steroidogenic genes under hypoxia. (A) Effects of zHIF-1α knockdown on the hypoxia marker genes IGFBP-1a and CITED-2 under hypoxia. Methods for zHIF-1α knockdown are detailed in Figure 3. (B) Effects of zHIF-1α knockdown on steroidogenic gene expression under hypoxia. Gene expression was quantified by real-time PCR and normalized against β-actin mRNA. Data are presented as the mean relative fold change ± SD with respect to the gene expression level in the hypoxia control (its expression level was arbitrarily set to 1) for each experiment. Expression levels significantly different from the control are indicated by asterisks (t test, n ⩾ 4, *P ⩽ .05, **P ⩽ .01, ***P ⩽ .001). mRNA indicates messenger RNA.

In the zHIF-1α morphant embryos under hypoxia, StAR was significantly downregulated by 2.8 ± 0.08 (P ⩽ .001) and 2.0 ± 0.01 (P ⩽ .001) fold at 24 and 48 hpf, respectively; CYP11a was upregulated by 3.2 ± 0.26 (P ⩽ .001) and 4.7 ± 0.86 (P ⩽ .001) fold at 24 and 48 hpf, respectively; CYP11b2 was downregulated by 1.4 ± 0.13 (P ⩽ .01) and 1.7 ± 0.09 (P ⩽ .001) fold at 24 and 48 hpf, respectively; 3β-HSD was upregulated by 2.0 ± 0.05 (P ⩽ .001) and 2.7 ± 0.34 (P ⩽ .001) fold at 24 and 48 hpf, respectively; HMGCR was downregulated by 1.8 ± 0.29 (P ⩽ .001) and 2.5 ± 0.05 (P ⩽ .001) fold at 24 and 48 hpf, respectively; and CYP19a was downregulated by 1.8 ± 0.01 (P ⩽ .001) fold at 24 hpf but was undetectable at 48 hpf; CYP19b was downregulated by 1.3 ± 0.11 (P ⩽ .05) and 1.6 ± 0.17 (P ⩽ .001) fold at 24 and 48 hpf, respectively; and CYP17a1 and 17β-HSD2 were not affected in the zHIF-1α morphants at any of the time points.

E2 and T levels

The concentrations of E2 and T in whole embryos varied under normoxic and hypoxic conditions and with or without overexpression and knockdown of zHIF-1α (Figure 5). The concentrations of E2 and T in whole embryos under hypoxia were 2.2-fold less and 1.8-fold greater than those of the normoxic controls (Figure 5A), which resulted in a 3.9-fold greater T/E2 ratio (from 3.0 ± 0.42 (normoxia) to 11.73 ± 1.28 (hypoxia) (Figure 5B). The concentrations of both E2 and T were lower in embryos overexpressing zHIF-1α under normoxic conditions and in zHIF-1α knockdown embryos under hypoxic conditions (Figure 5A). The T/E2 ratio in the zHIF-1α morphants under hypoxia was similar to that of the control embryos under normoxia. However, the T/E2 ratio in embryos overexpressing zHIF-1α under normoxia was 10.3-fold lower than those maintained under hypoxia (Figure 5B).

Figure 5

Effects of hypoxia and overexpression of zHIF-1α on E2 and T levels in zebrafish embryos. (A) Effects of hypoxia and zHIF-1α expression on E2 and T levels. Methods for hypoxic exposure and zHIF-1α overexpression and knockdown are detailed in Figures 1 to 3, respectively. Hormone levels in whole zebrafish embryos (5 replicates of 60 pooled embryos each) were measured by the ELISA method. The data are presented as the mean concentration of hormone (in picograms per gram of embryo wet weight) ± SD. The different letters above the bars indicate statistically significant differences between the indicated groups for the same hormone (P < .05, Student-Newman-Keuls). (B) Effects of hypoxia and zHIF-1α expression on the T/E2 ratios. Data are presented as the mean ± SD (P < .05, Student-Newman-Keuls). The different letters above the bars indicate statistically significant differences between the indicated groups.

Discussion

There is increasing evidence that the HIF-1 signaling system plays essential roles in embryonic development in vertebrates through the activation of genes that regulate energy metabolism and survival. Previous studies have reported the expression of HIF-1α mRNA in the brain, blood vessels, somites, notochord, retina and optic tectum in zebrafish embryos,32 and HIF-1α-deficient mouse embryos (HIF-1α−/−) show defects in the formation of blood vessels and the neural fold as well as malformations in the cardiovascular system.33,34 HIF-1 has also been demonstrated to affect the development of zebrafish embryos through the regulated expression of IGFBP-1.29 In this study, zHIF-1α knockdown using morpholinos resulted in various types of morphological abnormalities in developing zebrafish embryos (Figure 3A) which indicates that the expression of zHIF-1α is essential for the normal development of zebrafish embryos under normoxic conditions and is consistent with the results from previous studies conducted in mice.35 Recently, studies have shown that in addition to its canonical role as a transcription factor, HIF-1α has a nontranscriptional role as a signaling regulator of Notch signaling,36,37 which is involved in many cellular processes during embryonic development including somitogenesis,38 muscle tissue formation,39 blood vessel maturation,40 and heart formation.41 It is reasonable to speculate that the varied morphological abnormalities observed in zebrafish embryos could be attributed to the disruption of Notch signaling as a result of ectopic zHIF-1α expression.

There is ample evidence that hypoxia impairs the sexual and reproductive development of fish through perturbation of the HPG axis at multiple levels, including enzymes that control steroidogenesis.5,10 A qualitative comparison of the results from this study on the expression patterns of steroidogenic genes in zebrafish embryos exposed to hypoxia, embryos overexpressing zHIF-1α mRNA, and zHIF-1α knockdown embryos is provided in Table 2. Changes in the expression of genes encoding steroidogenic enzymes, including StAR, CYP11a, CYP11b2, 3β-HSD, HMGCR, CYP17a1, 17β-HSD2, CYP19a, and CYP19b, would likely affect the production of steroid hormones, such as testosterone, estradiol, aldosterone, and cortisol. Steroid hormones are crucial substances that play a number of important physiological roles. An important function of steroid hormones is coordinating the physiological and behavioral responses for specific biological purposes, such as reproduction. Mutations in StAR are known to cause a complete lack of steroidogenesis in the adrenals and gonads, which, in turn, leads to classical congenital lipoid adrenal hyperplasia in humans.42 Mutations in the CYP11a1 gene causes adrenal insufficiency in humans and rabbits,43,44 whereas CYP11a1-null mice are deficient in steroid synthesis ability.45 In zebrafish embryos in which CYP11a1 was knocked down, the embryos exhibited a shortened axis and defects in cell movement during epiboly.46 CYP11B2 catalyzes 3 sequential reactions that convert 11-deoxycorticosterone to aldosterone.47 Because the biosynthesis of aldosterone is solely catalyzed by CYP11B2,48 downregulation of CYP11b2 might result in the reduced production of aldosterone. 3β-HSD deficiency results in lower amounts of adrenal steroids including mineralocorticoids, glucocorticoids, and sex steroids.14,49 Knockdown studies of CYP17 showed a significant reduction in progesterone formation and de novo synthesis of steroids in mouse tumor Leydig cells.50 Inactivating mutations in CYP19 are rare in humans, but affected individuals are unable to synthesize endogenous estrogens, which results in virilization of both the fetus and the mother.51

Table 2

Comparison of steroidogenic gene expression patterns in hypoxic, zHIF-1α overexpressing, and zHIF-1α knockdown embryos.

The results of downregulating of StAR and CYP19a under hypoxia were consistent with the results of previous observations reported in mammalian cells52 and in the male and female gonads of zebrafish, which also correlated with significant reductions of estradiol in females and testosterone in males.53 A positive relationship between StAR expression and sex steroid production has also been reported in male rainbow trout, where sex steroid hormone secretion was increased by the elevated expression of StAR.54 Repression of StAR can lead to the decreased production of sex hormones, as StAR is a cholesterol-shuttling protein that facilitates the translocation of cholesterol from the outer to the inner mitochondrial membrane in endocrine tissues, and this represents a rate-limiting step in steroidogenesis.55 In rats under hypoxia, CYP11a expression was previously reported to be stimulated, whereas StAR expression was suppressed.56–58 Similar expression patterns of CYP11a and StAR were found in this study. Concentrations of E2 and T in zebrafish embryos were less under all treatment conditions, with the exception of T, which was upregulated under hypoxia. The findings are consistent with the reduced expression of StAR under all treatment conditions relative to that of the normoxia control. The increased concentration of T under hypoxia presumably reflects the reduced conversion of T to E2 due to the downregulation of CYP19a.

Specifically, we showed that hypoxia and zHIF-1α overexpression significantly downregulated the expression of StAR but upregulated the expression of CYP11a and 3β-HSD. The reversed expression of these 3 genes in the zHIF-1α knockdown embryos under normoxia strongly suggests that these genes are directly regulated by HIF-1. Moreover, computational analysis of the 5′ flanking regions of these 3 genes revealed several putative HREs (5′-A/GCGTG-3′). Within the first 3 kb of the 5′ flanking sequence, 4, 3, and 4 putative HREs (forward/reverse) were found in StAR, CYP11a, and 3β-HSD, respectively. These results indicate that zHIF-1α may directly regulate these genes via HRE binding in the promoter region.

CYP11b2 is the only gene that showed an opposite expression pattern in hypoxia and hypoxic zHIF-1α knockdown embryos. Considering the unaffected expression pattern of CYP11b2 under zHIF-1α overexpression and knockdown treatment under normoxia, these results imply that CYP11b2 is unlikely to be regulated by zHIF-1α directly, or the HREs residing in the CYP11b2 promoter may be occupied by other HIFs under hypoxia.

Conclusions

In this study, it is demonstrated that zHIF-1 regulates the expression of several steroidogenic enzyme genes in zebrafish embryos. Specifically, the regulated expression of StAR, CYP11a, 3β-HSD, and CYP19a likely affects the production levels of the steroid hormones, E2 and T. This study provides evidence that zHIF-1 has a role in the regulation of embryonic development and steroidogenesis in zebrafish.